呵呵，此程序是用来对matrix进行变换的，算法是勒让德变换。

integer n,np,indx(np)
        real*8 a(np,np),a2(np,np),y(np,np),prod(np,np),d
        character result*8
改成：
        integer n,np,
        integer,allocatable:: indx(:)
        real*8,allocatable:: a(:,:),a2(:,:),y(:,:),prod(:,:)
        real*8 d
        character result*8

        allocate(indx(np))
        allocate(a(np,np),a2(np,np),y(np,np),prod(np,np))

谢谢！

我稍后试试。不过我目前更改了矩阵大小，到没有出现内存分配不足的问题。
但是好像找不到要转换的矩阵。大家帮忙看一下程序。