% predprey_rh.m - feedback function for predator prey system
% RMM, 13 Sep 06
function rh=predprey_rh(t, x, par)
global predprey_K predprey_xd predprey_ud;

% Compute the feedback control law
rh = -predprey_K * ([x(1); x(2)] - predprey_xd) + predprey_ud;
if (rh < 0) rh = 0; end;
if (rh > 4*par) rh = 4*par; end;
