Scribd
Upload
23 views

Integración y Diferenciación Númerica: Del Centro Del Peru

This document describes the numerical integration of various functions using the trapezoidal rule and Simpson's rule. It begins by calculating the integral of f(x)=2+3x from 0 to 5 using the trapezoidal rule with n=1, 2, 4 points. It then calculates multiple integrals using Simpson's rule, including integrating a function of x, y and z individually and together. It concludes by calculating a double integral from -1 to 1 of the function x^2 - 2y^2 + xy^3 using Simpson's composite rule.

Uploaded by

leslie torres rivera
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
23 views

Integración y Diferenciación Númerica: Del Centro Del Peru

This document describes the numerical integration of various functions using the trapezoidal rule and Simpson's rule. It begins by calculating the integral of f(x)=2+3x from 0 to 5 using the trapezoidal rule with n=1, 2, 4 points. It then calculates multiple integrals using Simpson's rule, including integrating a function of x, y and z individually and together. It concludes by calculating a double integral from -1 to 1 of the function x^2 - 2y^2 + xy^3 using Simpson's composite rule.

Uploaded by

leslie torres rivera
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 9

“AÑO DE LA CONSOLIDACION DEL MAR DE GRAU”

UNIVERSIDAD NACIONAL

“UN INGENIERO QUÍMICO

DEL CENTRO DEL PERÚ UNA EMPRESA”

FACULTAD DE INGENIERÍA QUÍMICA


ESCUELA ACADÉMICA DE INGENIERÌA QUÌMICA

Integració n y diferenciació n
nú merica
CÁTEDRA : METODOS NUMERICOS

CATEDRÁTICO : Ing. EUFRACIO ARIAS, Wilder

ESTUDIANTE : TORRES RIVERA, Leslie

SEMESTRE : IV

SECCIÓN : “A”

Huancayo –
Perú

clc
clear
%f(x)=2+3*x
a=0;
b=5;
n1=1;
h=(b-a)/n1;
disp(h)
disp('h=')
f1=2+3*a
f2=2+3*b;
A1=h/2*(f1+f2);
disp(A1)
disp('A1=')
n2=2;
h2=(b-a)/n2;
disp(h2)
disp('h2=')
x0=a;
xo1=x0+h2;
xo2=b;
f0=f1;
f01=2+3*xo1;
f02=f2;
A2=h2/3*(f0+4*f01+f02);
disp(A2)
disp('A2=')
n3=3
h3=(b-a)/n3
disp(h3)
disp('h3=')
xo=a;
x1=xo+h3;
x2=x1+h3;
x3=b;
f00=f1;
f001=2+3*x1;
f002=3+3*x2;
f003=f2;
A3=3*h3/8*(f00+3*f001+3*f002+f003);
disp(A3)
disp('A3=')
n4=4
h4=(b-a)/n4
disp(h4)
disp('h4=')
xo=a;
x01=xo+h4;
x02=x01+h4;
x03=x2+h4;
x04=b;
f000=f1;
f0001=2+3*x01;
f0002=3+3*x02;
f0003=2+3*x03;
f0004=f2;
A4=2*h4/45*(7*f000+32*f0001+12*f0002+32*f0003+7*f0004);
disp(A4)
disp('A4=')
n5=5;
h5=(b-a)/n5;
disp(h5)
disp('h5=')
xo=a;
x01=xo+h4;
x02=x01+h4;
x03=x2+h4;
x04=x03+h4;
x05=b;
f0000=f1;
f00001=2+3*x01;
f00002=3+3*x02;
f00003=2+3*x03;
f00004=2+3*x04;
f00005=f2;
A5=5*h5/288*(19*f0000+75*f00001+50*f00002+50*f00003+75*f00004+19*f00005);
disp(A5)
disp('A5=')
n6=6;
h6=(b-a)/n6;
disp(h6)
disp('h6=')
xo=a;
x01=xo+h4;
x02=x01+h4;
x03=x2+h4;
x04=x03+h4;
x05=x04+h4;
x06=b;
f00000=f1;
f000001=2+3*x01;
f000002=3+3*x02;
f000003=2+3*x03;
f000004=2+3*x04;
f000005=2+3*x05;
f000006=f2;
A6=h6/140*(41*f00000+216*f000001+27*f000002+272*f000003+27*f000004+216*f0
00005+41*f000006);
disp(A6)
disp('A6=')

