Ejercicio: integral doble de una superficie

Comparación entre Simpson doble y Gauss-Legendre para una integral derivada de la semiesfera x2+y2+f2=9x^2+y^2+f^2=9.

Construir el integrando

EjemploSemiesfera en el cuadrado unidad

Sea f(x,y)=9x2y2f(x,y)=\sqrt{9-x^2-y^2} y R=[0,1]×[0,1]R=[0,1]\times[0,1]. Calcula la integral de fx2+fy2\sqrt{f_x^2+f_y^2} con Simpson n=m=8n=m=8 y Gauss-Legendre n=4n=4.

  1. Las derivadas parciales son:

    fx=x9x2y2,fy=y9x2y2\frac{\partial f}{\partial x}=-\frac{x}{\sqrt{9-x^2-y^2}},\qquad \frac{\partial f}{\partial y}=-\frac{y}{\sqrt{9-x^2-y^2}}
  2. Por tanto, el integrando se simplifica a:

    (fx)2+(fy)2=x2+y29x2y2\sqrt{\left(\frac{\partial f}{\partial x}\right)^2+\left(\frac{\partial f}{\partial y}\right)^2}=\sqrt{\frac{x^2+y^2}{9-x^2-y^2}}
  3. Con Simpson doble y n=m=8 se obtiene:

    IS=0.267814255559730I_S=0.267814255559730
  4. Con Gauss-Legendre n=4n=4, tras transformar x=u+12x=\frac{u+1}{2}, y=v+12y=\frac{v+1}{2}, se obtiene:

    IG=0.267770529696778I_G=0.267770529696778

Los dos valores son muy próximos; Gauss-Legendre usa muchos menos puntos porque elige nodos no equiespaciados con pesos óptimos.