목차
1. 이상기체 몰랄 볼륨 구하기
2. 실제기체 몰랄 볼륨 구하기 - 뉴틈랩슨 방법
2. 실제기체 몰랄 볼륨 구하기 - 뉴틈랩슨 방법
본문내용
- - - - - - - - - - - - - -");
do{
xi2 = xi - f(xi) / fp(xi);
error = xi2 - xi;
printf("\n %3d %6.4f %6.4f %6.4f %6.4f %6.4f ", count, xi, f(xi), fp(xi), f(xi2), error);
xi = xi2;
count++;
}while(fabs(error)>0.0001 || fabs(f(xi))>0.0001);
printf("\n\n\n\n -- - - - - - - - - - - - - - - - - - - -\n");
printf(" CO2 // T = 300K P = 10atm // Molal volume %6.4f l/mol \n",xi);
printf(" - - - - - - - - - - - - - - - - - - - - - - - - - - - \n");
}
double f(double x)
{
return (10+(3.592/x)/x)*(x-0.04267)-(0.082054*300);
}
double fp(double x)
{
return 10-(3.592/x)/x+((2*3.592*0.04267/x)/x)/x;
}
do{
xi2 = xi - f(xi) / fp(xi);
error = xi2 - xi;
printf("\n %3d %6.4f %6.4f %6.4f %6.4f %6.4f ", count, xi, f(xi), fp(xi), f(xi2), error);
xi = xi2;
count++;
}while(fabs(error)>0.0001 || fabs(f(xi))>0.0001);
printf("\n\n\n\n -- - - - - - - - - - - - - - - - - - - -\n");
printf(" CO2 // T = 300K P = 10atm // Molal volume %6.4f l/mol \n",xi);
printf(" - - - - - - - - - - - - - - - - - - - - - - - - - - - \n");
}
double f(double x)
{
return (10+(3.592/x)/x)*(x-0.04267)-(0.082054*300);
}
double fp(double x)
{
return 10-(3.592/x)/x+((2*3.592*0.04267/x)/x)/x;
}