#include <stdio.h>
void cblas_drotmg (double *d1, double *d2, double *x1, const double y1, double *param);

void drotmg_ (double *d1, double *d2, double *x1, const double* y1, double *param);

int main()
{
  int i;
  double d1,d2,x1,y1;
  double param[5];

  /*
  d1=2;
  d2=3;
  x1=5;
  y1=8;
  */

  d1=2;
  d2=2;
  x1=8;
  y1=8;

  //  cblas_drotmg(&d1, &d2, &x1, y1, param);
  drotmg_(&d1, &d2, &x1, &y1, param);
  
  printf("%lf\n", d1);
  printf("%lf\n", d2);
  printf("%lf\n", x1);

  for(i=0; i<5; i++) {
    printf("param[%d]=%lf\n", i, param[i]);
  }
  return 0;
}