void cmult_(x,y,z)
float *x, *y, *z;
{
  *z = *x * *y + 2.0;
  return;
}