
386 Programming for Chemical Engineers
static double delta; /* rho / rho_c */
static double tau; /* T_c / T */
static double dm1; /* delta - 1 */
static double d2; /* (delta - 1)^2 */
static double t1; /* 1 - tau */
static double t2; /* (1 - tau)^2 */
/* Temporary variables */
static double psi, theta, Delta, dc, de, te, ex;
static double phi_0, phi_r;
static double d_phi_d_tau_0, d_phi_d_tau_r;
static double d2_phi_d_tau2_0, d2_phi_d_tau2_r;
static double d_phi_d_delta_0, d_phi_d_delta_r;
static double d2_phi_d_delta2_0, d2_phi_d_delta2_r;
static double d2_phi_d_delta_d_tau_0, d2_phi_d_delta_d_
tau_r;
static double d_Delta_d_delta;
static double ...