jueves, 5 de abril de 2012

Ecuaciones Diferenciales ode45

Estoy trabajando con el comando ode45, el cual requiere la asignación de una función en un archivo .m. Me encuentro definiendo dicha función f.m donde aparece la ecuación diferencial a resolver por ode45, pero ahora necesito que en esta ecuación aparezca el valor de una variable de entrada definida en el Command Window principal, como por ejemplo 'a'.
function ode45caballero5
global a
a=input('ingrese el valor de a:');
% Ejemplo
% y'=a+2t+y 1<=t<=5
% y(1)=0.5
[T,Y]=ode45(@f,[1,5],0.5);
plot(T,Y)
function dydt=f(t,y)
global a
dydt=a+2*t+y;
EJECUCION EN MATLAB
>> ode45caballero5
ingrese el valor de a:10



function ode45caballero6
global a
a=input('ingrese el valor de a:');
[t,x]=ode45(@fint2,[0,5],[-4,0]);
figure(gcf)
plot(t,x(:,1));
end

function dydt=fint2(t,x)
global a
dydt=[x(2);a-x(2)];
end
EJECUCION EN MATLAB
>> ode45caballero6

ingrese el valor de a:6

José Jeremías Caballero
Servicios de Programación en MatLab
Lima - Perú
Para mayor informacion, contactarse a jjcc94@hotmail.com

No hay comentarios.:

Publicar un comentario