function  dy_dt = DglFun( t, y, k )
   dy_dt = - k * y; 
