func1 = '10*(-x+y)';
func2 = '24*x-y-x*z';
func3 = '-2.5*z+x*y';
[T,Y1,Y2,Y3]=ode45_3(func1, func2, func3, [0, 20], [5; 7; 9]);
