RESOLVER LAS SIGUIENTES EDOS EN MAT LAB
EJERCICIO 6.1
%% El comando dsolve permite hallar la solución de muchas edo’s
eq_dif:=(1+exp(x))*diff(y(x),x)=exp(x)/y(x);
dsolve({eq_dif,y(0)=1});
%% Además Maple nos permite saber el tipo de ecuación diferencial
%% con la que estamos tratando, cosa que es muy útil a la hora
%% de resolverlas a mano
with(DEtools);
odeadvisor(eq_dif);
eq_dif1:=(1+exp(x))*(diff(y(x),x)) - exp(x)/y(x)=0;
%% Además podemos forzar a que nos de mas informacion de como se resuelve
%% la ecuación diferencial antes de imponer la condicion inicial
dsolve(eq_dif, [separable], useInt);
value(%);
%% Finalmente tenemos que imponer la condición inicial
%% para determinar la constante _C1
subs(x=0,%);
subs(y(0)=1,%);
solve(%,_C1);
%% El comando DEplot permite dibujar el campo de velocidades y la
solución
DEplot (eq_dif, y(x), x=o..100,[{ y(o) = 1}], y= -1..15 linecolour = red,
color = blue, sepsize =.1,arrows = MEDIUN );
eq_dif, y(x), x=0..100, [[ y(0) = 1 = blue, stepsize=.1,arrows=MEDIUM );
EJERCICIO 6.2
eq_dif:=diff(y(x),x)-y(x)*cot(x)=2*x*sin(x);
dsolve(eq_dif);
with(DEtools):
odeadvisor(eq_dif);
%% Dibujamos la solución para y(1)=1
DEplot( eq_dif, y(x), x=0..2*Pi, [[ y(1) = 1 ]],y=-10..5,
linecolour=red,
color = blue, stepsize=.1,arrows=MEDIUM );
EJERCICIO 6.3
eq_dif:=diff(y(x),x)+2*x*y(x)=x*y(x)^4;
dsolve(eq_dif);
sol:=rhs(%[1]);
simplify(sol^3);
simplify(1/sol^3);
%% Vamos a ver que tipo de cuacion es
with(DEtools):
odeadvisor(eq_dif);
dsolve(eq_dif, [separable] , useInt);
value(%);
%% Si queremos resolverla como una Bernouilli:
sol2:=rhs(%[1]);
value(sol2);
simplify(expand(1/%^3));
%% Dibujamos la solución para y(0)=-1, y(0)=1, y(0)=2^(1/3) y para
y(0)=3/2
DEplot( eq_dif, y(x), x=0..8, [[y(0)=-
1],[y(0)=1],[y(0)=2^(1/3)],[y(0)=3/2]],
y=-1..4,linecolour=[red,green,yellow,pink],color = blue,
stepsze=.01,arrows=MEDIUM )
%% Notar que para y(0)=2^(1/3), la % El dibujo que hemos obtenido (curva
%% bfunción deja de ser constante. % Esto es debido a que las curvas
solucón se calculan
%% y no utilizando la solución analítica # Representemos las soluciones
analíticas
%%y(0)=-1 solución seria _C1=0 -> y(x) = 2^(1/3) amarilla) vemos que
antes de x=4 la numéricamente
-1.500000000
sol1:=subs(_C1=-1.5,sol);
%% y(0)=1
fsolve(subs(x=0,sol)=1,_C1);
sol2:=subs(_C1=0.5,sol);
%% y(0)=2^(1/3)
fsolve(subs(x=0,sol)=2^(1/3),_C1)
sol2:=subs(_C1=0.5,sol);
%% y(0)=2^(1/3)
fsolve(subs(x=0,sol)=2^(1/3),_C1);
sol3:=subs(_C1=-0,sol);
%% y(0)=3/2
fsolve(subs(x=0,sol)=3/2,_C1);
sol4:=subs(_C1=%,sol);
plot([sol1,sol2,sol3,sol4],x=0..8,y=-1..4,colour=[red,green,yellow,pink],
thickness=2,discont=true);
EJERCICIO 6.4
eq_dif:=diff(y(x),x,x)+4*diff(y(x),x)+4*y(x)=exp(2*x)*(16*x^2+16*x-14);
dsolve(eq_dif);
with(DEtools):
odeadvisor(eq_dif);
%% Como tenemos una ecuación de segundo orden nesecitamos 2 condiciones
%% adicionales para determinar la solución
%% Vamos a dibujar las soluciones asociadas a y(0)=-1,y’(0)=-2 y
%%y(0)=0, y’(0)=0
DEplot(eq_dif,y(x),x=0..1.5,y=-5..10,[[y(0)=-1,D(y)(0)=-
2],[y(0)=0,D(y)(0)=0]],
stepsize=0.01,arrows=medium,color = blue,linecolor=[red,blue]);
EJERCICIO 6.5
eq_dif:=diff(y(x),x,x)-3*diff(y(x),x)+2*y(x)=x^2*exp(x)+5*cos(x);
dsolve(eq_dif);
%% vamos a ver que tipo de ecuación es
[[ _ 2nd _ order , _ linear , _ nonhomogeneous ]]
with(DEtools):
odeadvisor(eq_dif);
dsolve(eq_dif, [linear], useInt);
EJERCICIO 6.6
eq_dif:=diff(y(x),x,x)+y(x)=1/cos(x);
dsolve(eq_dif);
%% Vamos a ver que tipo de ecuación es
with(DEtools):
odeadvisor(eq_dif);
dsolve(eq_dif, [linear], useInt);
[[ _ 2nd _ order , _ linear , _ nonhomogeneous ]]
EJERCICIO 6.7
sys_ode := diff(x(t),t)-3*x(t)+2*y(t)=0,diff(y(t),t)-2*x(t)+y(t)=0;
dsolve([sys_ode]);
%% Dibujamos las soluciones asociadas a [x(0)=1,y(0)=0] y
[x(0)=1.5,y(0)=2]
with(DEtools):
DEplot( [sys_ode], [x(t),y(t)], t=0..10, x=-5..5, y=-5..5,
[[x(0)=1.5,y(0)=2],[x(0)=1,y(0)=0]],stepsize=0.01,arrows=medium,
color = blue,linecolor=[red,blue]);
Fabian Loor.
