有哪些信誉好的足球投注网站- 1、本文档共24页，可阅读全部内容。
- 2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
《数值分析》上机习题.
《数值分析》上机习题
用Newton法求方程
X7 - 28X4 + 14=0
在（0.1，1.9）中的近似根（初始近似值取为区间端点，迭代6次或误差小于0.00001）。
#includestdio.h
#includemath.h
int main()
{ float x1,x,f1,f2;
static int count=0;
x1=0.1 ;//定义初始值
do
{
x=x1;
f1=x*(x*x*x*x*x*x-28*x*x*x)+14;
f2=7*x*x*x*x*x*x-4*28*x*x*x;//对函数f1求导
x1=x-f1/f2; count++; }
while(fabs(x1-x)=1e-5);
printf(%8.7f\n,x1); printf(%d\n,count);
return 0;}
2.已知函数值如下表：
x 1 2 3 4 5 f(x) 0 01.0986123 1.3862944 1.6094378 x 6 7 8 9 10 f(x) 1.7917595 1.9459101 2.079445 2.1972246 2.3025851 f’(x) f’(1)=1 f’(10)=0.1 试用三次样条插值求f(4.563)及f’(4.563)的近似值。
include stdio.h
#include math.h
#define N 11
main()
{double B[N+1][N+1],m,x,u[N+1],y[N+1],c[N+1],d[N+1];
double e[N+1]=
{2,0,46.5916738,8.3177664,9.6566268,10.750557,11.6754606,12.47667,13.1833476,13.8155106,14.0155106};
int i;
x=4.563;
B[0][0]=-2;
B[0][1]=-4;
B[N][N-1]=4;
B[N][N]=2;
for(i=1;iN;i++)
{B[i][i-1]=1;
B[i][i]=4;
B[i][i+1]=1;
}
u[0]=B[0][0];
y[0]=e[0];
for(i=1;iN;i++)
{m=B[i][i-1]/u[i-1];
u[i]=B[i][i]-m*B[i-1][i];
y[i]=e[i]-m*y[i-1];
}
c[N]=y[N]/u[N];
for(i=N-1;i=0;i--)
c[i]=(y[i]-B[i][i+1]*c[i+1])/u[i];
for(i=0;i12;i++)
{m=fabs(x-i);
if(m=2)
d[i]=0;
else if(m=1)
d[i]=0.5*fabs(pow(m,3))-m*m+2.0/3;
else
d[i]=(-1.0/6.0)*fabs(pow(m,3))+m*m-2*fabs(m)+4/3.0;
}
m=0;
for(i=0;i12;i++)
m=m+c[i]*d[i];
printf(f(%4.3f)=%f\n,x,m);
printf(f(4.563)=%lf\n,(c[4]-c[2])/2);
}
3.用Romberg算法求.
#include stdafx.h
#includestdio.h
#includemath.h
float f(float x)
{
float f=0.0;
f=pow(3.0,x)*pow(x,1.4)*(5*x+7)*sin(x*x);
return (f);
}
main()
{
int i=1,j,k,n=12;
float T[12],a=1.0,b=3.0,s=0.0;
T[0]=0.5*(b-a)*(f(a)+f(b));
for(j=1;jn-1;j++)
{ for(k=1;k=pow(2,j-1);k++)
s+=f(a+(2*k-1)*(b-a)/pow(2,j));
T[j]=0.5*(T[j-1]+(b-a)*s/pow(2,j-1));
s=0.0;
}
T[11]=(4*T[1]-T[0])/(float)3;
for(;fabs(T[11]-T[0])0.00001;i++)
{
T[0]=T[11];
for(j=1;jn-1-i;j++) T[j]=(pow(4,i)*T[j+1]-T[j])/(pow(4,i)-1);
T
文档评论（0）