- —————————————————————————————– 99/30/08 9/30/08) >> y_int=0.1; >> [t,y]=euler(‘test_example’,[0.1,9.0],y_int,200); >> y_init=[.9,1,1.1]; >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,1000); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,2500); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,3000); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,500); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,800); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,700); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,750); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,799); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,800); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,850); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,830); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,810); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,805); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,1500); >> plot(t,y) >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,1500); >> [t,y]=euler_system(‘lorenz_system’,[0.0,20.0],y_init,1000); >>———————————————————————————— we created m files g.m function xdot = g(t,x) xdot [...]
