<<Graphics`Arrow`

Kolme liigi mudel (a=b=c=d=e=f=g=1; ga=fb)

a = 1;  
b = 1; 
c = 1; 
d = 1;
e = 1;
f = 1;
g = 1; 
sol = NDSolve[
{x'[t] == a x[t] - b x[t] y[t],
 y'[t] == - c y[t] + d x[t] y[t] - e y[t] z[t],
 z'[t] == - f z[t] + g y[t] z[t],
 x[0] == 0.5, y[0] == 1, z[0] == 2 }, {z, y, x}, {t, 0, 16} ]

 Plot[Evaluate[{x[t], y[t], z[t]} /. sol], {t, 0, 16}, PlotStyle -> {RGBColor[1,0,0], 
 RGBColor[0,1,0], RGBColor[0,0,1]},
   PlotRange->All, AxesStyle -> Arrowheads[{0.025}]]

gr1 = ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. sol], {t, 0, 16},PlotRange->All,AxesLabel->{"x","y","z"}]

TableForm[Table[{t,(x/.sol[[1]])[t],(y/.sol[[1]])[t], (z/.sol[[1]])[t]},{t,0,16
}],
TableHeadings->{None,{"t", "x[t]", "y[t]", "z[t]"}}]

Kolme liigi mudel(a=b=c=d=e=f=1, g=0.88; ga<fb)

a = 1;  
b = 1; 
c = 1; 
d = 1;
e = 1;
f = 1;
g = 0.88; 
sol = NDSolve[
{x'[t] == a x[t] - b x[t] y[t],
 y'[t] == - c y[t] + d x[t] y[t] - e y[t] z[t],
 z'[t] == - f z[t] + g y[t] z[t],
 x[0] == 0.5, y[0] == 1, z[0] == 2 }, {z, y, x}, {t, 0, 96} ]

 Plot[Evaluate[{x[t], y[t], z[t]} /. sol], {t, 0, 96}, PlotStyle -> {RGBColor[1,0,0], 
 RGBColor[0,1,0], RGBColor[0,0,1]},
   PlotRange->All, AxesStyle -> Arrowheads[{0.025}]]

gr1 = ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. sol], {t, 0,96},PlotRange->All,PlotStyle->{RGBColor[0,0,0.7]},AxesLabel->{"x","y","z"}]

punktid=Table[{(x/.sol[[1]])[t],(y/.sol[[1]])[t], (z/.sol[[1]])[t]},{t,0,96}]

TableForm[Table[{t,(x/.sol[[1]])[t],(y/.sol[[1]])[t], (z/.sol[[1]])[t]},{t,0,96
}],
TableHeadings->{None,{"t", "x[t]", "y[t]", "z[t]"}}]

Kolme liigi mudel(a=b=c=d=e=f=1, g=1.6; ga>fb)

a = 1;  
b = 1; 
c = 1; 
d = 1;
e = 1;
f = 1;
g = 1.6; 
sol = NDSolve[
{x'[t] == a x[t] - b x[t] y[t],
 y'[t] == - c y[t] + d x[t] y[t] - e y[t] z[t],
 z'[t] == - f z[t] + g y[t] z[t],
 x[0] == 0.5, y[0] == 1, z[0] == 2 }, {z, y, x}, {t, 0, 16} ]

 Plot[Evaluate[{x[t], y[t], z[t]} /. sol], {t, 0, 16}, PlotStyle -> {RGBColor[1,0,0], 
 RGBColor[0,1,0], RGBColor[0,0,1]},
   PlotRange->All, AxesStyle -> Arrowheads[{0.025}]]

gr1 = ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. sol], {t, 0,16
},PlotRange->All,AxesLabel->{"x","y","z"}]

TableForm[Table[{t,(x/.sol[[1]])[t],(y/.sol[[1]])[t], (z/.sol[[1]])[t]},{t,0,16}],
TableHeadings->{None,{"t", "x[t]", "y[t]", "z[t]"}}]