《计算方法实验报告电气4班李亮亮学号:1216102043》由会员分享,可在线阅读,更多相关《计算方法实验报告电气4班李亮亮学号:1216102043(6页珍藏版)》请在金锄头文库上搜索。
1、计算方法实验报告姓名:李亮亮班级:电气4班学号:12161020431、 用追赶法解下列三对角方程组。(1)(2)解:(1)c语言:#include #define n 3main()double unn=0,lnn=0,yn=0,xn=0,an=0,bn=0,cn-1=0,s1=0,s3,s4,on-1=0,pn=0,qn=0,Ann=5,1,0,1,5,1,0,1,5,Bn=17,14,7;int i=0,j=0,k=0,m=0;for(i=0;i=0;k-)for(s4=0,j=k+1;j#define n 3main()double unn=0,lnn=0,yn=0,xn=0,an=0
2、,bn=0,cn-1=0,s1=0,s3,s4,on-1=0,pn=0,qn=0,Ann=4,-1,0,-1,4,-1,0,-1,4,Bn=2,4,10;int i=0,j=0,k=0,m=0;for(i=0;i=0;k-)for(s4=0,j=k+1;jmain()float x15=0,x25=0,x35=0;int i;for(i=1;i#define f1(x,y,z) (x)*(y)-(z)#define f2(x,y,z) (x)+(y)/(z)main()double x=0,y=1,z=2,y1=0,z1=0,h=0.1;int i;for(i=1;i=5;i+)printf(x=%f,i*h);y1=y+h*f1(i-1)*h,y,z);y=y+h/2*(f1(i-1)*h,y,z)+f1(i*h,y1,z);printf(%20f,y);z1=z+h*f2(i-1)*h,y,z);z=z+h/2*(f2(i-1)*h,y,z)+f2(i*h,y,z1);printf(%20fn,z);