La integral 'signointegraldosveces con una D de subexponente' f(x,y)dxdy siendo D(recuerdo que era el subexponente, vamos como el contorno) el círculo {(x,y) € R^2: x^2 + y^2 <=1} se puede tambien evaluar mediante 'signointegraldosveces con subexponente: [0,1]x[0,2*pi]'rho*f(rho*cos(theta),rho*sen(theta))drho*dtheta.
Hacer un programa en Matlab que calcule la integral anterior tomando como parámetros de entrada f. Compruebe su programa calculando el valor de la integral de f en cartesianas (también con Matlab) con algunas elecciones de f tomadas por usted.

function a=integraldoble(funcion)

rho=0:0.1:1;
theta=0:0.1:2*pi;
[R,T]=meshgrid(rho,theta);

x=R.*cos(T);
y=R.*sin(T);

f=@(x) funcion;

a= dblquad(@(x,y) funcion,0,1,0,2*pi)

end