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'.
José Jeremías Caballero
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
