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

#define N 70000

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

  static int i,info,n=N;
  static double b[N],c[N],b2[N],c2[N],work[4*N];

  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<N;i++){
    b2[i]=b[i];
  }

  for(i=0;i<N;i++){
    c2[i]=c[i];
  }

  dlasq1_(&n, b, c, work, &info);
        
  printf("%30.20g\n",b[N-1]);

  ddqds_(&n, b2, c2, work, &info);
        
  printf("%30.20g\n",b2[N-1]);

  return 1;
}