function [ dy ] = diffsyst( t,x )
dy=[    -10*x(1)+28*x(2);       ...
        10*x(1)-x(2)-x(1)*x(3); ...
        -8/3*x(3)+x(1)*x(2)     ...
   ];
end

