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

const int precision = 8;               
const int zmaxsize=10;
const int cmaxsize=0;

#include "genprec.h"
#include "cgenprec.h"
#include "apps.h"

int n,nmax;
int m,mmax;
int i,qqq;

realmp s[zmaxsizepad];
pademp d;                           //declare name of padeapproximant
realmp dnum[zmaxsizepad];     //declare name of numerator of pade
realmp dden[zmaxsizepad];     //declare name of denominator of pade

main()                    
  {n=1;
   m=1;
   cout<<'\n';
   cout<<"the example is (1,1) approximant to the series 1 + (1/2)x +(5/8)x^2"<<'\n';
   cout<<"for a continuation of this example see seriesx2.cpp"<<'\n';
   cout<<'\n';
   s[0]=1;            //1,1/2,-5/8
   s[1]=s[0]/2;
   s[2]=-(5*s[1]/4);    //the parentheses are necessary--?
   d<padeapprox(n,m,s);     //s is the series
   dnum<extractnum(n,d);
   dden<extractden(n,d);
   cout<<"the numerator of the Pade approximant is"<<'\n';
   printmp(dnum[0]);         
   printmp(dnum[1]);         
   cout<<'\n';
   cout<<"the denominator is"<<'\n';
   printmp(dden[0]);
   printmp(dden[1]);
   zmessage="  ";
   cerr<< zmessage << '\n';
   return 0;
  }