#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;

complexmp p; 
complexmp a[cmaxsizepad];
complexmp b[cmaxsizepad];
complexmp c[cmaxsizepad];     

int main()
  {int i;
   zsuppi();
   zsupe();                 //run this if e is needed
   zsuplogsp();             //run this if log base e of 10 is needed
   n=4;       
   a[0]=cinttomp(1,0);       //compute series 
   a[1]=cinttomp(2,1);
   for (i=2;i<n+1;i++)       //remaining terms in a are zero
     {a[i]=cinttomp(0,0);}            
   p=cinttomp(1,0);
   p=p/2;
   a<cstonth(n,a,p);            //sqrt(1+(2+i)x)

   for (i=0;i<n+1;i++) printmp(a[i]);
   zmessage="coeffs in sqrt(1+(2+i)x)";
   zbreak();

   b[0]=cinttomp(1,0);         //b=1+(1+i)x
   b[1]=cinttomp(1,1);
   for (i=2;i<n+1;i++)
     {b[i]=cinttomp(0,0);}
   p=-p;
   b<cstonth(n,b,p);            //1/sqrt(b)

   for (i=0;i<n+1;i++) printmp(b[i]);
   zmessage="coeffs in 1/sqrt(1+(1+i)x)";
   zbreak();
 
//now add and subtract these series, multiply one by a factor,
//and take clogder of one

   c<csadd(n,a,b);
   for (i=0;i<n+1;i++) printmp(c[i]);
   zmessage="coeffs in sum of these series";
   zbreak();

   c<cssub(n,a,b);
   for (i=0;i<n+1;i++) printmp(c[i]);
   zmessage="coeffs in difference of these series";
   zbreak();


   complexmp fac;
   fac<cinttomp(2,0);
   c<csmultfac(n,c,fac);
   for (i=0;i<n+1;i++) printmp(c[i]);
   zmessage="coeffs in that multiplied by 2,0";
   zbreak();

   c<cslogder(n,a);
   for (i=0;i<n;i++) printmp(c[i]);   //not n+1 because of der
   zmessage="coeff in log der of sqrt(1+(2+i)x)";
   zbreak();

   a[0]=cinttomp(1,0);       //this will be a 1st degree poly 
   a[1]=cinttomp(2,1);

   b[0]=cinttomp(1,0);         //b=1+(1+i)x also a 1st degree poly
   b[1]=cinttomp(1,1);

   c<cmultpolys(1,a,1,b);      //product of two 1st degree polys

   for (i=0;i<3;i++) printmp(c[i]);
   zmessage="product of polys 1+(2+i)x and 1+(1+i)x";
   cerr << zmessage <<'\n';
   return 0;
  }