double SF(double x, int n) {  //n为正整数
           double p=1,s=1;
           int i;
           for(i=1;i<=n;i ) {
               p*=x;
               s =p;
           }
           return s;
       }