h= 5
f1= 2
A1= 47.5000
h2=2.5000
A2=47.5000
n3 = 3
h3 =1.6667
A3=49.3750
n4 = 4
h4 = 1.2500
A4= 52.6111
h5=1
A5=61.9314
h6= 0.8333
A6= 72.2440
clc,clear memory
disp('************INTEGRACION MULTIPLE*************')
disp('************SIMSOMP 1/3 OMPUESTO *************')
T=[1200 900 600 300]
f=[165 150 130 108]
h1=(T(4)-T(1))/3
x=T(1):h1:T(4)
I1=(h1/2)*[f(1)+2*(f(2)+f(3))+f(4)]

************INTEGRACION MULTIPLE*************

************SIMSOMP 1/3 OMPUESTO *************

T = 1200 900 600 300

f = 165 150 130 108

h1 = -300

x = 1200 900 600 300

I1 = -124950
clc,clear memory
disp('************INTEGRACION MULTIPLE*************')
disp('************SIMSOMP 3/8 COMPUESTO*************')
T=[0 100 200 300 400 600 800 1000]
f=[73 67 62 55 48 40 36 35]
h1=(T(4)-T(1))/3
h2=(T(8)-T(4))/3
h3=(T(5)-T(4))/1
x=T(1):h1:T(4)
%simpson 3/8
I1=(3*h1/8)*[f(1)+3*f(2)+3*f(3)+f(4)]
I2=(3*h2/8)*[f(5)+3*f(6)+3*f(7)+f(8)]
%trapecio n=1
I3=(h3/2)*[f(4)+f(5)]
I=I1+I2+I3

************INTEGRACION MULTIPLE*************

************SIMSOMP 3/8 COMPUESTO*************


T=0 100 200 300 400 600 800 1000

f = 73 67 62 55 48 40 36 35

h1 =100

h2 =233.3333

h3 =100

x = 0 100 200 300

I1 = 1.9313e+04

I2 = 2.7213e+04

I3 = 5150

I = 51675

MÉTODO TRAPEZOIDAL
clc

disp('METODO TRAPEZOIDAL')
% integral de x
x0=3;
x1=4;
%n=1
h=x1-x0;
fx0=x0;
fx1=x1;
a1=h/2*(fx0+fx1);

%n=2
h2=(x1-x0)/2;
x20=3;
x21=x20+h2;
x22=4;
fx20=x20;
fx21=x21;
fx22=x22;
a2=h2/2*(fx20+2*fx21+fx22);
%n=4
h4=(x1-x0)/4;
x40=3;
x41=x40+h4;
x42=x41+h4;
x43=x42+h4;
x44=4;
fx40=x40;
fx41=x41;
fx42=x42;
fx43=x43;
fx44=x44;
ax4=h4/2*(fx40+2*fx41+2*fx42+2*fx43+fx44);
% integral de y
y0=1;
y1=2;
%n=1
hy=y1-y0;
fy0=y0;
fy1=y1;
ay1=hy/2*(fy0+fy1);

%n=2
hy2=(y1-y0)/2;
y20=1;
y21=y20+hy2;
y22=2;
fy20=y20;
fy21=y21;
fy22=y22;
ay2=hy2/2*(fy20+2*fy21+fy22);
%n=4
hy4=(y1-y0)/4;
y40=1;
y41=y40+hy4;
y42=y41+hy4;
y43=y42+hy4;
y44=2;
fy40=y40;
fy41=y41;
fy42=y42;
fy43=y43;
fy44=y44;
ay4=hy4/2*(fy40+2*fy41+2*fy42+2*fy43+fy44);
% integral de z
z0=2;
z1=3;
%n=1
hz=z1-z0;
fz0=z0^2;
fz1=z1^2;
az1=hz/2*(fz0+fz1);

%n=2
hz2=(z1-z0)/2;
z20=2;
z21=z20+hz2;
z22=3;
fz20=z20^2;
fz21=z21^2;
fz22=z22^2;
az2=hz2/2*(fz20+2*fz21+fz22);
%n=4
hz4=(z1-z0)/4;
z40=2;
z41=z40+hz4;
z42=z41+hz4;
z43=z42+hz4;
z44=3;
fz40=z40^2;
fz41=z41^2;
fz42=z42^2;
fz43=z43^2;
fz44=z44^2;
az4=hz4/2*(fz40+2*fz41+2*fz42+2*fz43+fz44);
I=ax4*ay4*az4

INTEGRAL MULTIPLE
Evalué la siguiente integral doble con el método de Simpson compuesto
1 2

