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

计算机习题 1、CO2气体服从范德华方程，其中CO2气体的范德华常数a = 0.364Pa·m6·mol－2 m3·mol－，、、分别为压强、温度、摩尔体积，计算、时，CO2气体的摩尔体积。 2、作业 分别用二分法、牛顿迭代法求解题1。 1）二分法，以为判据。 2）牛顿迭代法，以为判据。 #include stdio.h#include math.hfloat 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 stdio.h #include math.h 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、例题 将微分方程转变为差分方程 微分方程为，，，采用均匀步长，将微分方程转变为差分方程。 解：∵， ∴微分方程可化成以下差分方程： 即： 其中： 即： 边界条件可转变为以下形式： ，，， ，， ，， ，