相似理论作业代码

#include
#include
main（）
{
int i;
double dxk1l1k2l2k3l3k4l4z1z2x[100]y[100]z[100];

printf（“η         f（η）          f‘（η）
“）;

dx=0.1;
x[0]=0;
y[0]=0;
y[99]=1;
z1=1;
z2=-1;

loop:	if（y[99]-1>0）
{
z2=z[0];
}

else

z1=z[0];

z[0]=（z1+z2）/2.0;
for（i=0;i<=98;i++）
{
    k1=z[i]*dx;
	l1=（-1）*x[i]*z[i]*dx/2;
    k2=（z[i]+l1/2）*dx;
	l2=（-1）*（x[i]+dx/2）*（z[i]+l1/2）*dx/2;
    k3=（z[i]+l2/2）*dx;
	l3=（-1）*（x[i]+dx/2）*（z[i]+l2/2）*dx/2;
	k4=（z[i]+l3）*dx;
	l4=（-1）*（x[i]+dx）*（z[i]+l3）*dx/2;
	x[i+1]=x[i]+dx;
	y[i+1]=y[i]+（k1+2*k2+2*k3+k4）/6;
	z[i+1]=z[i]+（l1+2*l2+2*l3+l4）/6;
}
if （fabs（y[i]-1>0.00000001））



for（i=0;i<=99;i++）

printf（“%2g%15g%15g%15g
“x[i]y[i]z[i]）;

}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-12-28 21:11  我的
     目录           0  2012-12-28 20:14  我的2.2
     文件        1754  2012-05-28 19:41  我的2.22.2新.txt
     文件        1144  2012-05-25 23:00  我的2.2老师 2.2（C++）.txt
     目录           0  2012-12-28 20:14  我的2.2改
     文件         748  2012-06-24 18:12  我的2.2改2a.cpp
     目录           0  2012-12-28 20:14  我的2.2改Debug
     文件      188495  2012-06-24 18:14  我的2.2改Debuga.exe
     文件      185344  2012-06-24 18:14  我的2.2改Debuga.ilk
     文件        4856  2012-06-24 18:14  我的2.2改Debuga.obj
     文件      221396  2012-06-24 18:09  我的2.2改Debuga.pch
     文件      361472  2012-06-24 18:14  我的2.2改Debuga.pdb
     文件       33792  2012-06-24 18:14  我的2.2改Debugvc60.idb
     文件       45056  2012-06-24 18:14  我的2.2改Debugvc60.pdb
     文件         758  2012-06-24 18:14  我的2.2改a.cpp
     文件        3341  2012-06-24 18:14  我的2.2改a.dsp
     文件         510  2012-06-24 18:15  我的2.2改a.dsw
     文件       41984  2012-06-24 18:15  我的2.2改a.ncb
     文件       48640  2012-06-24 18:15  我的2.2改a.opt
     文件         723  2012-06-24 18:14  我的2.2改a.plg
     目录           0  2012-12-28 20:14  我的2.3完成
     目录           0  2012-12-28 20:14  我的2.3完成Debug
     文件      200788  2012-06-24 16:41  我的2.3完成Debuga.exe
     文件      198776  2012-06-24 16:41  我的2.3完成Debuga.ilk
     文件        4970  2012-06-24 16:41  我的2.3完成Debuga.obj
     文件      221396  2012-06-24 16:39  我的2.3完成Debuga.pch
     文件      443392  2012-06-24 16:41  我的2.3完成Debuga.pdb
     文件       41984  2012-06-24 16:41  我的2.3完成Debugvc60.idb
     文件       45056  2012-06-24 16:41  我的2.3完成Debugvc60.pdb
     文件         953  2012-06-24 16:41  我的2.3完成a.cpp
     文件        3341  2012-06-24 16:39  我的2.3完成a.dsp
............此处省略56个文件信息