xppoly.h

来自「基本矩阵运算 : + - *, power, transpose, trace,」· C头文件 代码 · 共 57 行

/* xppoly.h                     freeware                   [e-mail]  */

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */


/* -------------------------------------------------------------------------- */
void p_polyI(
PmF mpoly)
{
   int i;
double coeff;
   int p0 = 1;
   int p = mpoly->rows;

printf("   ");

for (   i = 0 ; i < mpoly->rows    ; i++)
{
 coeff = *(mpoly->pb+i *mpoly->cols);

 --p;

 if(coeff)
 {
  p0 = 0;
  if(p)
  {
     if(     coeff == 1)
       {
        printf("+ x**%d  ",p);
       }
     else if(coeff == -1)
       {
        printf("- x**%d  ",p);
       }
     else
       {
        printf("%+.0fx**%d  ",coeff,p);
       }
  }
  else
  {
     printf("%+.0f  ",coeff);
  }
 }

}

   if(p0){printf("   0");}

   printf("\n");
}