- #1

- 17

- 0

## Main Question or Discussion Point

Hi every one!

I would like to draw a double integral function in related to R and h parameters by below M-File but It does goes wrong!!!

Is there any one to correct it for me???

thank you

syms R h;

a1 = 0;

a2 = atan(R./(R+h));

r1 = h;

r2 = sqrt(R.^2+(R+h).^2);

integrand = @(r,a)(h.*sin(a)/((r.^2).*(r.^2+h.^2-2.*r.*h.*cos(a))));

f = quad2d(integrand,r1,r2,a1,a2);

ezsurf(f,[0.001,5]);

I would like to draw a double integral function in related to R and h parameters by below M-File but It does goes wrong!!!

Is there any one to correct it for me???

thank you

syms R h;

a1 = 0;

a2 = atan(R./(R+h));

r1 = h;

r2 = sqrt(R.^2+(R+h).^2);

integrand = @(r,a)(h.*sin(a)/((r.^2).*(r.^2+h.^2-2.*r.*h.*cos(a))));

f = quad2d(integrand,r1,r2,a1,a2);

ezsurf(f,[0.001,5]);