Fortran编程求助

一颗梅子在南京

对以下算法进行fortran编程：


其中xn是一维数组，也就是说每迭代计算一步，Y数组就扩大一维，对于这样的算法要怎么编程呢？求Fortran大神帮助，谢谢！（不知道我有没有将我的问题说清楚）

lqouc

嗯，不明白你最后是要一个算式还是要一个结果。
不过不管怎么说，这看起来都不难，为啥不先动手试试呢。

Delta

这样的算法是不科学

首先，数组的维数在声明以后是不可改变的，除了reshape可以，但大小一旦确定则绝不可更改，没试过用指针可不可以实现
其次，Fortran对数组维数是有限制的，77/95/03分别多少可以去相应的Fortran标准去查

所以，修改算法，看看用allocatable数组是否可同样实现，粗略看下，应该是可以的。
尽量降低数组维数是编程基本原则

一颗梅子在南京

Delta 发表于 2013-11-15 16:31

                               
登录/注册后可看大图

这样的算法是不科学

首先，数组的维数在声明以后是不可改变的，除了reshape可以，但大小一旦确 ...

我最大要算到17维，allocate不是只能定义每一维上不知道有多少个数的数组嘛，还能定义层不知道多少维数的数组？我的设想是每算一步，将算到的多维数组转化成二维的，不知道这个是否可行

Delta

声明两个可变大小数组，交替allocate、赋值、deallocate不知是否可行

虫儿飞

一种可行方案时，每算一步就将结果输出，下一步再读进来；每一步的计算可以想办法通过子程序的迭代进行，可以只需要一维数组。你这个程序，想一下子弄个数组一起算了，不太可能，fortran中，可变数组变的只是长度，而不是维数。

光头

IMPLICIT NONE
INTEGER::lobs,mobs
REAL,DIMENSION(lobs)::ry,rx
ry(1)=f(x(1))
DO mobs=2,n-1
   ry(mobs)=f(x(i))*ry(i-1)
ENDDO
END

光头

IMPLICIT NONE
INTEGER::lobs,mobs
REAL,DIMENSION(lobs)::ry,rx
ry(1)=f(x(1))
DO mobs=2,n
   ry(mobs)=f(x(i))*ry(i-1)
ENDDO
END

光头

第一个下标搞错了