#ifndef myfunc_h
#define myfunc_h

/* mycos1 and mysin1: ordinary taylor expansions of sine and cosine.
   End the expansion when a term is smaller than 10e-18 of the accumulated series
  Note: pi is given in 21 decimals here
*/

#define PI 3.14159265358979323846 

double mycos1(double x);
double mysin1(double x);
float mycos2(float x);
float mysin2(float x);
double fac(int i);
double mycos3(double x);
double mysin3(double x);
void powercount();

#endif

