UNIVERSIDAD MAYOR DE SAN ANDRES

FACULTAD DE INGENIERÍA
INGENIERÍA QUÍMICA, AMBIENTAL, ALIMENTOS Y PETROQUIMICA

PROGRAMACION Y ANALISIS NUMERICO (PRQ – 404)

DOCENTE: M. Sc. Ing. Roberto Parra Zeballos

AUXILIAR: Univ. Álvaro Abel Rivera Cárdenas
ESTUDIANTE: Univ. Flores Santos Ayrton Tadeo
GRUPO: A
GESTIÓN: 2022

TAREA No. 10
1.- Usando Excel, resuelva la ecuación que se presenta a continuación, de t = 0 a 3, con h = 0.1, con
los métodos de a) Euler modificado o Heun, y b) RK de segundo c) RK de cuarto orden. Realice
en un solo gráficos las tres soluciones.

a) Euler modificado o Heun

t0: 0 tf: 3
y0: 1
h: 0,1 f(t,y)=y*(sen(t))^3

n tn yn (yn+1)* tn+1 yn+1

0 0 1 1 0,1 1,00004975
1 0,1 1,00004975 1,00014926 0,2 1,00049163
2 0,2 1,00049163 1,00127615 0,3 1,00217596
3 0,3 1,00217596 1,00476242 0,4 1,00643595
4 0,4 1,00643595 1,01237936 0,5 1,01498563
5 0,5 1,01498563 1,02617031 0,6 1,02981453
6 0,6 1,02981453 1,04835324 0,7 1,05309832
7 0,7 1,05309832 1,08125405 0,8 1,0871335
8 0,8 1,0871335 1,12726517 0,9 1,13429035
9 0,9 1,13429035 1,18881004 1 1,19696623
10 1 1,19696623 1,26828425 1,1 1,27751248
11 1,1 1,27751248 1,36794017 1,2 1,37810461
12 1,2 1,37810461 1,48968415 1,3 1,50052871
13 1,3 1,50052871 1,63476753 1,4 1,64587021
14 1,4 1,64587021 1,80337691 1,5 1,81411648
15 1,5 1,81411648 1,99416822 1,6 2,00372327
16 1,6 2,00372327 2,20383939 1,7 2,21124079
 17 1,7 2,21124079 2,4268815 1,8 2,43113172
18 1,8 2,43113172 2,65566542 1,9 2,65591864
19 1,9 2,65591864 2,88098025 2 2,87674937
20 2 2,87674937 3,09303114 2,1 3,08436269
21 2,1 3,08436269 3,28274999 2,2 3,27030087
22 2,2 3,27030087 3,443132 2,3 3,42810444
23 2,3 3,42810444 3,57025729 2,4 3,55419503
24 2,4 3,55419503 3,66372836 2,5 3,64822846
25 2,5 3,64822846 3,72642974 2,6 3,71285332
26 2,6 3,71285332 3,76371576 2,7 3,75297479
27 2,7 3,75297479 3,78227143 2,8 3,77473215
28 2,8 3,77473215 3,78892189 2,9 3,78442142
29 2,9 3,78442142 3,78960407 3 3,78754526
30 3 3,78754526 3,78860971 3,1 3,7880911

b) RK de segundo orden
METODO RK-2

t0: 0 tf: 3
y0: 1
h: 0,1 dy/dt=f(t,y)=y*(sen(t))^3

t y k1 k2 Δy yn+
1
0 1 0 9,95011E-05 4,9751E-05 1,00004975
0,1 1,00004975 9,9506E-05 0,000784255 0,00044188 1,00049163
0,2 1,00049163 0,000784523 0,002584136 0,00168433 1,00217596
0,3 1,00217596 0,002586459 0,005933523 0,00425999 1,00643595
0,4 1,00643595 0,005943405 0,011155956 0,00854968 1,01498563
0,5 1,01498563 0,011184676 0,018473112 0,01482889 1,02981453
0,6 1,02981453 0,018538716 0,028028869 0,02328379 1,05309832
0,7 1,05309832 0,028155734 0,039914631 0,03403518 1,0871335
0,8 1,0871335 0,040131672 0,054182024 0,04715685 1,13429035
0,9 1,13429035 0,05451969 0,070832064 0,06267588 1,19696623
1 1,19696623 0,071318029 0,089774479 0,08054625 1,27751248
1,1 1,27751248 0,090427692 0,110756565 0,10059213 1,37810461
1,2 1,37810461 0,111579539 0,13326866 0,1224241 1,50052871
1,3 1,50052871 0,134238826 0,156444186 0,14534151 1,64587021
1,4 1,64587021 0,157506692 0,178985838 0,16824626 1,81411648
1,5 1,81411648 0,180051744 0,199161839 0,18960679 2,00372327
1,6 2,00372327 0,200116122 0,21491892 0,20751752 2,21124079
1,7 2,21124079 0,215640706 0,224141156 0,21989093 2,43113172
1,8 2,43113172 0,224533697 0,225040145 0,22478692 2,65591864
1,9 2,65591864 0,225061603 0,216599858 0,22083073 2,87674937
 2 2,87674937 0,216281769 0,198944861 0,20761332 3,08436269
