Scribd
Cargar
21 vistas6 páginas

Prac Final Integrales

Este documento presenta el cálculo numérico de integrales definidas utilizando los métodos del trapecio simple, Simpson 1/3 y Simpson 3/8. Se evalúan funciones como pi*(1+(x/2)^2)^2 y exp((-x^2)/2) en diferentes intervalos y se compara el valor numérico obtenido con el valor exacto para calcular el error.

Cargado por

Richard Matute
Derechos de autor
© Attribution Non-Commercial (BY-NC)
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
21 vistas6 páginas

Prac Final Integrales

Este documento presenta el cálculo numérico de integrales definidas utilizando los métodos del trapecio simple, Simpson 1/3 y Simpson 3/8. Se evalúan funciones como pi*(1+(x/2)^2)^2 y exp((-x^2)/2) en diferentes intervalos y se compara el valor numérico obtenido con el valor exacto para calcular el error.

Cargado por

Richard Matute
Derechos de autor
© Attribution Non-Commercial (BY-NC)
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

% 1)a

clc,clear,
a=0;b=2;
f=inline('pi*(1+(x/2)^2)^2');
I=((f(a)+f(b))/2)*(b-a);
Error=abs(11.7286-I);
fprintf('INTEGRAL= %f y el ERROR ES= %f',I,Error);

% 1)b
clc,clear,
a=0;b=2;
c=(a+b)/2;
h=(b-a)/2;
f=inline('pi*(1+(x/2)^2)^2');
I=(f(a)+4*f(c)+f(b))*(h)/3;
Error=abs(11.7286-I);
fprintf('INTEGRAL= %f y el ERROR ES= %f',I,Error);

% 1)c
clc,clear,
a=0;b=2;
h=(b-a)/3;
x1=a+h;x2=a+2*h;
f=inline('pi*(1+(x/2)^2)^2');
I=(f(a)+3*f(x1)+3*f(x2)+f(b))*(3*h/8);
Error=abs(11.7286-I);
fprintf('INTEGRAL= %f y el ERROR ES= %f',I,Error);

%2 a)
clc,clear,
f=inline('pi*(1+(x/2)^2)^2')
a=0;b=2;
s=0;
n=input('numero de divisiones=');
clc,
h=(b-a)/n;
for i=1:n-1
x(i)=a+h*i;
s=s+f(x(i));
end
I=(f(a)+f(b)+2*s)*(h/2);
Error=abs(11.7286-I);
fprintf('\nFUNCION= pi*(1+(x/2)^2)^2\n');
fprintf('\nNUMERO DE INTERVALOS= %d\n',n);
fprintf('\nAMPLITUD DEL INTERVALO= %f\n',h);
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

% 2)b
clc,clear,
f=inline('pi*(1+(x/2)^2)^2');
a=0;b=2;
s=0;

r=0;
n=input('numero de divisiones= ');
clc,
h=(b-a)/n;
for i=[Link]n-1
x(i)=a+h*i;
s=s+f(x(i));
end
for i=[Link]n-2
x(i)=a+h*i;
r=r+f(x(i));
end
I=(f(a)+f(b)+4*s+2*r)*(h/3);
Error=abs(11.7286-I);
fprintf('\nFUNCION= pi*(1+(x/2)^2)^2\n');
fprintf('\nNUMERO DE INTERVALOS= %d\n',n);
fprintf('\nAMPLITUD DEL INTERVALO= %f\n',h);
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

% 2)c
clc,clear,
f=inline('pi*(1+(x/2)^2)^2');
a=0;b=2;
s=0;
r=0;
n=input('numero de divisiones multiplo de 3=');
clc,
h=(b-a)/n;
for i=1:n-1
x(i)=a+h*i;
s=s+f(x(i));
end
for i=[Link]n
x(i)=a+h*i;
r=r+f(x(i));
end
I=(f(a)+f(b)+3*s-r)*((3*h)/8);
Error=abs(11.7286-I);
fprintf('\nFUNCION= pi*(1+(x/2)^2)^2\n');
fprintf('\nNUMERO DE INTERVALOS= %d\n',n);
fprintf('\nAMPLITUD DEL INTERVALO= %f\n',h);
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

%3 a);
clc,clear,
f=inline('exp((-x^2)/2)')
a=-0.8;b=1.5;
s=0;
n=input('numero de divisiones=');
clc,
h=(b-a)/n;
for i=1:n-1

x(i)=a+h*i;
s=s+f(x(i));
end
I=(f(a)+f(b)+2*s)*(h/2);
I=I*(1/(sqrt(2*pi)));
%VALOR EXACTO:
INT=inline(int(sym(f)));
D=INT(1.5)-INT(-0.8);
Error=abs(D-I);
fprintf('\nFUNCION= exp((-x^2)/2)\n');
fprintf('\nNUMERO DE INTERVALOS= %d\n',n);
fprintf('\nAMPLITUD DEL INTERVALO= %f\n',h);
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

