Write a local function that will return the final value of the double integral given in $(1).$ Start with

the function definition shown below. The function will receive $6$ input arguments:

$1.$ the original function f$($x$,$y$),$

$2.$ the lower bound of the double integral in the x$-$direction $($a$),$

$3.$ the upper bound of the double integral in the x$-$direction $($b$),$

$4.$ the lower bound of the double integral in the y$-$direction $($c$),$

$5.$ the upper bound of the double integral in the y$-$direction $($d$),$

$6.$ number of sub$-$intervals $($n$).$

We assume that the number of sub$-$intervals are the same in both x and y direction.

function $[$output$]=$ intdouble$($f$,$a$,$b$,$c$,$d$,$n$)$

$\%[$OUTPUT$]=$ INTDOUBLE$($F$,$A$,$B$,$C$,$D$,$N$)$ evaluates and returns the final

value of $\%$ the double integral of any functions $($f$)$

end