数值分析作业(c语言编程实现)

《数值分析作业(c语言编程实现)》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、1、（10分）用Romberg方法计算积分^的近似值，误差不超过10'include//includedoubleffdoublex){doubleans;ans=exp(x);returnans;}voidmain(){doublea=l,b=3,error=0.0001,t[20][20],h,c;inti,j,k,m,n;h=b-a;t[0][0]=h*(f(a)+f(b))/2;k=l;while⑴{t[O][k]=O;m=l;for(j=0;j

2、m*2;for(i=l;i<=m;i++)t[0][k]=t[0][k]+h*f(a+(i-0.5)*h);t[0][k]=(t[0][k]+t[0][k-l])/2;for(j=l;j<=k;j++){c=l;for(n=0;n;break;}else{h=h/2;k

3、++；}}}积分结果Iw17.3672552、（10分）用简单迭代法求解方程x3-ccw.v=0，要求include#includedoublef(doublet){doubleans;ans=pow(cos(t),1.0/3);returnans;}voidmain(){doublex=0,eslong=0.000001,xO;intN=20,i;printf("近似初值xO=%lf'»;for(i=0;i

4、d=%lf,,/i+l/x);if(fabs(x-xO)

5、2=0.865473xl3=0.865474x!4=0.865474得到近似结果为XW0.8654743、（10分）用四阶RK方法求初值问题的数值解。取步长A=0.1，计算结果取8位小数。yf=1+x20includedoublea=0,b=l,x,y=0,h=0.1,kl,k2,k3,k4;inti,N;doublef(doublet,doubles){doubleans;ans=l+t*t;returnans;}voidmain(){N=

6、(b-a)/h;x=a;printf("初值为(xO,yO)=(%.8f,%.8f)",x,y);for(i=0;i

7、EhEhEhEhEhEh日h士口Ldb口吉吉吉吉吉吉吉吉£-v=<=<=<=<=<=<=<=<=<

8、>0.100000000.200000000.300000000.400000000.500000000.600000000.700000000.800000000.900000000.100333330.202666670.309000000.421333330.541666670.672000000.814333330.970666671.14300000<1.00000000.1.333333