%% Plot trajectory of Lorenz system
 clear; clf;
 global a b c;
 a=10; b=28; c=8/3;
 u0 = zeros(3,1)+0.5; param=0; t=(0:0.02:100)';
 u=lsode('lorenz_odeOctave',u0,t);
 plot3(u(:,1),u(:,2),u(:,3))
