Thanks much better : The worse I have done day21of julia and it worked but
converted to mine it didn"t :
Has ode23 any importance ?
I have looked two days on most ODE help (there is more than the manual with
example) but I really got grasy.
Son now I have to finish phase portrait, done
Solved thanks to your help :)
using ODE,Plots
pyplot(size=(300,200),leg=false,
guidefont=font(7), titlefont=font(7));
function oscillator(t, y)
y[2] = - 3* + y[1] - y[2] / 10
end
oscillator(t, y) = [y[2], - 3* + y[1] - y[2] / 10]
initial = [1.0,0.0]; # Initial conditions
t = float([0:0.01:50]);
t,y = ode23(oscillator, initial, t)
y=hcat(y...).';
plot(t,y,title="Mouvement",
bg=RGB(.2,.2,.2),
xlabel ="Temps",ylabel = "Vitesse")
