阿达姆斯隐式法.docVIP

作业名称：常微分方程数值解法 完成日期：2012 年 4 月27日 【小组成员： 自评成绩选择95\90\85\80\75\70\65\60\50即可 】 班级 学号 姓名 主要工作说明 自评成绩 090713 座机电话号码02 王振 编写程序，及文档制作 90 090713 座机电话号码05 麻兴斌 编写流程图 0座机电话号码713 座机电话号码17 任庆伟 编写程序 阿达姆斯隐式法 #include #include float f float x,float y return y-2*x/y; //定义函数 int main printf 阿达姆斯隐式公式\n ; float h 0.05; //定义步长h float y[21],x[21],l[21];//定义数组 float k0,k1,k2,k3; int i,j; y[0] 1; for i 1;i 3;i++ //用四阶龙格-库塔方法计算y1,y2,y3 x[i-1] i-1 *h; k0 f x[i-1],y[i-1] ; k1 f x[i-1]+h/2,y[i-1]+h/2*k0 ; k2 f x[i-1]+h/2,y[i-1]+h/2*k1 ; k3 f x[i-1]+h,y[i-1]+h*k2 ; y[i] y[i-1]+h/6* k0+2*