∫ ∫ x 2−2 y 2+ xy 3 dxdy
−1 0

%evalue la siguiente integral doble


%n=2
y=0
clc
x0=0;
x1=2;
h2=(x1-x0)/2;
x20=0;
x21=x20+h2;
x22=2;
fx0=x20^2-2*y^2+x20*y^3;
fx1=x21^2-2*y^2+x21*y^3;
fx2=x22^2-2*y^2+x22*y^3;
a2=h2/3*(fx0+4*fx1+fx2);
%n=4
h4=(x1-x0)/4;
x40=0;
x41=x40+h4;
x42=x41+h4;
x43=x42+h4;
x44=2;
fx40=x40^2-2*y^2+x40*y^3;
fx41=x41^2-2*y^2+x41*y^3;
fx42=x42^2-2*y^2+x42*y^3;
fx43=x43^2-2*y^2+x43*y^3;
fx44=x44^2-2*y^2+x44*y^3;
a4=h4/3*(fx40+4*fx41+2*fx42+4*fx43+fx44);
%primera integracion
y0=-1;
y1=1;
%n=4
H=(y1-y0)/4;
y0=-1;
y1=y0+H;
y2=y1+H;
y3=y2+H;
y4=1;
fy0=x40^2-2*y0^2+x40*y0^3;
fy1=x40^2-2*y1^2+x40*y1^3;
fy2=x40^2-2*y2^2+x40*y2^3;
fy3=x40^2-2*y3^2+x40*y3^3;
fy4=x40^2-2*y4^2+x40*y4^3;
A1=H/3*(fy0+4*fy1+2*fy2+4*fy3+fy4);
%segunda integracion
y0=-1;
y1=1;
%n=4
H2=(y1-y0)/4;
y20=-1;
y21=y20+H2;
y22=y21+H2;
y23=y22+H2;
y24=1;
fy20=x41^2-2*y20^2+x41*y20^3;
fy21=x41^2-2*y21^2+x41*y21^3;
fy22=x41^2-2*y22^2+x41*y22^3;
fy23=x41^2-2*y23^2+x41*y23^3;
fy24=x41^2-2*y24^2+x41*y24^3;
A2=H2/3*(fy20+4*fy21+2*fy22+4*fy23+fy24);
%tercera integracion
y0=-1;
y1=1;
%n=4
H3=(y1-y0)/4;
y30=1;
y31=y30+H3;
y32=y31+H3;
y33=y22+H3;
y34=2;
fy30=x42^2-2*y30^2+x42*y30^3;
fy31=x42^2-2*y31^2+x42*y31^3;
fy32=x42^2-2*y32^2+x42*y32^3;
fy33=x42^2-2*y33^2+x42*y33^3;
fy34=x42^2-2*y34^2+x42*y34^3;
A3=H3/3*(fy30+4*fy31+2*fy32+4*fy33+fy34);
%cuarta integracion
y0=-1;
y1=1;
%n=4
H4=(y1-y0)/4;
y40=-1;
y41=y30+H4;
y42=y31+H4;
y43=y22+H4;
y44=1;
fy40=x43^2-2*y40^2+x43*y40^3;
fy41=x43^2-2*y41^2+x43*y41^3;
fy42=x43^2-2*y42^2+x43*y42^3;
fy43=x43^2-2*y43^2+x43*y43^3;
fy44=x43^2-2*y43^2+x43*y43^3;
A4=H4/3*(fy40+4*fy41+2*fy42+4*fy43+fy44);
%quinta integracion
y0=-1;
y1=1;
%n=4
H5=(y1-y0)/4;
y50=-1;
y51=y30+H5;
y52=y31+H5;
y53=y22+H5;
y54=1;
fy50=x44^2-2*y50^2+x44*y50^3;
fy51=x44^2-2*y51^2+x44*y51^3;
fy52=x44^2-2*y52^2+x44*y52^3;
fy53=x44^2-2*y53^2+x44*y53^3;
fy54=x44^2-2*y54^2+x44*y54^3;
A5=H5/3*(fy50+4*fy51+2*fy52+4*fy53+fy54);
%tomando en cuenta la integracion total
a4=h4/3*(fx40+4*fx41+2*fx42+4*fx43+fx44);
I=h4/3*A1+4/3*h4*A2+2/3*h4*A3+4/3*h4*A4+h4/3*A5
SALIDA:
I=
4.9653

You might also like