数值分析实验报告(同名8562).docx

文档介绍：数值分析实验报告(同名8562)
河北工业大学
《数值分析》课程实验
实 验 报 告
\n");
}
运行结果

f\n",N[0]);
printf("x= %f\n",N[1]);
printf("x= %f\n",N[2]);
}
float w(float X,int n,float x[])
{
float w=;
for (int i=0;i<n;i++)
w=w*(X-x[i]);
return w;
}
float quotient(int k,int i,float f[][M],float x[],float y[])
{
if(k==0)
f[0][i]=y[i];
else
f[k][i]=(quotient(k-1,i,f,x,y)-quotient(k-1,i-1,f,x,y))/(x[i]-x[i-k]);
return f[k][i];
}
float newton(float X,int n,float f[][M],float x[],float y[])
{
float N;
N=f[0][0];
for (int i=1;i<=n;i++)
N=N+w(X,i,x)*quotient(i,i,f,x,y);
return N;
}
运行结果
拉格朗日：
牛顿：
对算法的理解与感悟
实验三
实验要求
分别用复化梯形公式和复化辛卜生公式计算f(x)=sin(x)/x的积分，并与准确值比较判断精度。
实验目的
通过实际计算来真实体会复合梯形公式和复合辛普森公式两种数值求积方法的精确度。
算法描述
1、复合梯形公式
步骤一：确定函数的区间为【a,b】，并确定其节点的个数，算出每个小区间的步长h；
步骤二：根据确定的h,代入复合梯形公式，算出其对应的结果。
2、复合辛普森公式
步骤一：确定函数的区间为【a,b】，并确定其节点的个数，算出每个步长h；
步骤二：根据确定的h,代入复化辛普生公式，算出对应的积分结果。
源程序
#include<>
#include<>
void main()
{
int m,n,k,i;
float Tn,Sn,d1,d2,a=0,b=0,c=0;
float x[1000],y[1000],z[1000],w[1000];
{
printf("将复化梯形区间划分:");
scanf("%d",&n);
for(i=1;i<n;i++)
{
x[i]=float(*i);
y[i]=(float)(sin(x[i])/x[i]);
a=a+y[i];
}
Tn=(float)((+sin()/+2*a)/);
printf("复化梯形输出:%f\n",Tn);
printf("将复化辛卜生区间划分为:");
scanf("%d",&m);
}
for (k=1;k<m;k++)
{
z[k]=float(*k);
z[k+1]=float(*(k+1));
z[k+1/2]=float(z[k]+z[k+1])/2;
w[k]=(float)(sin(z[k])/z[k]);
w[k+1/2]=(float)(sin(z[k+1/2])/((z[k]+z[k+1])/2));
b=b+w[k];
c=c+w[k+1/2];
}
Sn=(float)((+2*b+4*c+sin()/)/);
printf("复化辛卜生输出:%f\n",Sn);
{
d1=(float)(fabs(-Tn));
d2=(float)(fabs(-Sn));
printf("复化梯形误差：%f\n",d1);
printf("复化辛卜生误差：%f\n",d2);
if(d1>d2)
printf("复化梯形求法精度低于复化辛卜生求法\n");
else if (d1==d2)
printf("复化梯形求法精度等于复化辛卜生求法\n");
else (d1<d2)