function dy=scm_fun(t,y)
dy = zeros(3,1);    % a column vector
% constant definition
tau_atm=2.82      % years
% tau_ocean

dy(1)= y(2) * y(3);
dy(2)= -y(1) * y(3);