东南大学数值分析资料报告上机实验题(下).doc

《东南大学数值分析资料报告上机实验题(下).doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数值分析上机报告：学号：2013年12月22日第四章38.（上机题）3次样条插值函数（1）编制求第一型3次样条插值函数的通用程序;（2）已知汽车曲线型值点的数据如下：0123456789102.513.304.044.705.225.545.785.405.575.705.80端点条件为=0.8，=0.2。用所编制程序求车门的3次样条插值函数S(x),并打印出S(i+0.5)(i=0,1,…9)。解：（1）#include#includefloatx1[100],f1[100],f2[99],f3[98],m[10

2、0],a[100][101],x,d[100];floatc[99],e[99],h[99],u[99],w[99],y_0,y_n,arr,s;inti,j,k,n,q;voidselectprint(floaty){if((y>0)&&(y!=1))cout<<"+"<

3、t)(f1[j-1]+c[j-1]*(x-x1[j-1])+m[j-1]/2.0*(x-x1[j-1])*(x-x1[j-1])+e[j-1]*(x-x1[j-1])*(x-x1[j-1])*(x-x1[j-1]));break;}returns;}voidmain(){do{cout<<"请输入n值:";cin>>n;if((n>100)

4、

5、(n<1))cout<<"请重新输入整数(1..100):"<100)

6、

7、(n<1));cout<<"请输入Xi(i=0,1,...,"<

8、i++)cin>>x1[i];cout<>f1[i];cout<>y_0>>y_n;cout<

9、(i=0;i

10、n][n-1]=1;for(i=1;i=0;i--)m[i]=(a[i][n+1]-a[i][i+1]*m[i+1])/a[i][

11、i];for(i=0;i

12、[i]);cout<<")";}if(m[i]!=0){selectprint(