#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]);

}

不知道，路过

这两行：
    h=1/float(n);   
    t=1/float(m);
改为：
    h=1/(float)n;   
    t=1/(float)m;
可以编译，但是
#define N 200
#define M 100000
......
    double u[N+1][M+1];
分配了占用约153MB内存的数组，运行成问题。

先谢谢你
但是求解的时候确实需要数组那么大的,是因为这个原因导致程序无法运行么?

你这个二维数组是有点大，用malloc试试看。

该数组在运行栈上分配，一般的可执行文件不会有那么大的栈。改为 malloc 分配吧。

看了看,不会用恩,能不能帮我把用malloc的地方转化一下

将
    double u[N+1][M+1];
 改为：
    double *u = ( double * )malloc( (N+1)*(M+1)* sizeof( double ) );
    if (u == NULL) {
        printf ("No enough memory\n"); // 分配内存失败
        return 1;
    }

main 函数结束前，应加一句：
free( u );
数组访问的写法可能要改一下。


[ 本帖最后由 top398 于 2014-5-15 14:11 编辑 ]

由于 malloc 分配的内存赋值给一个 double 指针，用下标方式存取指针指向的元素时，只能以一维方式存取，不能用 u[j][k] 的形式，因此应改为：
u[j * (M+1) + k]
的形式。

谢谢了...