%
%
function [f1,f2]=bobs_dop0_c(u,s,dt,t,flg)
%
%  s = [r q]'
%
r = s(1); q = s(2); 
%
if flg == 1
   f1 = [r-dt*tan(u); q+2*dt*sin(u)^2*tan(u)*(2*r-dt*tan(u))];  
elseif flg == 2
   f1 = 2*r^2+q; f2 = [4*r 1];
elseif flg == 3
    f1 = [1 0; 4*dt*sin(u)^2*tan(u) 1]; 
    f2 = [-dt*sec(u)^2; dt*4*r*(3*sin(u)^2+sin(u)^2*tan(u)^2)-4*dt*dt*(2*sin(u)^2*tan(u)+tan(u)^3*sin(u)^2)];
end;
      