2,1 3,08436269 0,198387303 0,173489057 0,18593818 3,27030087
2,2 3,27030087 0,172831137 0,142776 0,15780357 3,42810444
2,3 3,42810444 0,142152853 0,110028339 0,1260906 3,55419503
2,4 3,55419503 0,109533332 0,078533525 0,09403343 3,64822846
2,5 3,64822846 0,078201278 0,051048433 0,06462486 3,71285332
2,6 3,71285332 0,050862449 0,02938049 0,04012147 3,75297479
2,7 3,75297479 0,029296643 0,014218082 0,02175736 3,77473215
2,8 3,77473215 0,014189741 0,005188811 0,00968928 3,78442142
2,9 3,78442142 0,005182648 0,001065025 0,00312384 3,78754526
3 3,78754526 0,001064446 2,72367E-05 0,00054584 3,7880911

c) RK de cuarto orden
METODO RK-4
t0: 0 tf: 3
y0: 1
h: 0,1 dy/dt=f(t,y)=y*(sen(t))^3

t y k1 k2 k3 k4 yn+1
0 1 0 1,24844E-05 1,24845E-05 9,95023E-05 1,00002491
0,1 1,00002491 9,95036E-05 0,000333746 0,000333786 2,6174E-07 1,00026404
0,2 1,00026404 0,000784345 0,001515322 0,001515875 3,91327E-06 1,00140582
0,3 1,00140582 0,002584471 0,004042633 0,004045572 2,39243E-05 1,00453662
0,4 1,00453662 0,005932189 0,008291073 0,008300779 9,18857E-05 1,01107125
0,5 1,01107125 0,011141541 0,014517563 0,014541668 0,000264677 1,0226587
0,6 1,0226587 0,018409896 0,022871245 0,022920688 0,000626695 1,04109544
0,7 1,04109544 0,027834824 0,033413393 0,033501733 0,001287544 1,06825421
0,8 1,06825421 0,03943474 0,046134279 0,046276322 0,002376088 1,10602622
0,9 1,10602622 0,053161174 0,060956006 0,061165762 0,0040308 1,15626547
1 1,15626547 0,068892983 0,077714106 0,07800197 0,006384096 1,22071701
1,1 1,22071701 0,08640747 0,096115906 0,096485049 0,009536245 1,30090794
1,2 1,30090794 0,105329238 0,115680209 0,116122521 0,013514403 1,39798279
1,3 1,39798279 0,125064964 0,135670745 0,136163343 0,018216522 1,5124744
1,4 1,5124744 0,144740963 0,15504504 0,155549063 0,023350038 1,64402094
1,5 1,64402094 0,163169698 0,172448662 0,172912309 0,028390797 1,79106801
1,6 1,79106801 0,178877787 0,186287628 0,186654647 0,032602101 1,95062875
1,7 1,95062875 0,19022576 0,194902527 0,195125311 0,035152909 2,11820114
1,8 2,11820114 0,195632153 0,196841356 0,196895061 0,0353418 2,28794227
1,9 2,28794227 0,193879416 0,191186637 0,191078703 0,032868148 2,45315531
2 2,45315531 0,184434826 0,17785022 0,177620181 0,028026301 2,60705564
2,1 2,60705564 0,167686744 0,157730696 0,157438902 0,021691852 2,74367527
2,2 2,74367527 0,144999661 0,132654204 0,132363442 0,015059214 2,85869096
 2,3 2,85869096 0,118541043 0,105088567 0,104846326 0,009236884 2,94996558