% 3)b
clc,clear,
f=inline('exp((-x^2)/2)');
a=-0.8;b=1.5;
s=0;
r=0;
n=input('numero de divisiones= ');
clc,
h=(b-a)/n;
for i=[Link]n-1
x(i)=a+h*i;
s=s+f(x(i));
end
for i=[Link]n-2
x(i)=a+h*i;
r=r+f(x(i));
end
I=(f(a)+f(b)+4*s+2*r)*(h/3);
I=I*(1/(sqrt(2*pi)));
%VALOR EXACTO:
INT=inline(int(sym(f)));
D=INT(1.5)-INT(-0.8);
Error=abs(D-I);
fprintf('\nFUNCION= exp((-x^2)/2)\n');
fprintf('\nNUMERO DE INTERVALOS= %d\n',n);
fprintf('\nAMPLITUD DEL INTERVALO= %f\n',h);
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

% 3)c
clc,clear,
f=inline('exp((-x^2)/2)')
a=-0.8;b=1.5;

F=inline('exp((-((2.3*z+0.7)/2)^2)/2)');
w2=8/9;w1=5/9;w3=w1;
z1=-0.774596;z2=0;z3=-z1;
I=w1*F(z1)+w2*F(z2)+w3*F(z3);
I=I*(1/(sqrt(2*pi)));
%VALOR EXACTO:
INT=inline(int(sym(f)));
D=INT(1.5)-INT(-0.8);
Error=abs(D-I);
fprintf('\nFUNCION= exp((-x^2)/2)\n');
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

%4)INTEGRAL DOBLE
clc,clear,
F=inline('exp(x+y)');
F1=inline('exp(x)');F2=inline('exp(y)');
%integral en eje y: simpson 1/3
fprintf('\nintegral en eje y: simpson 1/3\n')
a=0;b=-4;
s=0;r=0;
n=input('numero de divisiones= ');
h=(b-a)/n;
for i=[Link]n-1
x(i)=a+h*i;
s=s+F2(x(i));
end
for i=[Link]n-2
x(i)=a+h*i;
r=r+F2(x(i));
end
Iy=(F2(a)+F2(b)+4*s+2*r)*(h/3);
%integral en eje x: cuadratura de Gauss n=3
fprintf('\nintegral en eje x: cuadratura de Gauss n=3\n')
c=1;d=3;
Fz=inline('exp(z+2)');
w2=8/9;w1=5/9;w3=w1;
z1=-0.774596;z2=0;z3=-z1;
Ix=w1*Fz(z1)+w2*Fz(z2)+w3*Fz(z3);
%SOLUCION
I=Ix*Iy;
fprintf('\n LA FUNCION ES: exp(x+y)');
fprintf('\nLA INTEGRAL DOBLE ES ES = %.10f',I);

1 ENTERO
clc,clear,
a=0;b=2;
f=inline('pi*(1+(x/2)^2)^2');
fprintf('\nLA FUNCION ES: pi*(1+(x/2)^2)^2\n\n');
I=((f(a)+f(b))/2)*(b-a);
Error=abs(11.7286-I);
fprintf('INTEGRAL TRAPECIO SIMPLE= %f y el ERROR ES= %f\n\n',I,Error);
% 1)b
c=(a+b)/2;
h=(b-a)/2;
f=inline('pi*(1+(x/2)^2)^2');
I=(f(a)+4*f(c)+f(b))*(h)/3;
Error=abs(11.7286-I);
fprintf('INTEGRAL SIMPSON 1/3= %f y el ERROR ES= %f\n\n',I,Error);
% 1)c
h=(b-a)/3;
x1=a+h;x2=a+2*h;
f=inline('pi*(1+(x/2)^2)^2');
I=(f(a)+3*f(x1)+3*f(x2)+f(b))*(3*h/8);
Error=abs(11.7286-I);
fprintf('INTEGRAL SIMPSON 3/8= %f y el ERROR ES= %f',I,Error);

2 ENTERO
%2 a)
clc,clear,
f=inline('pi*(1+(x/2)^2)^2');
a=0;b=2;
s=0;
n=input('numero de divisiones=');
clc,
h=(b-a)/n;
for i=1:n-1
x(i)=a+h*i;
s=s+f(x(i));
end
I=(f(a)+f(b)+2*s)*(h/2);
Error=abs(11.7286-I);
fprintf('\nFUNCION=(1+(x/2)^2)\n');
fprintf('\nNUMERO DE INTERVALOS= %d\n',n);
fprintf('\nAMPLITUD DEL INTERVALO= %f\n\n',h);
fprintf('\nTRAPECIO EXTENDIDO\n');
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);
% 2)b
s=0;r=0;
for i=[Link]n-1

x(i)=a+h*i;
s=s+f(x(i));
end
for i=[Link]n-2
x(i)=a+h*i;
r=r+f(x(i));
end
I=(f(a)+f(b)+4*s+2*r)*(h/3);
Error=abs(11.7286-I);
fprintf('\nSIMPSON 1/3\n');
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

% 2)c
s=0;r=0;
for i=1:n-1
x(i)=a+h*i;
s=s+f(x(i));
end
for i=[Link]n
x(i)=a+h*i;
r=r+f(x(i));
end
I=(f(a)+f(b)+3*s-r)*((3*h)/8);
Error=abs(11.7286-I);
fprintf('\nSIMPSON 3/8\n');
fprintf('\nLA INTEGRAL = %f\n',I);
fprintf('\nEL ERROR ES = %f\n',Error);

También podría gustarte