problem.c

来自「这里是有关材料力学的工程分析软件,是学习工程力学的学生的一个不可少的工具,欢迎大」· C语言 代码 · 共 36 行

#define  n   6
main ()
{int i;
float a[n],b[n],c[n],d[n],p[n],q[n],x[n];
printf("qing shu ru a[]\n");
for(i=0;i<n;i++)
 scanf("%f",&a[i]);
printf("qing shu ru b[]\n");
for(i=0;i<n;i++)
 scanf("%f",&b[i]);
printf("qing shu ru c[]\n");
for(i=0;i<0;i++)
 scanf("%f",&c[i]);
printf("qing shu ru d[]\n");
for(i=0;i<n;i++)
 scanf("%f",&d[i]);
for(i=0;i<n;i++)
 if(i==0)
   {p[0]=d[0]/b[0];
    q[0]=c[0]/b[0];
    }
  else
   {p[i]=(d[i]-a[i]*p[i-1])/(b[i]-a[i]*q[i-1]);
    q[i]=c[i]/(b[i]-a[i]*q[i]);
    }
 for(i=n;i>=0;i--)
  if(i==n)
   x[i]=(d[i]-a[i]*p[n-1])/(b[i]-a[i]*q[i-1]);
  if(i==(n-1))
   x[i]=d[i]/a[i]-b[i]*x[i+1]/a[i];
 if(i<(n-1))
  x[i]=d[i]/a[i]-b[i]*x[i+1]/a[i]-c[i]*x[i+2]/a[i];
for(i=0;i<n;i++)
 printf("x[%d]=%f\n",i+1,x[i]);
}