CALCULO 2

PROGRAMA: INTEGRAL

#include"conio.h"
#include"iostream"
#include"stdlib"
#include"math"
float x[15],y[15];
float A=0,B=0,C=0,D=0,E=0;
float F=0,G=0,H=0,I=0,J=0,K=0,L=0;
float M=0,N=0,O=0,P=0,solucion;
int i,j;

float f(float X,float Y){ return (4-(1/9)*(X*X)-(1/16)*(Y*Y));}
main(){system("pause");
float a,b,c,d,h,k,n=2,m=2;
 cout<<"                   Metodo de integracion multiple      nn b   d   nõ   õ  F(x,y) dxdyn a   c  ";
       cout<<"ndame el intervalo (a): ";
       cin>>a;
       cout<<"dame el intervalo (b): ";
       cin>>b;
       cout<<"dame el intervalo (c): ";
       cin>>c;
       cout<<"dame el intervalo (d): ";
 cin>>d;

h=(b-a)/(2*n);k=(d-c)/(2*m);

for(i=0;i<=(2*n);i++){
 x[i]=a+(i*h);
}
for(j=0;j<=(2*m);j++){
 y[j]=c+(j*k);
}

 A=f(x[0],y[0]);
 for(i=1;i<=(n-1);i++){
 B=B+f(x[2*i],y[0]); }
 for(i=1;i<=n;i++){
 C=C+f(x[2*i-1],y[0]);}
 D=f(x[2*n],y[0]);
 for(j=1;j<=(m-1);j++){
 E=E+f(x[0],y[2*j]);}
 for(j=1;j<=(m-1);j++){
 for(i=1;i<=(n-1);i++){
 F=F+f(x[2*i],y[2*j]);}}
 for(j=1;j<=m-1;j++){
 for(i=1;i<=n;i++){
 G=G+f(x[2*i-1],y[2*j]);}}
 for(j=1;j<=m-1;j++){
 H=H+f(x[2*n],y[2*j]);}
 for(j=1;j<=m;j++){
 I=I+f(x[0],y[2*j-1]);}
 for(j=1;j<=m;j++){
 for(i=1;i<=n-1;i++){
 J=J+f(x[2*i],y[2*j-1]);}}
 for(j=1;j<=m;j++){
 for(i=1;i<=n;i++){
 K=K+f(x[2*i-1],y[2*j-1]);}}
 for(j=1;j<=n;j++){
 L=L+f(x[2*n],y[2*j-1]);}
 M=f(x[0],y[2*m]);
 for(i=1;i<=(n-1);i++){
 N=N+f(x[2*i],y[2*m]);}
 for(i=1;i<=n;i++){
 O=O+f(x[2*i-1],y[2*m]);}
 P=f(x[2*n],y[2*m]);
solucion/*1*/=(G*8+H*2+I*4+J*8+K*16+L*4+M+N*2+O*4+P);
solucion/*2*/=(A+B*2+C*4+D+E*2+F*4+solucion);
solucion/*3*/=((h*k)/9)*(solucion);
cout<<"nõ   õ (4-(1/9)(X^2)-(1/16)(Y^2));  dxdy= "<<solucion;
getch();
}