function Xpr = Sprott_G(t,X)

%%% Created by ODESOLVE
%
%Input:
%   Independent variable: t
%   Dependent variable vector: X = [x, y, z]
%
%Output:
%   Vector of derivatives: Xpr
%
%To solve the system, execute
%      [t,u] = ode45(@Sprott_G, tspan, u0)
%   tspan is the solution interval
%   u0 is the initial vector


Xpr = [2*X(1)/5+X(3);...
       X(1)*X(3)-X(2);...
       -X(1)+X(2);];