function robeq=robeq(x);
   global ak;
   robeq=(ak-1)*(phi(x)+x*cumnord(x))-ak*x
end;