%The phase error for Runge_Kutta 1 is given by 
%
%erp=1-arctan(wh)/wh