麻烦各位帮我看一下,运行时说程序错误,无法运行,是不是我漏掉那个地方了
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define EPSILON 1.0E-12
#define N 200
#define M 100000
void main()
{
int j,k,m,n;
float h,t;
double u[N+1][M+1];
double f[N+1];
printf("输入空间上的剖分数n\n");
scanf("%d",&n);
printf("输入时间上的剖分数m\n");
scanf("%d",&m);
h=1/float(n);
t=1/float(m);
printf("h的值为%1.11f\n",h);
printf("t的值为%1.11f\n",t);
for(j=0;j<n+1;j++)
{
u[j][0]=sin(180*j*h);
f[j]=0;
}
for(k=1;k<m+1;k++)
{
u[0][k]=0; u[n][k]=0;
for(j=1;j<n;j++)
{
u[j][k]=(1-2*t/(h*h))*u[j][k-1]+t/(h*h)*(u[j+1][k-1]+u[j-1][k-1])+t*f[j];
f[j]=j*h*k*t;
}
}
printf("输出最后每个插值点的解\n");
for(k=1;k<m+1;k++)
for(j=1;j<n;j++)
printf("%2.10f\t",u[j][k]);
}
