计算机计算方法上机.docx

用二分法求方程
#include<>
#include<>
doublef(doublex){return(x*x*x*x-3*x+1);}
main()
{
doublea=,b=;
dn;j++)
{if(maxerr<fabs(y[j]-x0[j]))maxerr=fabs(y[j]-x0[j]);
}
printf("xA%d=",N);
for(i=0;i<n;i++)printf("%f",x[i]);
printf("%f",maxerr);printf("\n");
}while();
}
Lagrange插值:课本习题五第1题(P132)
#include<>
#include<>
#definen4main()
{
intk,j;
doublexx[n]={,,,};doubleyy[n]={,,,};
doublet,y=0,x=;
printf("插值节点为:\n");
for(k=0;k<n;k++)
{printf("(%f,%f)\n",xx[k],yy[k]);
}for(k=0;k<n;k++)
{t=1;
for(j=0;j<n;j++)
{if(j!=k)
t=t*(x-xx[j])/(xx[k]-xx[j]);
}
y=y+t*yy[k];
}
printf("%f\n",y);
}
U"D^IVIyDocuments\D&&ktop\Debug\Textl.&xe"
亘亘［
骷值节点为：
<>
<,>
<,>
<,>
酉此4点得至[的多坝式仕1-1300^值为0-121406
PressanyF;eytocontinue
分别利用梯形公式、新浦生公式和科特斯公
式计算定积分。(习题七第1题(P180))
#include<>
#include<>
doublef(doublex)
{
return(exp(-x));
}
doubletrapezoid(doublea,doubleb)
return(I1=(b-a)*(f(a)+f(b))/2);
}
doublesimpson(doublea,doubleb)
{
return(I2=(b-a)*(f(a)+4*f((a+b)/2)+f(b))/6);
}
doublecotes(doublea,doubleb)
{return(I3=(b-a)*(7*f(a)+32*f(3*a/4+b/4)+12*f((a+b)/2)+32*f(a/4+3*b/4)+7*f(b))
/90);
voidmain()
{
doublea=0,b=1;
doubleI1,I2,I3;
I1=trapezoid(a,b);
I2=simpson(a,b);
I3=c