Mathc complexes/a144

Bibliothèque

Installer ce fichier dans votre répertoire de travail.

	wzst.h

/* ------------------------------------ */
/*  Save as :   wzst.h                  */
/* ------------------------------------ */
/* ------------------------------------ */
nb_Z pow_Z(
int     n,
nb_Z    A
)
{
nb_Z  A_n = A; 	
int         i;

  if(!n) A_n = i_Z(1,0);
     
  else
     
     for(i=n-1; i; i--)

             A_n = mul_Z(A_n,A);
      
  return(A_n);
}
/* ------------------------------------ */
/* ------------------------------------ */
nb_Z sin_Z(
nb_Z A
)
{	
nb_Z  A_n;
nb_Z   sA;
nb_Z sinA = i_Z(0,0);

int n = 12; 
int i;
int j;  
	
  for(i=1.,j=0.; i <= n; i+=2.,j++)
   {
     A_n  = pow_Z(i, A);         
     sA   = smul_Z(pow((-1),j)/factorial(i), A_n);
     sinA = add_Z(sinA,sA);       
   }
 
  return (sinA);	
}
/* ------------------------------------ */
nb_Z cos_Z(
nb_Z A
)
{	
nb_Z   A_n;
nb_Z   sA ;
nb_Z cosA = i_Z(0,0);


int n = 12; 
int i; 
int j; 
	
  for(i=0.,j=0.; i <= n; i+=2.,j++)
   {
     A_n  = pow_Z(i, A);         
     sA   = smul_Z(pow((-1),j)/factorial(i), A_n);
     cosA = add_Z(cosA,sA);  
   }
   
  return (cosA);	
}
/* ------------------------------------ */
/* ------------------------------------ */
nb_Z sinh_Z (
nb_Z A
)
{	
nb_Z  A_n;
nb_Z   sA;
nb_Z sinhA = i_Z(0,0);

int n = 12; 
int i; 
	
  for(i = 1.; i <= n; i+=2.)
   {
       A_n = pow_Z(i, A);         
        sA = smul_Z(1./factorial(i), A_n);
     sinhA = add_Z(sinhA,sA);  
   }
   
  return (sinhA);	
}
/* ------------------------------------ */
nb_Z cosh_Z (
nb_Z A
)
{	
nb_Z   A_n;
nb_Z    sA;
nb_Z coshA = i_Z(0,0);

int n = 12; 
int i; 
	
  for(i = 0.; i <= n; i+=2.)
   {
       A_n = pow_Z(i, A);         
        sA = smul_Z(1./factorial(i), A_n);
     coshA = add_Z(coshA,sA);  
   }
   
  return (coshA);	
}
/* ------------------------------------ */
/* ------------------------------------ */
nb_Z exp_Z(
nb_Z A
)
{	
nb_Z   A_n;
nb_Z   sA ;
nb_Z expA = i_Z(0,0);

int n = 10; 
int i; 
	
  for(i=0.; i <= n; i++)
   {
     A_n  = pow_Z(i, A);         
     sA   = smul_Z(1./factorial(i), A_n);
     expA = add_Z(expA,sA);  
   }
   
  return (expA);	
}
/* ------------------------------------ */
/* ------------------------------------ */