zernike代码

#include  < stdio.h >
#include  < string .h >
#include  < math.h >
#include  " cv.h "
#include  " highgui.h "  

#define  MaxP(x,y) (float)(x>y?x:y)
#define  MinP(x,y) (float)(x<y?x:y)
#define  PI  3.14
#define   MAXFACT    50
unsigned  int  fact_array[MAXFACT];


IplImage *     oriImg; 
IplImage *     grayImg;

unsigned  int   fact( int   n)
{        
     if (fact_array[n] != 0 )
         return  fact_array[n];
    fact_array[ 0 ]  =   1 ;
    fact_array[ 1 ]  =   1 ;
     int   index     =  n;
    
    fact_array[n]  =  n;
    index -- ;

  while (fact_array[index] == 0 ){
     fact_array[n]  =  fact_array[n]  *  index;
     index -- ;         
     }
    fact_array[n]  =   fact_array[n] * fact_array[index];
   
    return  fact_array[n];

}

double  ZernikePolynomial( int  n, int  m, double  r)
{
        
     if (n > 50 )
    {
    printf( " I can only deal with factorial lower than 50  " );
     return   0 ;}

     int  s;
     double  sum