#include <imsl.h>
extern double ak;
extern double phi(double f);
double robeq(double f)
{
  double r=(ak-1)*(phi(f)+f*imsl_d_normal_cdf(f))-ak*f;
  return r;
}