function xprime = f(t,x)
global a
% 
% right hand of an equation
%
xprime = -2*t.*x + sin(t); 