In[]:=
MakePlot[{α_,β_,γ_},{a_,b_,c_},h_,max_Integer]:=Module[{u,if},f[{x_,y_,z_}]:={z+(y-α)*x,1-β*y-x^2,-x-γ*z};​​u[0]={a,b,c};​​Do[u[n+1]=u[n]+h*f[u[n]+h/2*f[u[n]]],{n,0,max}];​​if=Interpolation[Table[{n,u[n]},{n,0,max/h}]];​​Plot[if[t],{t,0,200},PlotLegends->Automatic]]
In[]:=
plot1=MakePlot[{0.9,0.2,1.2},{1,3,2},0.1,200]​​plot2=MakePlot[{0.9,0.2,0.4},{1,3,2},0.1,200]​​
Out[]=
50
100
150
200
-2
-1
1
2
3
Out[]=
50
100
150
200
-1
1
2
3
In[]:=
Show[{plot1,plot2}]
Out[]=
50
100
150
200
-2
-1
1
2
3