/*     ATK IV - Numerical Programming (2006)     */
/*      Excercise 4.3 -  Pade Approximant        */

#include<stdio.h>
#include<math.h>
#define NRANSI
#include<nr.h>
#define Pi 3.141592

int main(void){
  double coef[7]={0,1.000,0,0.333,0,0.133};
  double x=-Pi;
  double step = 0.1;
  int i, n=3;
  float resid;
  FILE *file;
  file=fopen("print.txt","w");

  pade(coef,n,&resid);
  
  for(i=0;i<=200;i++){
    fprintf(file,"%lf %lf %lf\n",x,ratval(x,coef,3,3),tan(x));
      x+=2.0*Pi/200.0;
    }
  fclose(file);
  return 0;
  
}