Mathc complexes/062

Installer ce fichier dans votre répertoire de travail.

	d.h

/* ------------------------------------ */
/*  Save as :   d.h                     */
/* ------------------------------------ */
void i_A_b_with_XY_mZ(
double **XY,
double **A,
double **b
)
{
int r     = R1;
int c     = C1;
int power = 0;

  for(r=R1; r<XY[R_SIZE][C0]; r++)
    {
     power=rsize_Z(XY)-R1;                  
     
     for(c=C1; c<A[C_SIZE][C0]; c++,c++)
     
        A[r][c]=pow(XY[r][C1],power--);

     b[r][C1]=XY[r][C3];
    }
}
/* --------------------------------- */
void p_eq_poly_mZ(
double **Ab
)
{
int      r = R1;
int  power = rsize_Z(Ab)-R1;           
int     cL = csize_Z(Ab)*C2-C1;

  printf("  y = ");
  
   for(r=R1;r<Ab[R_SIZE][C0];r++)
   
     if(Ab[r][cL])
     {
            if(!power)   printf(" %+.3f",     Ab[r][cL]);
       else if(power==1){printf(" %+.3fx",    Ab[r][cL]);power--;}
       else              printf(" %+.3fx**%d",Ab[r][cL], power--);
      }
      else power--;
      
  printf("\n\n\n");
}
/* --------------------------------- */
void verify_X_mZ(
double **Ab,
double x
)
{
int     r = R1;
int power = rsize_Z(Ab)-R1;
int    cL = csize_Z(Ab)*C2-C1;
double  y = 0.;

    for(;r<Ab[R_SIZE][C0];r++)
    
       y+= Ab[r][cL]*pow(x,power--);

    printf(" With x =  %+.3f,       y = %+.3f \n",x,y);
}
/* --------------------------------- */
/* --------------------------------- */

Déclaration des fichiers h.