2,4 2,94996558 0,090912163 0,077703289 0,077531966 0,004902635 3,0176798
2,5 3,0176798 0,064685208 0,05290137 0,052799178 0,002182679 3,06405796
2,6 3,06405796 0,041974589 0,032446797 0,032396693 0,00077489 3,09279737
2,7 3,09279737 0,024143136 0,017261423 0,017242294 0,000200463 3,10835588
2,8 3,10835588 0,01168474 0,007398781 0,007393689 3,14734E-05 3,1152394
2,9 3,1152394 0,004266224 0,002152507 0,002151777 1,88389E-06 3,11738552
3 3,11738552 0,000876105 0,000238567 0,000238543 5,34603E-09 3,11769057

2.- Formule algoritmo, DF y codifique en MATLAB, una function para resolver cualquier EDO
usando el método de Runge-Kutta de 4to. Orden. Pruebe su programa (RK4.m) con la siguiente
ecuación diferencial:

Resuelva para el intervalo t = [0,10], utilizando un paso de h = 0.5. Realice una gráfica de la
solución.

PSEUDOCÓDIGO
PROGRAMA EN MATLAB

Respuesta es 82,0016
3.- Para simular una población se utiliza el modelo logístico:

donde p = población, kgm = tasa máxima de crecimiento en condiciones ilimitadas, y pmáx es la capacidad de carga. Simule
la
población mundial entre 1950 y 2000, utilizando el método numérico Heun de un solo paso programado en MATLAB.
Para la simulación, utilice las siguientes condiciones iniciales y valores de parámetros: p0 (en 1950) = 2555 millones de
personas, kgm = 0.026/año, y pmáx = 12 000 millones de personas. Haga que la función genere salidas que correspondan
a las fechas de los datos siguientes de población. Desarrolle una gráfica comparativa de la simulación frente a los
siguientes datos reales.

PSEUDOCÓDIGO DIAGRAMA DE FLUJO

PROGRAMA EN MATLAB
function [x,y] = Heun1(f,V,y0,n)
%MÉTODO DE Heun de un solo paso
para resolver dy/dx = f(x,y)
% entradas
% f: dy/dx
% V: vector de los valores inicial
y final de la variable
independiente
% p0: valor inicial de la variable
idepediente
% n: numero de subintervalos
% salidas
 % t: vector de la variable
independiente
% p: vector de
la variable
independiente
(vector
solución) x0 =
V(1);
xf = V(2);
h=(xf-x0)/n;
x=( x0:h:xf)';
nt=length(x);
y=y0*ones(nt,1);
for k = 2:nt
y(k) = y(k-1) + (h/2)*(f(x(k-
1),y(k-1))+f(x(k),y(k)));
end
fprintf('Año\n')
fprintf('%4.0f',x')
fprintf('\nPoblación\n')
fprintf('%4.1f',y')
fprintf('\n')
end
4.- Un tanque esférico tiene un orificio circular en el fondo a través del cual fluye líquido (véase la figura P25.24). La tasa
de flujo a través del agujero se calcula como:

donde Qsal = flujo de salida (m3/s), C = coeficiente obtenido en forma empírica, A = área del orificio (m2), g = constante
gravitacional (= 9.78 m/s2) y H = profundidad del líquido dentro del tanque. Emplee el programa RK4.m (Problema 2) a
fin de determinar cuánto tiempo tomaría que el agua fluyera por completo de un tanque de 3 m de diámetro con altura
inicial de 2.75 m, use un paso de 0.1. Observe que el orificio tiene un diámetro de 3 cm y C = 0.55. Realice una grafica de
Qsal vs. H.

PSEUDOCÓDIGO DIAGRAMA DE FLUJO

PROGRAMA EN MATLAB
5.- El siguiente es una ecuación diferencial de segundo orden con valor inicial:

Observe que w = 1. Descomponga la ecuación en dos ecuaciones diferenciales de primer orden. Después de la
descomposición, usando el metodo de Euler Modificado (Heun) resuelva el sistema de t = 0 a 15 con h = 0.5, y grafique
sus resultados.

PSEUDOCÓDIGO DIAGRAMA DE FLUJO

PROGRAMA EN MATLAB

6.- Resuelva el problema siguiente de valor inicial, en el intervalo de t = 2 a t = 3.

Utilice el método sin autoinicio de Heun con tamaño de paso de 0.5 y condiciones iniciales de y(1.5) = 5.800007 y y(2.0)
=
4.762673. Itere el corrector a s = 0.1%. Calcule los errores relativos porcentuales verdaderos t de sus resultados, con base
en la solución analítica.

PSEUDOCÓDIGO
PROGRAMA EN MATLAB
function [x,y] =
HeunSA(f,V,y0,y_1,h,tol) x0=V(1);
xf=V(2);
x=(x0:h:xf)';
nx = length(x);
y = y0*ones(nx,1);
B=[y_1;y];
for k = 1:nx-1
yp = B(k) + f(x(k),y(k))*2*h; error=1;
while error>tol
y(k+1)=y(k)+(h/2)*(f(x(k),y(k))+f(x(k+1),yp)); error=abs((y(k+1)-
yp)/y(k+1))*100;
yp=y(k+1);
end
end
end
7.- Resuelva el problema siguiente de valor inicial, de t = 4 a 5:

Use un tamaño de paso de 0.5 y valores iniciales de y(2.5) = 0.48, y(3) = 0.333333, y(3.5) = 0.244898, y y(4) = 0.1875.
Obtenga
sus soluciones con las técnicas siguientes: a) método de Heun sin autoinicio (s = 1%), y b) método de Adams de cuarto
orden
(s = 0.01%). [Nota: las respuestas correctas que se obtienen de forma analítica son y(4.5) = 0.148148 y y(5) = 0.12.]
Calcule los
errores relativos porcentuales verdaderos t de sus resultados.

