I'm trying to solve a 2 dimensional (f(x,y))
numerical integral in FORTRAN. I learned that there is a standard form for function and boundary levels determination. My function has other variables like below:
Function=f(a1,a2,a3,c0,X,Y)
real(kind=8) :: a1,a2,a3,X,Y,c0(6,6),f(3,3,3,3)
I wonder how can I declare the a1,a2,a3,c0
and f
(the function output) globally to use the standard form for numerical integration (f(x,y))
?