function dy_dt = D1( t, y, par ) 
   dy_dt = 2 * t; 
end