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

#define N 64
#define ZERO 0.0
#define ONE 1.0

extern double dlamch_();

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

  static int i,j,info,n=N,ldu=N;
  static double u[N*N],work[N*N],b[N],c[N],work2[8*N];

  for(i=0;i<N;i++){
    b[i]=sqrt(i+1);
  }

  for(i=0;i<N-1;i++){
    c[i]=sqrt(i+1);
  }
  c[N-1]=ZERO;

  for(i=0;i<N;i++){
    for(j=0;j<N;j++){
      u[i*N+j]=ZERO;
    }
  }

  for(i=0;i<N;i++){
    u[i*N+i]=ONE;
  }

  doqds2_("L", &n, b, c, u, &ldu, work, work2, &info);

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

  for(i=0;i<N;i++){
    printf("%30.20g\n",b[i]*b[i]);
  }

  for(i=0;i<N;i++){
    printf("%30.20g\n",u[N*i]);
  }

  return 0;
}