c语言作业的分析,数值分析作业(C语言编程实现

《数值分析作业(C语言编程实现》由会员分享，可在线阅读，更多相关《数值分析作业(C语言编程实现(12页珍藏版)》请在人人文库网上搜索。

1、include #include double f(double x)double ans; ans=exp(x); return ans;void main()double a=1,b=3,error=0.0001,t2020,h,c;int i,j,k,m,n;h=b-a;t00=h*(f(a)+f(b)/2;k=1;while(1)t0k=0;m=1;for(j=0;j#include double f(double t)double ans;ans=pow(cos(t),1.0/3);return ans;void main()double x=0,eslong=0.,x0;int N。

2、=20,i;printf(n近似初值 x0 = %lfn,x);for(i=0;i#include double a=0,b=1,x,y=0,h=0.1,k1,k2,k3,k4;int i,N;double f(double t,double s)double ans;ans=1+t*t;return ans;void main()N=(b-a)/h;x=a;printf(n 初值为 (x0,y0) = ( %.8f , %.8f )n,x,y);for(i=0;ivoid main()double datax4=1.2,2.9,4.6,5.8,datay10=14.84,33.71,58.3。

3、6,79.24,l3,x=1.5,y;int i,j;y=0;for(i=0;ivoid main()double datay9=11.7,14.87,21.44,31.39,44.73,61.46,81.57,105.11,131.91;int m=2,i,j,k;double p,data94,a34,datax9=1.2,2.3,3.4,4.5,5.6,6.7,7.8,8.9,10.0;for(i=0;i#include void main()double a34=1,-2,-1,3,-2,10,-3,15,-1,-2,5,10,mov,comp;int i,j,k,nrow;for(i。

4、=0;i=comp)nrow=k;comp=fabs(aki);for(j=0;j#include void main() double a37=1,-2,-1,3,-2,10,-3,15,-1,-2,5,10,error=0.,norm;int N=423,i,j,k;a04=0,a14=0,a24=0;for(i=0;i3;i+) /把a矩阵转化为b矩阵/ai6=aii;for(j=0;j3;j+)aij=-aij/ai6;ai3=ai3/ai6;aii=0;printf(化为b矩阵如下n);for(i=0;i3;i+)printf(%f %f %f %fn,ai0,ai1,ai2,ai3。

5、);for(i=1;iN;i+)for(j=0;j3;j+)aj5=0;for(k=0;k3;k+)aj5=ak4*ajk+aj5;aj5=aj5+aj3;norm=0;for(k=0;k3;k+)norm=norm+fabs(ak4-ak5);if(normerror)break;elsefor(k=0;k3;k+)ak4=ak5;if(normerror)printf(计算结果为n);for(i=0;i3;i+)printf( x%d = %.3fn,i+1,ai5);elseprintf(迭代失败n);现场考试题目1#include stdio.h#include math.hdoubl。

6、e f(double x)double ans;ans=exp(x);return(ans);void main()double a=-1,b=1,error=0.0001,m=1,h,T0,T,F;int k;h=(b-a)/2;T0=h*(f(a)+f(b);while(1)F=0;for(k=1;k=pow(2.0,m-1);k+)F=F+f(a+(2*k-1)*h);T=T0/2+h*F;if(fabs(T-T0)error)break;m+;h=h/2;T0=T;printf(积分结果为I %fn,T);题目2#include stdio.hdouble f(double t,double s)double ans;ans=1+t*t;return(ans);void main()double a=0,b=1,h=0.2,x0=0,y0=0,x,k1,k2,k3,y;int N,n;N=(b-a)/h;for(n=1;n=N;n+)x=x0+h;k1=f(x0,y0);k2=f(x0+h/2,y0+h/2*k1);k3=f(x0+h,y0-h*k1+2*h*k2);y=y0+h/6*(k1+4*k2+k3);printf(第%d次输出结果为(%.8f,%.8f)n,n,x,y);x0=x;y0=y。