ndist=NormalDistribution[0,1];
phi[x_] := PDF[ndist,x];
robeq[x_] :=(ak-1)*(phi[x]+x*CDF[ndist,x])-ak*x;