c X(X+Y)DXDY
READ(*,*) A1,A2,B1,B2,N
X=A1
Y=A2
H1=(B1-A1)/N
H2=(B2-A2)/N
F1=X*(x+y)
DO 20,J=1,N
DO 10,I=1,N
S1=F1*H1
S=S1+S
X=X+H1
F1=X*(x+y)
10 CONTINUE
Y=Y+H2
S2=S*H2
SS=S2+SS
X=A1
S1=0.
S=0.
F1=X*(x+y)
20 CONTINUE
WRITE(*,*) SS
END
我初学fortran，上面是一段双重积分x*(x+y)的程序，运行结果时错误，但我发现不出
错在哪里，忘哪位高手大哥给与指点！！
谢了！
或者，有关于双重积分的程序给我也可以！！
谢了！！
