文档介绍:√(100)=10,√(121)=11,√(144)=12,用拉格朗日插值法求√(115)。
C语言代码如下
/* Note:Your choice is C IDE */
#include ""
#defin++)
scanf("%f",&a[i][j]);
scanf("%f",&b[i]);
}
printf("你输入的原始增广矩阵为:\n");
for(i=0;i<N;i++){
for(j=0;j<N;j++)
printf("%f\t",a[i][j]);
printf("%f\n",b[i]);
}
for(k=0;k<N;k++){
akk=a[k][k];
for(j=k;j<N;j++){
a[k][j]=a[k][j]/akk;
}
b[k]=b[k]/akk;
for(i=k+1;i<N;i++){
aik=a[i][k];
for(j=k;j<N;j++){
a[i][j]=a[i][j]-aik*a[k][j];
}
b[i]=b[i]-aik*b[k];
}
printf("第%d次消元:\n",k);
for(i=0;i<N;i++){
for(j=0;j<N;j++)
printf("%f\t",a[i][j]);
printf("%f\n",b[i]);
}
}
for(i=2;i>=0;i--){
s=;
for(j=2;j>i;j--)
s=s+a[i][j]*b[j];
b[i]=b[i]-s;
}
for(i=0;i<N;i++){
printf("x[%d]=%f\t",i,b[i]);
}
}
输出结果如下:
雅克比迭代求方程组10x1-x2-2x3=
-x1+10x2-2x3=
-x1-x2+5x3=
程序c语言如下:
/* Note:Your choice is C IDE */
#include ""
#define N 10
#define n 3
double func(double a,double b,double c,double x,double y)
{
double fun;
fun=a*x+b*y+c;
return fun;
}
void main()
{
double func(double,double,double,double,double);
double x1[N],x2[N],x3[N];
double ann;
double a[n][n+1]={10,-1,-2,,-1,10,-2,8