Introduction - If you have any usage issues, please Google them yourself
Res = 1000*(1 + 0.0038623139728*Tem - 0.00000065314932626*Tem*Tem);
fprintf(fp,"%6.1f\t%18f\t",Tem,Res);
Tem1 = sqrt( A*A + 4.0*B*(Res/Res_0 -1.0) ) - A;
Tem1 = 1000 * Tem1/(2*B);
fprintf(fp," %8.4f\t",Tem1);
Tem2 = 0.2712508997 * Res - 274.4423762855;
fprintf(fp," %8.4f\n",Tem2);