Your first line should be @pyimport scipy.integrate as integrate. On Mon, Jun 20, 2016 at 5:32 PM, Henri Girard <henri.gir...@gmail.com> wrote:
> I am trying to convert from python to julia, but I don't know how to use > y1=integrate.odeint(f,t), apparently I should have add a derivativ ? > > --python---------------------- > from scipy.integrate import odeint > def f(y,t): > if t<25: > return [y[1],-8*y[0]+0.1*y[1]] > elif 25<t<45: > return [y[1],-8*y[0]] > else: > return [y[1],-8*y[0]-0.1*y[1]] > t=np.linspace(0,35,1000) > # start from y=0.01, y’=0 > y1=odeint(f,[0.01,0],t) > plt.plot(t,y1[:,0]) > plt.title("Evolution temporelle") > plt.xlabel("Temps t (s)") > plt.ylabel("Intensite i (A)") > plt.plot(0,0.01,"ro") > #plt.savefig("RLC-demarrage.eps") > plt.show() > ----------ijulia---------------- > @pyimport scipy.integrate as odeint > function f(y,t) > if t<25 > return [y[1],-8*y[0]+0.1*y[1]] > elseif 25<t<45 > return [y[1],-8*y[0]] > else > return [y[1],-8*y[0]-0.1*y[1]] > end > end; > t=linspace(0,35,1000) > y=linspace(0.01,0) > y1=integrate.odeint(f,y,t) >