function X = BSCgeneric(t,x)

% The first example of specific equations undergoing the catastrophe 
% (Blue Sky Catastrophe) was given by N. Gavrilov and A. Shilnikov 
% [Gavrilov and Shilnikov, 2000; L. Shilnikov et al, 2001]
%
%   dx/dt = x(2+mu-beta((x^2)+(y^2)))+(z^2)+(y^2)+2*y
%   dy/dt = -(z^3)-(1+y)((z^2)+(y^2)+2y)-4x+ mu y
%   dz/dt = (1+y)(z^2)+(x^2)-eta.


global mu beta eta 

X=[ x(1)*(2+mu-beta*((x(1)^2)+(x(2)^2)))+(x(3)^2)+(x(2)^2)+2*x(2); 
-(x(3)^3)-(1+x(2))*((x(3)^2)+(x(2)^2)+2*x(2))-4*x(1)+mu*x(2); 
(1+x(2))*(x(3)^2)+(x(1)^2)-eta ];



% "Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009
% Copyright(c).