东南大学计算方法与实习实验报告Word格式.docx
- 文档编号:13324810
- 上传时间:2022-10-09
- 格式:DOCX
- 页数:18
- 大小:129.05KB
东南大学计算方法与实习实验报告Word格式.docx
《东南大学计算方法与实习实验报告Word格式.docx》由会员分享,可在线阅读,更多相关《东南大学计算方法与实习实验报告Word格式.docx(18页珍藏版)》请在冰豆网上搜索。
。
3)按两种顺序分别计算S1000,S10000,S30000,并指出有效位数。
●实验代码
C语言程序如下:
#include<
iostream>
cmath>
usingnamespacestd;
intmain(){
floatSn=0;
intN;
cin>
>
N;
for(floatj=2;
j<
=N;
j++){
Sn=1/(j*j-1)+Sn;
}
cout<
<
"
从小到大计算的结果为"
Sn<
endl;
Sn=0;
for(j=N;
j>
=2;
j--){
从大到小计算的结果为"
return0;
}
●运行窗口
实习题二
1、用牛顿法求下列方程的根:
1)
C语言程序代码如下:
#defineN100
#defineeps1e-6
#defineeta1e-8
floatNewton(floatf(float),floatfl(float),floatx0){
floatx1,d;
intk=0;
do{
x1=x0-f(x0)/fl(x0);
if(k++>
N||fabs(fl(x1))<
eps){cout<
发散"
break;
d=fabs(x1)<
1?
x1-x0:
(x1-x0)/x1;
x0=x1;
cout<
x="
x0<
}while(fabs(d)>
eps&
&
fabs(f(x1))>
eta);
returnx1;
floatf(floatx){
returnx+log10(x)-2;
floatfl(floatx){
return1+1/x;
voidmain(){
floatx0,y0;
x0;
y0=Newton(f,fl,x0);
方程的根为"
y0<
实习题三
1、用列主元消去法解方程组:
1)
C语言程序代码如下:
#include<
voidColPivot(float*c,intn,floatx[])
{
inti,j,t,k;
floatp;
for(i=0;
i<
=n-2;
i++){
k=i;
for(j=i+1;
=n-1;
j++)
if(fabs(*(c+j*(n+1)+i))>
(fabs(*(c+k*(n+1)+i))))k=j;
if(k!
=j)
for(j=i;
=n;
p=*(c+i*(n+1)+j);
*(c+i*(n+1)+j)=*(c+k*(n+1)+j);
*(c+k*(n+1)+j)=p;
}
for(j=i+1;
p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
for(t=i;
t<
t++)*(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
}
for(i=n-1;
i>
=0;
i--){
for(j=n-1;
=i+1;
j--)
(*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
voidColPivot(float*,int,float[]);
inti;
floatx[4];
floatc[4][5]={1,1,0,3,4,
2,1,-1,1,1,
3,-1,-1,3,-3,
-1,2,3,-1,4};
ColPivot(c[0],4,x);
for(i=0;
=3;
i++)printf("
[x%d]=%f\n"
i,x[i]);
4、编写用追赶法解三对角线性方程组的程序,并解下列方程组:
2)
,其中
A10x10=-41
1-41
...
1-4
b=-27
-15
…
-15
floatx[10];
floatc[10][11]={-4,1,0,0,0,0,0,0,0,0,-27,
1,-4,1,0,0,0,0,0,0,0,-15,
0,1,-4,1,0,0,0,0,0,0,-15,
0,0,1,-4,1,0,0,0,0,0,-15,
0,0,0,1,-4,1,0,0,0,0,-15,
0,0,0,0,1,-4,1,0,0,0,-15,
0,0,0,0,0,1,-4,1,0,0,-15,
0,0,0,0,0,0,1,-4,1,0,-15,
0,0,0,0,0,0,0,1,-4,1,-15,
0,0,0,0,0,0,0,0,1,-4,-15};
ColPivot(c[0],10,x);
=9;
实习题四
2、按下列数据
Xi
0.30
0.42
0.50
0.58
0.66
0.72
Yi
1.04403
1.08462
1.11803
1.15603
1.19817
1.23223
作5次插值,并求X1=0.46,X2=0.55,X3=0.60时的函数近似值。
cstdio>
#defineN5
voidDifference(floatx[],floaty[],intn){
float*f=newfloat[n+1];
intk,i;
for(k=1;
k<
k++){
f[0]=y[k];
k;
i++)
f[i+1]=(f[i]-y[i])/(x[k]-x[i]);
y[k]=f[k];
deletef;
return;
floata,b,c,varx=0.46,vary=0.55,varz=0.60;
floatx[N+1]={0.30,0.42,0.50,0.58,0.66,0.72};
floaty[N+1]={1.044
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 东南大学 计算方法 实习 实验 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)