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，想請教板上前輩指導一下雙重積分要怎麼寫，
: 自己有翻書找到梯形法跟辛普森法，但實在是看不懂，能不能舉個例子
: 或哪邊能夠參考可以告訴我，拜託了！
: 手機發文傷眼抱歉！
:

有上色!!