计算方法与实习上机题答案Word格式文档下载.docx
- 文档编号:14084167
- 上传时间:2022-10-18
- 格式:DOCX
- 页数:36
- 大小:276.33KB
计算方法与实习上机题答案Word格式文档下载.docx
《计算方法与实习上机题答案Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《计算方法与实习上机题答案Word格式文档下载.docx(36页珍藏版)》请在冰豆网上搜索。
sum[%d]=%f\n"
}
结果:
(2)逆序计算
voidmain()
{
intn=10000;
while
(1)
{
if(n%1000==0)
printf("
sum[%d]=%-30f"
if(n<
=1)break;
n--;
}
2已知连分数
利用下面的方法计算f:
写一个程序,读入n,,计算并打印f
inti=0,n;
floata[1024],b[1024],d[1024];
pleaseinputn,n="
);
scanf("
%d"
&
n);
\npleaseinputa[1]toa[n]:
\n"
for(i=1;
i<
=n;
i++)
printf("
a[%d]="
i);
scanf("
%f"
a[i]);
\npleaseinputb[0]tob[n]:
for(i=0;
b[%d]="
b[i]);
d[n]=b[n];
for(i=n-1;
i>
=0;
i--)
d[i]=b[i]+a[i+1]/d[i+1];
\nf=%f\n"
d[0]);
3给出一个有效的算法和一个无效的算法计算积分
main()
doubley_0=(1/4.0)*log(5),y_1;
doubley_2=(1.0/55.0+1.0/11.0)/2,y_3;
intn=1,m=10;
有效算法输出结果:
y[0]=%-20f"
y_0);
y_1=1.0/(4*n)+y_0/(-4.0);
y[%d]=%-20f"
n,y_1);
=10)
break;
y_0=y_1;
if(n%3==0)
}
\n无效算法的输出结果:
y[10]=%-20f"
y_2);
y_3=1.0/n-4.0*y_2;
m-1,y_3);
if(m<
=1)
y_2=y_3;
m--;
if(m%2==0)printf("
4设,已知其精确值为
(1)编制按从小到大顺序计算的程序
(2)编制按从小达到的顺序计算的程序
(3)按两种顺序分别计算,并指出有效位数
intN;
doubleSN[30000];
SN[30000]=(3.0/2.0-1.0/30000.0-1/30001.0)/2.0;
for(N=30000;
N>
=2;
N--)
SN[N-1]=SN[N]-1.0/(N*N-1);
printf("
从大到小顺序计算:
\nSN[1000]=%f\nSN[10000]=%f\nSN[30000]=%f\n"
SN[1000],SN[10000],SN[30000]);
SN[2]=(3.0/2-1.0/2.0-1/3.0)/2.0;
for(N=3;
N<
=30000;
N++)
SN[N]=SN[N-1]+1.0/(N*N-1);
从小到大顺序计算:
实习题2
1.用牛顿法求下列方程的根
#include<
typedeffloat(*p)(float);
floatff1(floatx)
returnx*x-exp(x);
floatff2(floatx)
returnx*exp(x)-1;
floatff3(floatx)
returnlog(x)+x-2;
floatanswer(float(*p)(float))
intk=2;
floatm=1,n=-1,x2,a,b,c;
if(p==ff3)n=2;
x[0]=%.4f,x[1]=%.4f,"
m,n);
while
(1)
if(fabs(m-n)<
1e-4)break;
a=p(n)*(n-m);
b=p(n)-p(m);
c=a/b;
x2=n-c;
m=n;
n=x2;
x[%d]=%.4f,"
k,x2);
k++;
if(k%3==0)printf("
if(k%3!
=0)printf("
iterationtimes:
%d,roots:
%.4f\n"
k-2,n);
return0;
main()
x*x-exp(x),\n"
answer(ff1);
x*exp(x)-1,\n"
answer(ff2);
lg(x)+x-2,\n"
answer(ff3);
2.编写一个割线法的程序,求解上述各方程
floatgexian(float,float);
floatf(float);
inti,j;
floatx1=2.2;
floatx2=2,x3;
i);
if(i==1)
x1);
elseif(i==2)
x2);
else
for(j=3;
j<
=i;
j++)
{
x3=gexian(x1,x2);
x1=x2;
x2=x3;
gexian(x1,x2));
}
floatf(floatx)
return(x*x-exp(x));
floatgexian(floatx1,floatx2)
return(x2-(f(x2)/(f(x2)-f(x1)))*(x2-x1));
实习题3
1用列主元消去法解下列方程组;
源程序:
voidColPivot(float*,int,float[]);
voidColPivot(float*c,intn,floatx[])
inti,j,t,k;
floatp;
=n-2;
k=i;
for(j=i+1;
=n-1;
if(fabs(*(c+j*(n+1)+i))>
(fabs(*(c+k*(n+1)+i))))
k=j;
if(k!
=i)
for(j=i;
{
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(j=n-1;
j>
=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));
voidmain()
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);
=3;
x[%d]=%f\n"
i,x[i]);
第
(1)题
第
(2)题
2、
floata[4][5]={48,-24,0,-12,4,-24,24,12,12,4,0,6,20,2,-2,-6,6,2,16,-2};
voidDirectLU(float*,int,float[]);
DirectLU(a[0],4,x);
i++)
voidDirect
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算方法 实习 上机 答案
![提示](https://static.bdocx.com/images/bang_tan.gif)