program maximoFuncion

real,parameter :: pi=3.1416
real           :: deltaX , deltaY
real           :: x,y,xm,ym,maximo,funcion
integer        :: nrIntervalosX = 10 , nrIntervalosY = 10

deltaX  = pi*0.1
deltaY  = pi*0.1
maximo  = 0
xm      = 0
ym      = 0

x = 0
do i = 1,nrIntervalosX+1
   y = 0
   do j = 1,nrIntervalosY+1
      funcion = sin(x+y)            
      if (funcion >= maximo) then
          maximo = funcion
          xm = x
          ym = y
      end if
      y = y + deltaY
   end do
   x = x + deltaX
end do

print *,"el maximo (xm,ym,f(xm,ym)) es: ",xm,ym,maximo

end program maximoFuncion