Sample Program
H-3
static double temp_correction(double t)
{
return (t + 273)/(T_cal + 273);
}
/* ------------------------------------------------------------------------- */
static double inverse_cal(double cppm)
{
/* returns the mV that would be generated if in absolute mode with
cppm in the sample cell
*/
double mv_old = 0, mv, c;
/* crude first guess
*/
if (water_flag)
mv = cppm/75 * 5000;
else
mv = cppm/3000 * 5000;
while (fabs(mv - mv_old) > 0.5) {
mv_old = mv;
c = cal(mv_old);
mv = mv_old - (c - cppm) / (cal_slope(mv_old));
}
return mv;
}
/* ------------------------------------------------------------------------- */
static double cal_slope(double mv)
{
/* returns the slope of the cal function
*/
return A_cal + (2*B_cal + (3*C_cal + (4*D_cal + 5*E_cal*mv)*mv)*mv)*mv;
}
/* ------------------------------------------------------------------------- */
static double chi(double mfw)
{
