function ohat = freq(a,algor)
%
%  This function computes the frequency of an
%  ODE solver
%
f = ff2(a,1,algor);
ohat = atan2(imag(f),real(f));
%
return
