RE: ode45 for coupled system (2nd order ODE)

[Joza]

Here is my code anyway:

function du = f1(t, y) 
        dy = zeros(2,1); 
        dy(1) = y(2);
        dy(2) = -sin(y(1));
end

...

[t, y1] = ode45(@f1, [0, 100], [1, 0]);

This solves the system for y1. But I would like to also have dy1/dt .... how
can I obtain that with ode45?



--
View this message in context: 
http://octave.1599824.n4.nabble.com/ode45-for-coupled-system-2nd-order-ODE-tp4646324p4646344.html
Sent from the Octave - General mailing list archive at Nabble.com.