PSEUDOCÓDIGO DIAGRAMA DE FLUJO

PROGRAMA EN MATLAB

8.- Solucione el problema siguiente de valor inicial, de t = 1.5 a t = 2.5

Use el método de Adams de cuarto orden. Emplee un tamaño de paso de 0.5 y el método de RK de cuarto orden para
pronosticar
los valores de inicio si y(0) = 2.

PSEUDOCÓDIGO DIAGRAMA DE FLUJO

PROGRAMA EN MATLAB

9.- Desarrolle un programa amigable para el usuario para el método de Heun sin autoinicio con modificador predictor.
Emplee
el método de RK de cuarto orden para calcular valores de inicio. Pruebe el programa con la repetición del problema 8.

10.- Un reactor de tanque agitado continuo (CSRT) es alimentado con un reactivo A de concentración CA [mol/m3] que
reacciona en forma elemental para dar uno B según la reacción: 𝐴 → 𝐵, con una expresión cinética de primer orden −𝑟𝐴 =
𝑘𝐶𝐴 . El volumen del reactor es V=4 m3, el flujo volumétrico de alimentación F=2.778*10 -4 m3/s, la concentración inicial
de A en el tiempo t =0 es CAo= 20,000 mol/m3 y la constante de la velocidad de reacción de primer orden es k=1.111*10 -4
/s. Programe un GUIDE para encontrar el comportamiento dinámico del reactor con respecto al reactivo A, hasta que la
concentración alcance un nivel estable, realice gráficos que muestren este comportamiento.

La ecuación diferencial a resolver es:

 𝑑 𝐶𝐴 𝐹
= (𝐶𝐴𝑜 − 𝐶𝐴 ) − 𝑘𝐶𝐴
𝑑𝑡 𝑉

En esta ecuación el volumen del líquido en el reactor es mantenido constante, en otras palabras el cambio de volumen como
resultado de la reacción química es despreciable. Sólo la concentración 𝐶𝐴 es función del tiempo, como lo es la temperatura.

Sustituyendo valores en esta ecuación diferencial se obtiene:

𝑑 𝐶𝐴
= 0.25 (20,000 − 𝐶𝐴 ) − 0.4 𝐶𝐴 ↪ 𝑃𝑎𝑟𝑎 𝑡 = 0, 𝐶𝐴 = 20,000
𝑑𝑡

Resuelva con el método de multipaso de Milne usando MATLAB (pag. 788 del libro base: Chapra). Realice una gráfica
de CA vs. T.

PSEUDOCÓDIGO DIAGRAMA DE FLUJO

PROGRAMA EN MATLAB