EJERCICIOS EDO
ESTUDIANTE: NUÑEZ HUAMANI ANTHONY
CURSO: MÉTODOS NUMÉRICOS
G.H.: 01L
25.1 RESUELVA EN FORMA ANALÍTICA EL PROBLEMA DE VALORES
INICIALES SIGUIENTE, EN EL INTERVALO DE x=0 A 2:
dy
yx 2 1.1y
dx
DONDE y(0)=1. GRAFIQUE LA SOLUCIÓN.
𝑑𝑦
∫ = ∫ 𝑥 2 𝑑𝑥 − 1.1 ∫ 𝑑𝑥
𝑦
𝑥3
ln 𝑦 = − 1.1𝑥 + 𝐶
3
CONSIDERAMOS EL VALOR DE C=0 INICIALMENTE
𝑥3
𝑦 = 𝑒 ( 3 −1.1𝑥)
EN MATLAB:
SCRIPT
%%PROB_25_1
clc
clear all
x=0:0.25:2;
ya=exp((x.^3/3)-1.1*x);
plot(x,ya);
grid on
GRÁFICA:
25.2 UTILICE EL MÉTODO DE EULER CON h=0.5 Y 0.25, PARA
RESOLVER EL PROBLEMA 25.1. GRAFIQUE LOS RESULTADOS EN LA
MISMA GRÁFICA PARA COMPARAR EN FORMA VISUAL LA
EXACTITUD DE LOS DOS TAMAÑOS DE PESO.
SCRIPT
CON h=0.5
%PROB_25_2_A
a= input('ingresar el valor de t1: ');
b= input('ingresar el valor de tN: ');
h=0.5;
N=[(b-a)/h];
y(1)=1;
x=a:h:b;
for i=1:N
dy=y(i)*(x(i).^2)-1.1*y(i);
y(i+1)=y(i)+h*dy;
end
y(i+1)=y(i)+h*dy
plot(x,y,'r*');
grid on
ya=exp((x.^3/3)-1.1*x)
plot(x,ya,'k*')
grid on
legend('Euler','Analítico')
GRÁFICA
COMMAND WINDOW
ingresar el valor de t1: 0
ingresar el valor de tN: 2
y=
1.0000 0.4500 0.2587 0.2458 0.3872
ya =
1.0000 0.6015 0.4646 0.5916 1.5947
SCRIPT
CON h=0.25
%PROB_25_2_B
a= input('ingresar el valor de t1: ');
b= input('ingresar el valor de tN: ');
h=0.25;
N=[(b-a)/h];
y(1)=1;
x=a:h:b;
for i=1:N
dy=y(i)*(x(i).^2)-1.1*y(i);
y(i+1)=y(i)+h*dy;
end
y(i+1)=y(i)+h*dy
plot(x,y,'r*');
grid on
ya=exp((x.^3/3)-1.1*x)
plot(x,ya,'k*')
grid on
legend('Euler','Analítico')
GRÁFICA
COMMAND WINDOW
ingresar el valor de t1: 0
ingresar el valor de tN: 2
y=
Columns 1 through 5
1.0000 0.7250 0.5370 0.4229 0.3660
Columns 6 through 9
0.3569 0.3981 0.5126 0.7641
ya =
Columns 1 through 5
1.0000 0.7635 0.6015 0.5044 0.4646
Columns 6 through 9
0.4848 0.5916 0.8706 1.5947
28.7 LAS ECUACIONES SIGUIENTES DEFINEN LA CONCENTRACIÓN DE
3 REACTIVOS:
dca
10ca cc cb
dt
dcb
10ca cc cb
dt
dcc
10ca cc cb 2cc
dt
SI LAS CONDICIONES INICIALES SON DE ca=50, cb=0 y cc=40,
ENCUEENTRE LAS CONCENTRACIONES PARA LOS TIEMPOS DE 0 a 3s.
%PROB_28_7
clc
clear all
a= input('ingresar el valor de t1: ');
b= input('ingresar el valor de tN: ');
h=3;
N=[(b-a)/h];
ca(1)=50;
cb(1)=0;
cc(1)=40;
t=a:h:b;
for i=1:N
dca=-10*ca(i)*cc(i)+cb(i);
dcb=10*ca(i)*cc(i)-cb(i);
dcc=-10*ca(i)*cc(i)-2*cc(i);
ca(i+1)=ca(i)+h*dca;
cb(i+1)=cb(i)+h*dcb;
cc(i+1)=cc(i)+h*dcc;
end
ca
cb
cc
COMMAND WINDOW
ingresar el valor de t1: 0
ingresar el valor de tN: 3
ca =
50 -59950
cb =
0 60000
cc =
40 -60200
>>
28.28 PARA UN CIRCUITO SENCILLO RL, LA LEY DE KIRCHOFF DEL
VOLTAJE REQUIERE QUE (SI SE CUMPLE LA LEY DE OHM)
dI
L RI 0
dt
DONDE I= CORRIENTE, R=1.5 E I(0)=0.5. RESUELVA ESTE PROBLEMA
EN FORMA ANALÍTICA Y CON ALGÚN MÉTODO NUMÉRICO.
PRESENTE SUS RESULTADOS EN FROMA GRÁFICA.
SCRIPT
%PROB_28_28
clc
clear all
a= input('ingresar el valor de t1: ');
b= input('ingresar el valor de tN: ');
h=0.1;
N=[(b-a)/h];
I(1)=0.5;
t=a:h:b;
for i=1:N
dI=-1.5*I(i);
I(i+1)=I(i)+h*dI;
end
I(i+1)=I(i)+h*dI
plot(t,I,'r*')
grid on
Ia= exp(-1.5*t)
plot(t,Ia,'k*')
grid on
legend('Euler','Analítico')
COMMAND WINDOW
ingresar el valor de t1: 0
ingresar el valor de tN: 3
I=
Columns 1 through 5
0.5000 0.4250 0.3612 0.3071 0.2610
Columns 6 through 10
0.2219 0.1886 0.1603 0.1362 0.1158
Columns 11 through 15
0.0984 0.0837 0.0711 0.0605 0.0514
Columns 16 through 20
0.0437 0.0371 0.0316 0.0268 0.0228
Columns 21 through 25
0.0194 0.0165 0.0140 0.0119 0.0101
Columns 26 through 30
0.0086 0.0073 0.0062 0.0053 0.0045
Column 31
0.0038
Ia =
Columns 1 through 5
1.0000 0.8607 0.7408 0.6376 0.5488
Columns 6 through 10
0.4724 0.4066 0.3499 0.3012 0.2592
Columns 11 through 15
0.2231 0.1920 0.1653 0.1423 0.1225
Columns 16 through 20
0.1054 0.0907 0.0781 0.0672 0.0578
Columns 21 through 25
0.0498 0.0429 0.0369 0.0317 0.0273
Columns 26 through 30
0.0235 0.0202 0.0174 0.0150 0.0129
Column 31
0.0111
GRÁFICA
%% EJERCICIOS HECHOS CON EL OBJETIVO DE COMPLETAR DE LO HECHO
EN CLASE
