數(shù)值分析作業(yè)(c語言編程實(shí)現(xiàn))

《數(shù)值分析作業(yè)(c語言編程實(shí)現(xiàn))》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。

1、1、（10分）用Romberg方法計(jì)算積分^的近似值，誤差不超過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、++；}}}積分結(jié)果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得到近似結(jié)果為XW0.8654743、（10分）用四階RK方法求初值問題的數(shù)值解。取步長A=0.1，計(jì)算結(jié)果取8位小數(shù)。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