• 0回复贴,共1
#include<stdio.h>
#include<math.h>
int main()
{
int i;
double k,l,tmp;
double h=0.1,x[100],y[100],f[100],p[100];
for(i=0;(0<=i)&&(i<=9);i++)
{
x[0]=0;
x[i+1]=x[i]+h;
}
for(i=0;(0<=i)&&(i<=10);i++)
{
y[0]=1;
y[i+1]=y[i]+0.1*(y[i]-(x[i]*2)/y[i]);
}
for(i=0;(0<=i)&&(i<=10);i++)
{
f[0]=1,x[0]=0;
k=f[i]-(x[i]*2)/f[i];
l=f[i]+h*k-(x[i]+h)*2/(f[i]+h*k);
tmp=(k+l)*h/2;
f[i+1]=f[i]+tmp;
}
for(i=0;(0<=i)&&(i<=10);i++)
{
p[i]=sqrt(1+x[i]*2);
printf("%2.1f %10f %10f %10f\n",x[i],y[i],f[i],p[i]);
}
return 0;
}


来自Android客户端1楼2018-05-29 11:28回复