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

int precision = 8;               
const int zmaxsize=40;
const int cmaxsize=40;

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

int n;
int nmax;

realmp p;
realmp a[zmaxsizepad];
realmp b[zmaxsizepad];
realmp c[zmaxsizepad];
complexmp cc[zmaxsizepad];

int main()
  {int i;
   zsupe();                 //run this if e is needed
   zsuplogsp();             //run this if log base e of 10 is needed
   n=38;
                           //compute series for sqrt(1+x)
   b[0]=1;
   b[1]=1;
   for (i=2;i<n+1;i++)     //remaining terms in a,b are zero
     {b[i]=0;}             //end set up
   p=1;                  
   p=p/2;
   c<stonth(n,b,p);           //b^(1/2)
   c<slogder(n,c);            //take log derivative
   for (i=0;i<n;i++)           //not n+1 because derivative
     {cout<<i<<'\n';           //has decreased order by 1
      printmp(c[i]);
      zmessage="should be terms in expansion of (1/2)/(1+x)";
      zbreak();
     }

//same thing another way:

   p=-1;
   p=p/2;
   c<stonth(n,b,p);   //1/the fcn; i.e. (1+x)^(-1/2)
   p=1;
   p=p/2;              //the fcn: i.e. (1+x)^(1/2)
   a<stonth(n,b,p);
   a<sder(n,a);       //derivative of the fcn
   c<smult(n-1,a,c);      // df/dx/f
   for (i=0;i<n;i++)           //not n+1 because derivative
     {cout<<i<<'\n';           //has decreased order by 1
      printmp(c[i]);
      zmessage="same thing another way";
      zbreak();
     }

   realmp fac;
   fac=2;
   c<smultfac(n-1,c,fac);
   for (i=0;i<n;i++)
     {cout<<i<<'\n'; 
      printmp(c[i]);
      zmessage="2 x last result";
      zbreak();
     }

   cc<scontoc(n-1,c);     //convert that to complex a complex series

   for (i=0;i<n;i++)
     {cout<<i<<'\n'; 
      printmp(cc[i]);
      zmessage="that series converted to complex";
      zbreak();
     }

   fac=1;
   fac=fac/10;
   fac=seval(n-1,c,fac);    //1/(1-x) for x=0.1
   printmp(fac);
   zmessage="series for 1/(1-x) evaluated at 0.1";
   cerr << zmessage << '\n';

   return 0;
  }