clear; a = 0.4; b = 0.7; c = 0.0; d = 1.0 function dx = silnikov(t,x) dx(1) = x(2); dx(2) = x(3); dx(3) = - a * x(3) - x(2) + b * x(1) * (1 - c * x(1) - d * x(1) ^ 2); endfunction T = linspace(0,200,20000); x0 = 0.1; y0 = 0.1; z0 = 0.2; X0 = [x0; y0; z0]; X = ode(X0,0,T,silnikov); plot3d3(X(1,:),X(2,:),X(3,:));