for i=1:3 for j=1:3 for k=1:3 [T,Y]=ode45('g',[0,2], [0.01*(i-2),0.01*(j-2),0.01*(k-2)]); plot3(Y(:,1),Y(:,2),Y(:,3)); hold on end end end