[求助]一段fortran程序 [复制链接]

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)的程序，运行结果时错误，但我发现不出
错在哪里，忘哪位高手大哥给与指点！！
谢了！
或者，有关于双重积分的程序给我也可以！！
谢了！！

 

这么麻烦？？？
来个简单的吧
  programm main
c      X(X+Y)DXDY
      READ(*,*) A1,A2,B1,B2,N
      X=A1
      Y=A2
      H1=(B1-A1)/N
      H2=(B2-A2)/N
    summe=0
    do i=0,N-1
     do j=0,N-1
       summe=summe+(x+i*H1)*(x+i*H1+y+j*H2)*H1*H2
     end do
    end do
    write(*,*) summe
    end

如果这样写的话会更精确！！！（只对楼主的这个函数， 其它的函数要具体对待）
    programm main
c      X(X+Y)DXDY
      READ(*,*) A1,A2,B1,B2,N
      X=A1
      Y=A2
      H1=(B1-A1)/N
      H2=(B2-A2)/N
    summe=0
    do i=1,N
     do j=1,N
       summe=summe+(x+i*H1-H1/2)*(x+i*H1-H1/2+y+j*H2-H2/2)*H1*H2
     end do
    end do
    write(*,*) summe
    end

x＝1，10
y＝1，10
    x*（x+y）dx dy
正确答案是：



5447.25

谢谢大家
都是高手阿