program doubleint
implicit none
integer,parameter :: N=100
real,parameter :: dx=0.01, dy=0.01
real :: f,x,y,s
integer :: i,j
s=0.0
do i=1,N
do j=1,N
x=i*dx
y=j*dy
s=s+f(x,y)*dx*dy
enddo
enddo
print*,s
end program doubleint
real function f(x,y)
implicit none
real :: x,y
f=sin(x)+cos(y)
end function f
※ 引述《CYVS (Satolbuki)》之铭言：
: 小弟我以前没接触过fortran，想请教板上前辈指导一下双重积分要怎么写，
: 自己有翻书找到梯形法跟辛普森法，但实在是看不懂，能不能举个例子
: 或哪边能够参考可以告诉我，拜托了！
: 手机发文伤眼抱歉！
:

有上色!!