微分6.c

使用五点差分格式计算一个微分方程在边值问题上的解

#include<math.h>
#include<stdio.h>
main()
{
 float a[14][15],b;
 int i,j,k;
 a[1][0]=0.0; a[1][1]=16.0; a[1][2]=0.0; a[1][3]=0.0; a[1][4]=0.0; a[1][5]=16.0;
 a[1][6]=-64.0; a[1][7]=16.0; a[1][8]=0.0; a[1][9]=0.0; a[1][10]=0.0; a[1][11]=16.0;
 a[1][12]=0.0; a[1][13]=0.0; a[2][0]=0.0; a[2][1]=0.0; a[2][2]=16.0; a[2][3]=0.0;
 a[2][4]=0.0; a[2][5]=0.0; a[2][6]=16.0; a[2][7]=-64.0; a[2][8]=16.0; a[2][9]=0.0;
 a[2][10]=0.0; a[2][11]=0.0; a[2][12]=16.0; a[2][13]=0.0; a[13][0]=0.0; a[13][1]=0.0;
 a[13][2]=0.0; a[13][3]=16.0; a[13][4]=0.0; a[13][5]=0.0; a[13][6]=0.0; a[13][7]=16.0;
 a[13][8]=-64.0; a[13][9]=16.0; a[13][10]=0.0; a[13][11]=0.0; a[13][12]=0.0; a[13][13]=16.0;
 a[10][0]=0.0; a[10][1]=0.0; a[10][2]=0.0; a[10][3]=0.0; a[10][4]=0.0; a[10][5]=8.0; a[10][6]=0.0;
a[10][7]=0.0; a[10][8]=0.0; a[10][9]=0.0; a[10][10]=-32.0; a[10][11]=8.0; a[10][12]=0.0; a[10][13]=0.0;
a[0][0]=8.0; a[0][1]=0.0; a[0][2]=0.0; a[0][3]=0.0; a[0][4]=0.0; a[0][5]=-32.0;  a[0][6]=8.0;
 a[0][7]=0.0; a[0][8]=0.0; a[0][9]=0.0; a[0][10]=8.0; a[0][11]=0.0; a[0][12]=0.0; a[0][13]=0.0;
 a[5][0]=-32.0; a[5][1]=8.0; a[5][2]=0.0; a[5][3]=0.0; a[5][4]=0.0; a[5][5]=8.0; a[5][6]=0.0;
 a[5][7]=0.0; a[5][8]=0.0; a[5][9]=0.0; a[5][10]=0.0; a[5][11]=0.0; a[5][12]=0.0; a[5][13]=0.0;
 a[11][0]=0.0; a[11][1]=0.0; a[11][2]=0.0; a[11][3]=0.0; a[11][4]=0.0; a[11][5]=0.0; a[11][6]=8.0;
a[11][7]=0.0; a[11][8]=0.0; a[11][9]=0.0; a[11][10]=8.0; a[11][11]=-32.0; a[11][12]=8.0; a[11][13]=0.0;
 a[6][0]=8.0; a[6][1]=-32.0; a[6][2]=8.0; a[6][3]=0.0; a[6][4]=0.0; a[6][5]=0.0 ;  a[6][6]=8.0;
 a[6][7]=0.0; a[6][8]=0.0; a[6][9]=0.0; a[6][10]=0.0; a[6][11]=0.0; a[6][12]=0.0; a[6][13]=0.0;
 a[12][0]=0.0; a[12][1]=0.0; a[12][2]=0.0; a[12][3]=0.0; a[12][4]=0.0; a[12][5]=0.0; a[12][6]=0.0;
a[12][7]=8.0; a[12][8]=0.0; a[12][9]=0.0; a[12][10]=0.0; a[12][11]=8.0; a[12][12]=-32.0; a[12][13]=8.0;
 a[7][0]=0.0; a[7][1]=8.0; a[7][2]=-32.0; a[7][3]=8.0; a[7][4]=0.0; a[7][5]=0.0;  a[7][6]=0.0;
 a[7][7]=8.0; a[7][8]=0.0; a[7][9]=0.0; a[7][10]=0.0; a[7][11]=0.0; a[7][12]=0.0; a[7][13]=0.0;
 a[8][0]=0.0; a[8][1]=0.0; a[8][2]=0.0; a[8][3]=0.0; a[8][4]=0.0; a[8][5]=0.0; a[8][6]=0.0;
 a[8][7]=0.0; a[8][8]=8.0; a[8][9]=0.0; a[8][10]=0.0; a[8][11]=0.0; a[8][12]=8.0;
 a[8][13]=-32.0;a[4][0]=0.0;a[4][1]=0.0;a[4][2]=8.0;a[4][3]=-32.0;a[4][4]=8.0;
 a[4][5]=0.0;a[4][6]=0.0;a[4][7]=0.0;a[4][8]=8.0;a[4][9]=0.0;a[4][10]=0.0;
 a[4][11]=0.0;a[4][12]=0.0;a[4][13]=0.0;a[9][0]=0.0;a[9][1]=0.0;a[9][2]=0.0;
 a[9][3]=0.0;a[9][4]=8.0;a[9][5]=0.0;a[9][6]=0.0;a[9][7]=0.0;a[9][8]=8.0;
 a[9][9]=-32.0;a[9][10]=0.0;a[9][11]=0.0;a[9][12]=0.0;a[9][13]=0.0;a[3][0]=0.0;
 a[3][1]=0.0;a[3][2]=0.0;a[3][3]=8.0;a[3][4]=-32.0;a[3][5]=0.0;a[3][6]=0.0;
 a[3][7]=0.0; a[3][8]=0.0; a[3][9]=8.0; a[3][10]=0.0; a[3][11]=0.0; a[3][12]=0.0;
 a[3][13]=0.0; a[1][14]=0.0; a[2][14]=0.0; a[13][14]=0.0; a[10][14]=6.0; a[0][14]=4.0; a[5][14]=6.0;
 a[11][14]=0.0; a[6][14]=0.0; a[12][14]=-2.0; a[7][14]=-2.0; a[8][14]=-4.0-2.0*sqrt(7);
 a[4][14]=-4.0; a[9][14]=-6.0-4.0*sqrt(3); a[3][14]=-6.0-2.0*sqrt(15);
 for(k=0;k<=12;k++)
  {
   for(i=k+1;i<=13;i++)
    {
     a[i][k]=a[i][k]/a[k][k];
     for(j=k+1;j<=14;j++)
      {
       a[i][j]=a[i][j]-a[i][k]*a[k][j];
      }
    }
  }
 for(i=13;i>=0;i--)
  {
   for(j=i+1;j<=13;j++)
    {
     b=a[i][j]*a[j][14];
    }
   a[i][14]=(a[i][14]-b)/a[i][i];
  }
  printf("计算结果为：\n"); 
 for(i=0;i<=13;i++)
  {printf("对应值：%d     ",i+1);
   printf("%f\n ",     a[i][14]);
  }
  /*system("pause");*/
}