#include<stdio.h>
#include<stdlib.h>

#define N 600

int main(int argc,char *argv[]){

  static int i,info,n=N,n2=2*N;
  static double b[N],c[N];

  static char RANGE='A',JOBZ='V';
  static int m,nzc=2*N,il,iu,tryrac=1;
  static int isuppz[2*2*N],ldz=2*N,iwork[10*2*N],liwork=10*2*N,lwork=18*2*N;
  static double abstol,d[2*N],e[2*N],w[2*N],vl,vu,*z,work[18*2*N];

  z=malloc(sizeof(double)*(2*N)*(2*N));

  if (z==NULL) return 0;

  srand(1);

  for(i=0;i<N;i++){
    b[i]=rand()/(double)RAND_MAX;
    if (rand() % 2==0){
      b[i]=-b[i];
    }
  }
  
  for(i=0;i<N-1;i++){
    c[i]=rand()/(double)RAND_MAX;
    if (rand() % 2==0){
      c[i]=-c[i];
    }
  }
  c[N-1]=0;
      
  for(i=0;i<2*N;i++){
    d[i]=0;
  }

  for(i=0;i<N;i++){
    e[2*i]=b[i];
    e[2*i+1]=c[i];
  }

  dstemr_(&JOBZ,&RANGE,&n2,d,e,&vl,&vu,&il,&iu,&m,w,
	  z,&ldz,&nzc,isuppz,&tryrac,work,&lwork,iwork,&liwork,&info);

  printf("INFO=%d\n",info);

  return 0;
}