Hola: Tengo que resolver el siguiente PVI en matlab mediate elcomando ode45 para resolver el problema en el intervalo [0,1] .
y’’(t) + y(t) – sin(t) = 0, y(0) = 0, y’(0)=0.
El problema es que el comando tiene la sintaxis:
[T, Y]=ode45(f, tspan,y_0).
se me ocurrió hacer:
syms y t a
Dy=a;
y=int(a);
f=inline(vectorize(sin(t)-a))
[T,A]=ode45(f,[0,1],0);
pero ahora tengo que deshacer el cambio de variable, y tengo
y'= a;...... y(0)=0, pero ahora a es la matriz A..... entonces como puedo hacer¿? porque si pongo:
for i=1:41
[T,Y]=ode45(A(i),[0,1],0);
end
(41 es el tamaño de A en filas)
me da error...............
¿Alguien podría echarme una mano con esto?
Muchas gracias.
Iria.
