function gamma_dot = fcn(p,phi_phi_t)
%#codegen
R4 = 5000;
beta = 0.1;
if abs(p) > R4
    gamma_d = 0;
else
    gamma_d = beta*p-(p*phi_phi_t*p)/(1+phi_phi_t);
end

gamma_dot = gamma_d;