计算机在材料科学中的应用计算题,计算机在材料科学中的应用习题.doc

计算机在材料科学中的应用习题

计算机习题

1、CO2气体服从范德华方程，其中CO2气体的范德华常数a = 0.364Pa·m6·mol－2 m3·mol－，、、分别为压强、温度、摩尔体积，计算、时，CO2气体的摩尔体积。

2、作业

分别用二分法、牛顿迭代法求解题1。

1)二分法，以为判据。

2)牛顿迭代法，以为判据。

#include #include float f(float x){

return 2000000*x*x*x-2*42.7*x*x+0.364*x-0.364*0.0000427-8.314*350*x*x;}void main(){float a=0,b=1,c;c=(a+b)/2;while(fabs(f(c))>0.00001 && fabs(a-b)>0.00001){if(f(c)*f(b)<0) a=c;if(f(a)*f(c)<0) b=c;c=(a+b)/2;}printf ("%f\n",c);}

2)牛顿迭代法，以为判据。

#include

#include

float f(float x) /* 定义函数f(x) */

{ return((-3*x+4)*x-5)*x+6; }

float f1(float x) /* 定义函数f(x)的导数 */

{ return (-9*x+8)*x-5; }

void main()

{ float eps,x0,x1=1.0;

printf("input eps:\n");

scanf("%f",&eps); /* 输入容许误差 */

do

{ x0=x1; /* 准备下一次迭代的初值 */

x1=x0-f(x0)/f1(x0); /* 牛顿迭代 */

}while(fabs(x1-x0)>eps); /*当满足精度,输出近似根*/

printf("x=%f\n",x1);

}

10、例题

将微分方程转变为差分方程

微分方程为，，，采用均匀步长，将微分方程转变为差分方程。

解：∵，

∴微分方程可化成以下差分方程：

即：

其中：

即：

边界条件可转变为以下形式：

，，，

，，

，，

，