大二下学期数值计算实验报告答案C++.docVIP

PAGE

PAGE1

#includeiostream

#includecmath

#defineM100

usingnamespacestd;

intmain(){

inti,j,n,k,t;

floata[M][M],x[M],temp,sum=0;

cout输入方程的阶数:endl;

cinn;

cout输入方程矩阵:endl;

for(i=0;in;i++) //输入矩阵

for(j=0;jn+1;j++)

cina[i][j];

for(k=0;kn+1;k++){

t=k;

for(i=t+1;in;i++) //选出最大元素

if((fabs(a[t][k]))(fabs(a[i][k])))

t=i;

if(t!=k) //交换行元素

for(i=k+1;in+1;i++){

temp=a[k][i];

a[k][i]=a[t][i];

a[t][i]=temp;

}

for(i=k+1;in;i++){ //消元

for(j=0;jn+1;j++)

a[i][j]-=(a[i][k]/a[k][k])*a[k][j];

}

}

x[n-1]=a[n-1][n]/a[n-1][n-1]; //解出Xn

for(i=n-2;i=0;i--){ //解出其它解

for(j=i;jn;j++){

if(j==i)continue;

sum+=a[i][j]*x[j];

}

x[i]=(a[i][n]-sum)/a[i][i];

}

printf(最终解得:\n);

for(i=0;in;i++) //输出结果

coutx[i]=x[i]endl;

return0;

}

#includeiostream

#includestdlib.h

usingnamespacestd;

intmain()

{

intn;

do{

cout输入三对角方程组的变量的数目N:;

cinn;

}while(n1);

cout请依次输入三对角方程组每行的数据(0元素除外):\n;

doublea[100],b[100],c[100],f[100];

cin*b*c*f;

for(inti=1;in-1;++i)cin*(a+i-1)*(b+i)*(c+i)*(f+i);

cin*(a+i-1)*(b+i)*(f+i);

for(i=0;in-1;++i){

*(c+i)/=(*(b+i));

*(b+i+1)-=(*(a+i))*(*(c+i));

}

*b=(*f)/(*b);

for(i=1;in;++i)

*(b+i)=((*(f+i))-(*(a+i-1))*(*(b+i-1)))/(*(b+i));

*(a+n-1)=*(b+n-1);

for(i=n-2;i=0;--i)

*(a+i)=*(b+i)-(*(c+i))*(*(a+i+1));

cout方程的解如下:\n;

for(i=1;in;++i)coutX[i]:*(a+i-1)endl;

return0;

}

#includeiostream

#includecmath

usingnamespacestd;

voidmain()

{

inti,j,n,x;

doublem=0,p,a[100],b[100];

cout插入节点数：; //得到节点数

cinn;

cout请输入节点x值：endl; //得到节点值，即X

for(i=0;in;i++)

cina[i];

cout请输入f(x)值：endl; //得到Y

for(i=0;in;i++)

cinb[i];

cout请输入x1：endl; //取未知X

cinx;

for(i=0;in;i++){ //用拉格朗日插值法求解

p=1;

for(j=0;jn;j++){

if(i!=j)

p*=(x-a[j])/(a[i]-a[j]);

else

continue;

}

m+=p*b[i];

}

coutmendl; //输出结果的Y

}

#includeiostream

usingnamespacestd;

intmain()

{

intn;

doublex;

doublea[100],f[