计算实习报告怎么写？( 七 ) 「问人人知识网」

{
for(j=3;j{
x3=gexian(x1,x2);
x1=x2;
x2=x3;
}
printf("%f",gexian(x1,x2)) ;
}} float f(float x) { return (x*x-exp(x)); }
float gexian(float x1,float x2) { return (x2-(f(x2)/(f(x2)-f(x1)))*(x2-x1)); } （3）实验结果及分析
四、总结
了解和学习了二分法和牛顿迭代法的思想以及程序设计的方法，比较了迭代法和牛顿法的特点：牛顿法收敛速度较快，但对初值选取要求较高；割线法计算量少。
习题三：
3线性方程组数值解法
一、实验目的
（1）熟悉求解线性方程组的有关理论和方法；
（2）会编制列主元消去法，LU分解法，雅可比及高斯-赛德尔迭代法的程序；（3）通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。
二、实验内容
1、用列主元消去法解方程组
2、用LU分解法解方程组
三、实验步骤、程序设计、实验结果及分析
1、用列主元消去法解方程组（1）程序设计 #include #include void ColPivot(float*,int,float[]); void ColPivot(float*c,int n,float x[]) {int i,j,t,k;float p;for(i=0;ik=i;
for(j=i+1;j if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
k=j;
if(k!=i)
for(j=i;j {
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;j{
p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
for(t=i;t*(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
}}for(i=n-1;i>=0;i--){
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));} } void main() {int i;float x[4];float c[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;iprintf("x[%d]=%f\\\\n",i,x[i]); } （2）实验结果及分析
(1)题
(2)题
2、用LU分解法解方程组（1）程序设计 #include void main() {float x[4];int i;float a[4][5]={48,-24,0,-12,4,
-24,24,12,12,4,
0,6,20,2,-2,
-6,6,2,16,-2};void DirectLU(float*,int,float[]);DirectLU(a[0],4,x);for(i=0;ifor(i=r;rfor(k=0;k*(u+r*(n+1)+i)-=*(u+r*(n+1)+k)*(*(u+k*(n+1)+i));
for(i=r+1;i {
for(k=0;k *(u+i*(n+1)+r)-=*(u+i*(n+1)+k)*(*(u+k*(n+1)+r));
*(u+i*(n+1)+r)/=*(u+r*(n+1)+r);
}
}for(i=n-1;i>=0;i--){
for(r=n-1;r>=i+1;r--)
*(u+i*(n+1)+n)-=*(u+i*(n+1)+r)*x[r];
x[i]=*(u+i*(n+1)+n)/(*(u+i*(n+1)+i));} }
四、总结
掌握了用列主元消去法和LU分解法求解方程组程序编写的技巧。
习题四：
4插值法
一、实验目的
（1）熟悉拉格朗日插值法多项式和牛顿插值多项式，注意其不同点；（2）掌握三次样条插值解决一些实际问题。
二、实验内容
1、按所给数据做二次插值，并求给定点的函数值
2、按所给数据做五次插值，并求给定点的函数值
3、牛顿前插公式计算函数值
三、实验步骤、程序设计、实验结果及分析
1、二次插值（1）程序设计 #include float Lagrange(float x[],float y[],float xx,int n)