function FFit=PolyFun(x,y,RotAng,term)
   
    p=sqrt(x.^2+y.^2);
    A=atan2(x,y)+RotAng;
    leng=length(x);
    

FFit(:,1)=zeros(leng,1)+1;
FFit(:,2)=4^(1/2)*(p).*cos(A);
FFit(:,3)=4^(1/2)*(p).*sin(A);
FFit(:,4)=3^(1/2)*(2*p.^2 - 1);
FFit(:,5)=6^(1/2)*(p.^2).*sin(2*A);
FFit(:,6)=6^(1/2)*(p.^2).* cos(2*A);
FFit(:,7)=8^(1/2)*(3*p.^3 - 2*p).* sin(A);
FFit(:,8)=8^(1/2)*(3*p.^3 - 2*p).* cos(A);
FFit(:,9)=8^(1/2)*(p.^3).* sin(3*A);
FFit(:,10)=8^(1/2)*(p.^3).* cos(3*A);
FFit(:,11)=5^(1/2)*(6*p.^4 - 6*p.^2 + 1);
FFit(:,12)=10^(1/2)*(4*p.^4 - 3*p.^2).*cos(2*A);
FFit(:,13)=10^(1/2)*(4*p.^4 - 3*p.^2).*sin(2*A);
FFit(:,14)=10^(1/2)*(p.^4).*cos(4*A);
FFit(:,15)=10^(1/2)*(p.^4).*sin(4*A);
FFit(:,16)=12^(1/2)*(10*p.^5 - 12*p.^3 + 3*p).*cos(A);
FFit(:,17)=12^(1/2)*(10*p.^5 - 12*p.^3 + 3*p).*sin(A);
FFit(:,18)=12^(1/2)*(5*p.^5 - 4*p.^3).*cos(3*A);
FFit(:,19)=12^(1/2)*(5*p.^5 - 4*p.^3).*sin(3*A);
FFit(:,20)=12^(1/2)*(p.^5).*cos(5*A);
FFit(:,21)=12^(1/2)*(p.^5).*sin(5*A);
FFit(:,22)=7^(1/2)*(20*p.^6 - 30*p.^4 + 12*p.^2 - 1);
FFit(:,23)=14^(1/2)*(15*p.^6 - 20*p.^4 + 6*p.^2).*sin(2*A);
FFit(:,24)=14^(1/2)*(15*p.^6 - 20*p.^4 + 6*p.^2).*cos(2*A);
FFit(:,25)=14^(1/2)*(6*p.^6 - 5*p.^4).*sin(4*A);
FFit(:,26)=14^(1/2)*(6*p.^6 - 5*p.^4).*cos(4*A);
FFit(:,27)=14^(1/2)*(p.^6).*sin(6*A);
FFit(:,28)=14^(1/2)*(p.^6).*cos(6*A);
FFit(:,29)=16^(1/2)*(35*p.^7 - 60*p.^5 + 30*p.^3 - 4*p).*sin(A);
FFit(:,30)=16^(1/2)*(35*p.^7 - 60*p.^5 + 30*p.^3 - 4*p).*cos(A);
FFit(:,31)=16^(1/2)*(21*p.^7 - 30*p.^5 + 10*p.^3).*sin(3*A);
FFit(:,32)=16^(1/2)*(21*p.^7 - 30*p.^5 + 10*p.^3).*cos(3*A);
FFit(:,33)=16^(1/2)*(7*p.^7 - 6*p.^5).*sin(5*A);
FFit(:,34)=16^(1/2)*(7*p.^7 - 6*p.^5).*cos(5*A);
FFit(:,35)=16^(1/2)*(p.^7).*sin(7*A);
FFit(:,36)=16^(1/2)*(p.^7).*cos(7*A);
FFit(:,37)= 9^(1/2)*(70*p.^8 - 140*p.^6 + 90*p.^4 - 20*p.^2 + 1);



% 
% FFit(:,38)=18^(1/2)*(56*p.^8 - 105*p.^6 + 60*p.^4 - 10*p.^2).* cos(2*A);
% FFit(:,39)=18^(1/2)*(56*p.^8 - 105*p.^6 + 60*p.^4 - 10*p.^2).* sin(2*A);
% FFit(:,40)=18^(1/2)*(28*p.^8 - 42*p.^6 + 15*p.^4).* cos(4*A);
% FFit(:,41)=18^(1/2)*(28*p.^8 - 42*p.^6 + 15*p.^4).* sin(4*A);
% FFit(:,42)=18^(1/2)*(8*p.^8 - 7*p.^6).* cos(6*A);
% FFit(:,43)=18^(1/2)*(8*p.^8 - 7*p.^6).* sin(6*A);
% FFit(:,44)=18^(1/2)*(p.^8).* cos(8*A);
% FFit(:,45)=18^(1/2)*(p.^8).* sin(8*A);
% FFit(:,46)=20^(1/2)*(126*p.^9 - 280*p.^7 + 210*p.^5 - 60*p.^3 + 5*p).* cos(A);
% FFit(:,47)=20^(1/2)*(126*p.^9 - 280*p.^7 + 210*p.^5 - 60*p.^3 + 5*p).* sin(A);
% FFit(:,48)=20^(1/2)*(84*p.^9 - 168*p.^7 + 105*p.^5 - 20*p.^3).* cos(3*A);
% FFit(:,49)=20^(1/2)*(84*p.^9 - 168*p.^7 + 105*p.^5 - 20*p.^3).* sin(3*A);
% FFit(:,50)= 20^(1/2)*(36*p.^9 - 56*p.^7 + 21*p.^5).* cos(5*A);
% FFit(:,51)= 20^(1/2)*(36*p.^9 - 56*p.^7 + 21*p.^5).* sin(5*A);
% FFit(:,52)=20^(1/2)*(9*p.^9 - 8*p.^7).* cos(7*A);
% FFit(:,53)=20^(1/2)*(9*p.^9 - 8*p.^7).* sin(7*A);
% FFit(:,54)=20^(1/2)*(p.^9).* cos(9*A);
% FFit(:,55)=20^(1/2)*(p.^9).* sin(9*A);
% FFit(:,56)=11^(1/2)*(252*p.^10 - 630*p.^8 + 560*p.^6 - 210*p.^4 + 30*p.^2 - 1);
% FFit(:,57)= 22^(1/2)*(210*p.^10 - 504*p.^8 + 420*p.^6 - 140*p.^4 + 15*p.^2).* sin(2*A);
% FFit(:,58)=22^(1/2)*(210*p.^10 - 504*p.^8 + 420*p.^6 - 140*p.^4 + 15*p.^2).* cos(2*A);
% FFit(:,59)=22^(1/2)*(120*p.^10 - 252*p.^8 + 168*p.^6 - 35*p.^4).* sin(4*A);
% FFit(:,60)=22^(1/2)*(120*p.^10 - 252*p.^8 + 168*p.^6 - 35*p.^4).* cos(4*A);
% FFit(:,61)=22^(1/2)*(45*p.^10 - 72*p.^8 + 28*p.^6).* sin(6*A);
% FFit(:,62)=22^(1/2)*(45*p.^10 - 72*p.^8 + 28*p.^6).* cos(6*A);
% FFit(:,63)=22^(1/2)*(10*p.^10 - 9*p.^8).* sin(8*A);
% FFit(:,64)=22^(1/2)*(10*p.^10 - 9*p.^8).* cos(8*A);
% FFit(:,65)=22^(1/2)*(p.^10).* sin(10*A);
% FFit(:,66)=22^(1/2)*(p.^10).* cos(10*A);
% FFit(:,67)=24^(1/2)*(462*p.^11 - 1260*p.^9 + 1260*p.^7 - 560*p.^5 + 105*p.^3 - 6*p).* sin(A);
% FFit(:,68)=24^(1/2)*(462*p.^11 - 1260*p.^9 + 1260*p.^7 - 560*p.^5 + 105*p.^3 - 6*p).* cos(A);
% FFit(:,69)=24^(1/2)*(330*p.^11 - 840*p.^9 + 756*p.^7 - 280*p.^5 + 35*p.^3).* sin(3*A);
% FFit(:,70)=24^(1/2)*(330*p.^11 - 840*p.^9 + 756*p.^7 - 280*p.^5 + 35*p.^3).* cos(3*A);
% FFit(:,71)=24^(1/2)*(165*p.^11 - 360*p.^9 + 252*p.^7 - 56*p.^5).* sin(5*A);
% FFit(:,72)=24^(1/2)*(165*p.^11 - 360*p.^9 + 252*p.^7 - 56*p.^5).* cos(5*A);
% FFit(:,73)=24^(1/2)*(55*p.^11 - 90*p.^9 + 36*p.^7).* sin(7*A);
% FFit(:,74)= 24^(1/2)*(55*p.^11 - 90*p.^9 + 36*p.^7).* cos(7*A);
% FFit(:,75)=24^(1/2)*(11*p.^11 - 10*p.^9).* sin(9*A);
% FFit(:,76)=24^(1/2)*(11*p.^11 - 10*p.^9).* cos(9*A);
% FFit(:,77)=24^(1/2)*(p.^11).* sin(11*A);
% FFit(:,78)=24^(1/2)*(p.^11).* cos(11*A);
% FFit(:,79)=13^(1/2)*(924*p.^12 - 2772*p.^10 + 3150*p.^8 - 1680*p.^6 + 420*p.^4 - 42*p.^2 + 1);
% FFit(:,80)=26^(1/2)*(792*p.^12 - 2310*p.^10 + 2520*p.^8 - 1260*p.^6 + 280*p.^4 - 21*p.^2).* cos(2*A);
% % 
% % 
% FFit(:,81)=26^(1/2)*(792*p.^12 - 2310*p.^10 + 2520*p.^8 - 1260*p.^6 + 280*p.^4 - 21*p.^2).* sin(2*A);
% FFit(:,82)=26^(1/2)*(495*p.^12 - 1320*p.^10 + 1260*p.^8 - 504*p.^6 + 70*p.^4).* cos(4*A);
% FFit(:,83)=26^(1/2)*(495*p.^12 - 1320*p.^10 + 1260*p.^8 - 504*p.^6 + 70*p.^4).* sin(4*A);
% FFit(:,84)=26^(1/2)*(220*p.^12 - 495*p.^10 + 360*p.^8 - 84*p.^6).* cos(6*A);
% FFit(:,85)=26^(1/2)*(220*p.^12 - 495*p.^10 + 360*p.^8 - 84*p.^6).* sin(6*A);
% FFit(:,86)=26^(1/2)*(66*p.^12 - 110*p.^10 + 45*p.^8).* cos(8*A);
% FFit(:,87)=26^(1/2)*(66*p.^12 - 110*p.^10 + 45*p.^8).* sin(8*A);
% FFit(:,88)=26^(1/2)*(12*p.^12 - 11*p.^10).* cos(10*A);
% FFit(:,89)=26^(1/2)*(12*p.^12 - 11*p.^10).* sin(10*A);
% FFit(:,90)=26^(1/2)*(p.^12).* cos(12*A);
% FFit(:,91)=26^(1/2)*(p.^12).* sin(12*A);
% FFit(:,92)=28^(1/2)*(1716*p.^13 - 5544*p.^11 + 6930*p.^9 - 4200*p.^7 + 1260*p.^5 - 168*p.^3 + 7*p).* cos(A);
% FFit(:,93)=28^(1/2)*(1716*p.^13 - 5544*p.^11 + 6930*p.^9 - 4200*p.^7 + 1260*p.^5 - 168*p.^3 + 7*p).* sin(A);
% FFit(:,94)=28^(1/2)*(1287*p.^13 - 3960*p.^11 + 4620*p.^9 - 2520*p.^7 + 630*p.^5 - 56*p.^3).* cos(3*A);
% FFit(:,95)=28^(1/2)*(1287*p.^13 - 3960*p.^11 + 4620*p.^9 - 2520*p.^7 + 630*p.^5 - 56*p.^3).* sin(3*A);
% FFit(:,96)=28^(1/2)*(715*p.^13 - 1980*p.^11 + 1980*p.^9 - 840*p.^7 + 126*p.^5).* cos(5*A);
% FFit(:,97)=28^(1/2)*(715*p.^13 - 1980*p.^11 + 1980*p.^9 - 840*p.^7 + 126*p.^5).* sin(5*A);
% FFit(:,98)=28^(1/2)*(286*p.^13 - 660*p.^11 + 495*p.^9 - 120*p.^7).* cos(7*A);
% FFit(:,99)=28^(1/2)*(286*p.^13 - 660*p.^11 + 495*p.^9 - 120*p.^7).* sin(7*A);
% FFit(:,100)=28^(1/2)*(78*p.^13 - 132*p.^11 + 55*p.^9).* cos(9*A);
% FFit(:,101)=28^(1/2)*(78*p.^13 - 132*p.^11 + 55*p.^9).* sin(9*A);
% FFit(:,102)=28^(1/2)*(13*p.^13 - 12*p.^11).* cos(11*A);
% FFit(:,103)=28^(1/2)*(13*p.^13 - 12*p.^11).* sin(11*A);
% FFit(:,104)=28^(1/2)*(p.^13).* cos(13*A);
% FFit(:,105)=28^(1/2)*(p.^13).* sin(13*A);
% FFit(:,106)=15^(1/2)*(3432*p.^14 - 12012*p.^12 + 16632*p.^10 - 11550*p.^8 + 4200*p.^6 - 756*p.^4 + 56*p.^2 - 1);
% FFit(:,107)=30^(1/2)*(3003*p.^14 - 10296*p.^12 + 13860*p.^10 - 9240*p.^8 + 3150*p.^6 - 504*p.^4 + 28*p.^2).* sin(2*A);
% FFit(:,108)=30^(1/2)*(3003*p.^14 - 10296*p.^12 + 13860*p.^10 - 9240*p.^8 + 3150*p.^6 - 504*p.^4 + 28*p.^2).* cos(2*A);
% FFit(:,109)=30^(1/2)*(2002*p.^14 - 6435*p.^12 + 7920*p.^10 - 4620*p.^8 + 1260*p.^6 - 126*p.^4).* sin(4*A);
% FFit(:,110)=30^(1/2)*(2002*p.^14 - 6435*p.^12 + 7920*p.^10 - 4620*p.^8 + 1260*p.^6 - 126*p.^4).* cos(4*A);
% FFit(:,111)=30^(1/2)*(1001*p.^14 - 2860*p.^12 + 2970*p.^10 - 1320*p.^8 + 210*p.^6).* sin(6*A);
% FFit(:,112)=30^(1/2)*(1001*p.^14 - 2860*p.^12 + 2970*p.^10 - 1320*p.^8 + 210*p.^6).* cos(6*A);
% FFit(:,113)=30^(1/2)*(364*p.^14 - 858*p.^12 + 660*p.^10 - 165*p.^8).* sin(8*A);
% FFit(:,114)=30^(1/2)*(364*p.^14 - 858*p.^12 + 660*p.^10 - 165*p.^8).* cos(8*A);
% FFit(:,115)=30^(1/2)*(91*p.^14 - 156*p.^12 + 66*p.^10).* sin(10*A);
% FFit(:,116)=30^(1/2)*(91*p.^14 - 156*p.^12 + 66*p.^10).* cos(10*A);
% FFit(:,117)=30^(1/2)*(14*p.^14 - 13*p.^12).* sin(12*A);
% FFit(:,118)=30^(1/2)*(14*p.^14 - 13*p.^12).* cos(12*A);
% FFit(:,119)=30^(1/2)*(p.^14).* sin(14*A);
% FFit(:,120)=30^(1/2)*(p.^14).* cos(14*A);
% FFit(:,121)=32^(1/2)*(6435*p.^15 - 24024*p.^13 + 36036*p.^11 - 27720*p.^9 + 11550*p.^7 - 2520*p.^5 + 252*p.^3 - 8*p).* sin(A);
% FFit(:,122)=32^(1/2)*(6435*p.^15 - 24024*p.^13 + 36036*p.^11 - 27720*p.^9 + 11550*p.^7 - 2520*p.^5 + 252*p.^3 - 8*p).* cos(A);
% FFit(:,123)=32^(1/2)*(5005*p.^15 - 18018*p.^13 + 25740*p.^11 - 18480*p.^9 + 6930*p.^7 - 1260*p.^5 + 84*p.^3).* sin(3*A);
% FFit(:,124)=32^(1/2)*(5005*p.^15 - 18018*p.^13 + 25740*p.^11 - 18480*p.^9 + 6930*p.^7 - 1260*p.^5 + 84*p.^3).* cos(3*A);
% FFit(:,125)=32^(1/2)*(3003*p.^15 - 10010*p.^13 + 12870*p.^11 - 7920*p.^9 + 2310*p.^7 - 252*p.^5).* sin(5*A);
% FFit(:,126)=32^(1/2)*(3003*p.^15 - 10010*p.^13 + 12870*p.^11 - 7920*p.^9 + 2310*p.^7 - 252*p.^5).* cos(5*A);
% FFit(:,127)=32^(1/2)*(1365*p.^15 - 4004*p.^13 + 4290*p.^11 - 1980*p.^9 + 330*p.^7).* sin(7*A);
% FFit(:,128)=32^(1/2)*(1365*p.^15 - 4004*p.^13 + 4290*p.^11 - 1980*p.^9 + 330*p.^7).* cos(7*A);
% FFit(:,129)=32^(1/2)*(455*p.^15 - 1092*p.^13 + 858*p.^11 - 220*p.^9).* sin(9*A);
% FFit(:,130)=32^(1/2)*(455*p.^15 - 1092*p.^13 + 858*p.^11 - 220*p.^9).* cos(9*A);
% FFit(:,131)=32^(1/2)*(105*p.^15 - 182*p.^13 + 78*p.^11).* sin(11*A);
% FFit(:,132)=32^(1/2)*(105*p.^15 - 182*p.^13 + 78*p.^11).* cos(11*A);
% FFit(:,133)=32^(1/2)*(15*p.^15 - 14*p.^13).* sin(13*A);
% FFit(:,134)=32^(1/2)*(15*p.^15 - 14*p.^13).* cos(13*A);
% FFit(:,135)=32^(1/2)*(p.^15).* sin(15*A);
% FFit(:,136)=32^(1/2)*(p.^15).* cos(15*A);
% FFit(:,137)=17^(1/2)*(12870*p.^16 - 51480*p.^14 + 84084*p.^12 - 72072*p.^10 + 34650*p.^8 - 9240*p.^6 + 1260*p.^4 - 72*p.^2 + 1);
% FFit(:,138)=34^(1/2)*(11440*p.^16 - 45045*p.^14 + 72072*p.^12 - 60060*p.^10 + 27720*p.^8 - 6930*p.^6 + 840*p.^4 - 36*p.^2).* cos(2*A);
% FFit(:,139)=34^(1/2)*(11440*p.^16 - 45045*p.^14 + 72072*p.^12 - 60060*p.^10 + 27720*p.^8 - 6930*p.^6 + 840*p.^4 - 36*p.^2).* sin(2*A);
% FFit(:,140)=34^(1/2)*(8008*p.^16 - 30030*p.^14 + 45045*p.^12 - 34320*p.^10 + 13860*p.^8 - 2772*p.^6 + 210*p.^4).* cos(4*A);
% FFit(:,141)=34^(1/2)*(8008*p.^16 - 30030*p.^14 + 45045*p.^12 - 34320*p.^10 + 13860*p.^8 - 2772*p.^6 + 210*p.^4).* sin(4*A);
% FFit(:,142)=34^(1/2)*(4368*p.^16 - 15015*p.^14 + 20020*p.^12 - 12870*p.^10 + 3960*p.^8 - 462*p.^6).* cos(6*A);
% FFit(:,143)=34^(1/2)*(4368*p.^16 - 15015*p.^14 + 20020*p.^12 - 12870*p.^10 + 3960*p.^8 - 462*p.^6).* sin(6*A);
% FFit(:,144)=34^(1/2)*(1820*p.^16 - 5460*p.^14 + 6006*p.^12 - 2860*p.^10 + 495*p.^8).* cos(8*A);
% FFit(:,145)=34^(1/2)*(1820*p.^16 - 5460*p.^14 + 6006*p.^12 - 2860*p.^10 + 495*p.^8).* sin(8*A);
% FFit(:,146)=34^(1/2)*(560*p.^16 - 1365*p.^14 + 1092*p.^12 - 286*p.^10).* cos(10*A);
% FFit(:,147)=34^(1/2)*(560*p.^16 - 1365*p.^14 + 1092*p.^12 - 286*p.^10).* sin(10*A);
% FFit(:,148)=34^(1/2)*(120*p.^16 - 210*p.^14 + 91*p.^12).* cos(12*A);
% FFit(:,149)=34^(1/2)*(120*p.^16 - 210*p.^14 + 91*p.^12).* sin(12*A);
% FFit(:,150)=34^(1/2)*(16*p.^16 - 15*p.^14).* cos(14*A);
% FFit(:,151)=34^(1/2)*(16*p.^16 - 15*p.^14).* sin(14*A);
% FFit(:,152)=34^(1/2)*(p.^16).* cos(16*A);
% FFit(:,153)=34^(1/2)*(p.^16).* sin(16*A);
% FFit(:,154)=36^(1/2)*(24310*p.^17 - 102960*p.^15 + 180180*p.^13 - 168168*p.^11 + 90090*p.^9 - 27720*p.^7 + 4620*p.^5 - 360*p.^3 + 9*p).* cos(A);
% FFit(:,155)=36^(1/2)*(24310*p.^17 - 102960*p.^15 + 180180*p.^13 - 168168*p.^11 + 90090*p.^9 - 27720*p.^7 + 4620*p.^5 - 360*p.^3 + 9*p).* sin(A);
% FFit(:,156)=36^(1/2)*(19448*p.^17 - 80080*p.^15 + 135135*p.^13 - 120120*p.^11 + 60060*p.^9 - 16632*p.^7 + 2310*p.^5 - 120*p.^3).* cos(3*A);
% FFit(:,157)=36^(1/2)*(19448*p.^17 - 80080*p.^15 + 135135*p.^13 - 120120*p.^11 + 60060*p.^9 - 16632*p.^7 + 2310*p.^5 - 120*p.^3).* sin(3*A);
% FFit(:,158)=36^(1/2)*(12376*p.^17 - 48048*p.^15 + 75075*p.^13 - 60060*p.^11 + 25740*p.^9 - 5544*p.^7 + 462*p.^5).* cos(5*A);
% FFit(:,159)=36^(1/2)*(12376*p.^17 - 48048*p.^15 + 75075*p.^13 - 60060*p.^11 + 25740*p.^9 - 5544*p.^7 + 462*p.^5).* sin(5*A);
% FFit(:,160)=36^(1/2)*(6188*p.^17 - 21840*p.^15 + 30030*p.^13 - 20020*p.^11 + 6435*p.^9 - 792*p.^7).* cos(7*A);
% FFit(:,161)=36^(1/2)*(6188*p.^17 - 21840*p.^15 + 30030*p.^13 - 20020*p.^11 + 6435*p.^9 - 792*p.^7).* sin(7*A);
% FFit(:,162)=36^(1/2)*(2380*p.^17 - 7280*p.^15 + 8190*p.^13 - 4004*p.^11 + 715*p.^9).* cos(9*A);
% FFit(:,163)=36^(1/2)*(2380*p.^17 - 7280*p.^15 + 8190*p.^13 - 4004*p.^11 + 715*p.^9).* sin(9*A);
% FFit(:,164)=36^(1/2)*(680*p.^17 - 1680*p.^15 + 1365*p.^13 - 364*p.^11).* cos(11*A);
% FFit(:,165)=36^(1/2)*(680*p.^17 - 1680*p.^15 + 1365*p.^13 - 364*p.^11).* sin(11*A);
% FFit(:,166)=36^(1/2)*(136*p.^17 - 240*p.^15 + 105*p.^13).* cos(13*A);
% FFit(:,167)=36^(1/2)*(136*p.^17 - 240*p.^15 + 105*p.^13).* sin(13*A);
% FFit(:,168)=36^(1/2)*(17*p.^17 - 16*p.^15).* cos(15*A);
% FFit(:,169)=36^(1/2)*(17*p.^17 - 16*p.^15).* sin(15*A);
% FFit(:,170)=36^(1/2)*(p.^17).* cos(17*A);
% FFit(:,171)=36^(1/2)*(p.^17).* sin(17*A);
% FFit(:,172)=19^(1/2)*(48620*p.^18 - 218790*p.^16 + 411840*p.^14 - 420420*p.^12 + 252252*p.^10 - 90090*p.^8 + 18480*p.^6 - 1980*p.^4 + 90*p.^2 - 1);
% FFit(:,173)=38^(1/2)*(43758*p.^18 - 194480*p.^16 + 360360*p.^14 - 360360*p.^12 + 210210*p.^10 - 72072*p.^8 + 13860*p.^6 - 1320*p.^4 + 45*p.^2).* sin(2*A);
% FFit(:,174)=38^(1/2)*(43758*p.^18 - 194480*p.^16 + 360360*p.^14 - 360360*p.^12 + 210210*p.^10 - 72072*p.^8 + 13860*p.^6 - 1320*p.^4 + 45*p.^2).* cos(2*A);
% FFit(:,175)=38^(1/2)*(31824*p.^18 - 136136*p.^16 + 240240*p.^14 - 225225*p.^12 + 120120*p.^10 - 36036*p.^8 + 5544*p.^6 - 330*p.^4).* sin(4*A);
% FFit(:,176)=38^(1/2)*(31824*p.^18 - 136136*p.^16 + 240240*p.^14 - 225225*p.^12 + 120120*p.^10 - 36036*p.^8 + 5544*p.^6 - 330*p.^4).* cos(4*A);
% FFit(:,177)=38^(1/2)*(18564*p.^18 - 74256*p.^16 + 120120*p.^14 - 100100*p.^12 + 45045*p.^10 - 10296*p.^8 + 924*p.^6).* sin(6*A);
% FFit(:,178)=38^(1/2)*(18564*p.^18 - 74256*p.^16 + 120120*p.^14 - 100100*p.^12 + 45045*p.^10 - 10296*p.^8 + 924*p.^6).* cos(6*A);
% FFit(:,179)=38^(1/2)*(8568*p.^18 - 30940*p.^16 + 43680*p.^14 - 30030*p.^12 + 10010*p.^10 - 1287*p.^8).* sin(8*A);
% FFit(:,180)=38^(1/2)*(8568*p.^18 - 30940*p.^16 + 43680*p.^14 - 30030*p.^12 + 10010*p.^10 - 1287*p.^8).* cos(8*A);
% FFit(:,181)=38^(1/2)*(3060*p.^18 - 9520*p.^16 + 10920*p.^14 - 5460*p.^12 + 1001*p.^10).* sin(10*A);
% FFit(:,182)=38^(1/2)*(3060*p.^18 - 9520*p.^16 + 10920*p.^14 - 5460*p.^12 + 1001*p.^10).* cos(10*A);
% FFit(:,183)=38^(1/2)*(816*p.^18 - 2040*p.^16 + 1680*p.^14 - 455*p.^12).* sin(12*A);
% FFit(:,184)=38^(1/2)*(816*p.^18 - 2040*p.^16 + 1680*p.^14 - 455*p.^12).* cos(12*A);
% FFit(:,185)=38^(1/2)*(153*p.^18 - 272*p.^16 + 120*p.^14).* sin(14*A);
% FFit(:,186)=38^(1/2)*(153*p.^18 - 272*p.^16 + 120*p.^14).* cos(14*A);
% FFit(:,187)=38^(1/2)*(18*p.^18 - 17*p.^16).* sin(16*A);
% FFit(:,188)=38^(1/2)*(18*p.^18 - 17*p.^16).* cos(16*A);
% FFit(:,189)=38^(1/2)*(p.^18).* sin(18*A);
% FFit(:,190)=38^(1/2)*(p.^18).* cos(18*A);
% FFit(:,191)=40^(1/2)*(92378*p.^19 - 437580*p.^17 + 875160*p.^15 - 960960*p.^13 + 630630*p.^11 - 252252*p.^9 + 60060*p.^7 - 7920*p.^5 + 495*p.^3 - 10*p).* sin(A);
% FFit(:,192)=40^(1/2)*(92378*p.^19 - 437580*p.^17 + 875160*p.^15 - 960960*p.^13 + 630630*p.^11 - 252252*p.^9 + 60060*p.^7 - 7920*p.^5 + 495*p.^3 - 10*p).* cos(A);
% FFit(:,193)=40^(1/2)*(75582*p.^19 - 350064*p.^17 + 680680*p.^15 - 720720*p.^13 + 450450*p.^11 - 168168*p.^9 + 36036*p.^7 - 3960*p.^5 + 165*p.^3).* sin(3*A);
% FFit(:,194)=40^(1/2)*(75582*p.^19 - 350064*p.^17 + 680680*p.^15 - 720720*p.^13 + 450450*p.^11 - 168168*p.^9 + 36036*p.^7 - 3960*p.^5 + 165*p.^3).* cos(3*A);
% FFit(:,195)=40^(1/2)*(50388*p.^19 - 222768*p.^17 + 408408*p.^15 - 400400*p.^13 + 225225*p.^11 - 72072*p.^9 + 12012*p.^7 - 792*p.^5).* sin(5*A);
% FFit(:,196)=40^(1/2)*(50388*p.^19 - 222768*p.^17 + 408408*p.^15 - 400400*p.^13 + 225225*p.^11 - 72072*p.^9 + 12012*p.^7 - 792*p.^5).* cos(5*A);
% FFit(:,197)=40^(1/2)*(27132*p.^19 - 111384*p.^17 + 185640*p.^15 - 160160*p.^13 + 75075*p.^11 - 18018*p.^9 + 1716*p.^7).* sin(7*A);
% FFit(:,198)=40^(1/2)*(27132*p.^19 - 111384*p.^17 + 185640*p.^15 - 160160*p.^13 + 75075*p.^11 - 18018*p.^9 + 1716*p.^7).* cos(7*A);
% FFit(:,199)=40^(1/2)*(11628*p.^19 - 42840*p.^17 + 61880*p.^15 - 43680*p.^13 + 15015*p.^11 - 2002*p.^9).* sin(9*A);
% FFit(:,200)=40^(1/2)*(11628*p.^19 - 42840*p.^17 + 61880*p.^15 - 43680*p.^13 + 15015*p.^11 - 2002*p.^9).* cos(9*A);
% FFit(:,201)=40^(1/2)*(3876*p.^19 - 12240*p.^17 + 14280*p.^15 - 7280*p.^13 + 1365*p.^11).* sin(11*A);
% FFit(:,202)=40^(1/2)*(3876*p.^19 - 12240*p.^17 + 14280*p.^15 - 7280*p.^13 + 1365*p.^11).* cos(11*A);
% FFit(:,203)=40^(1/2)*(969*p.^19 - 2448*p.^17 + 2040*p.^15 - 560*p.^13).* sin(13*A);
% FFit(:,204)=40^(1/2)*(969*p.^19 - 2448*p.^17 + 2040*p.^15 - 560*p.^13).* cos(13*A);
% FFit(:,205)=40^(1/2)*(171*p.^19 - 306*p.^17 + 136*p.^15).* sin(15*A);
% FFit(:,206)=40^(1/2)*(171*p.^19 - 306*p.^17 + 136*p.^15).* cos(15*A);
% FFit(:,207)=40^(1/2)*(19*p.^19 - 18*p.^17).* sin(17*A);
% FFit(:,208)=40^(1/2)*(19*p.^19 - 18*p.^17).* cos(17*A);
% FFit(:,209)=40^(1/2)*(p.^19).* sin(19*A);
% FFit(:,210)=40^(1/2)*(p.^19).* cos(19*A);
% FFit(:,211)=21^(1/2)*(184756*p.^20 - 923780*p.^18 + 1969110*p.^16 - 2333760*p.^14 + 1681680*p.^12 - 756756*p.^10 + 210210*p.^8 - 34320*p.^6 + 2970*p.^4 - 110*p.^2 + 1);
% FFit(:,212)=42^(1/2)*(167960*p.^20 - 831402*p.^18 + 1750320*p.^16 - 2042040*p.^14 + 1441440*p.^12 - 630630*p.^10 + 168168*p.^8 - 25740*p.^6 + 1980*p.^4 - 55*p.^2).* cos(2*A);
% FFit(:,213)=42^(1/2)*(167960*p.^20 - 831402*p.^18 + 1750320*p.^16 - 2042040*p.^14 + 1441440*p.^12 - 630630*p.^10 + 168168*p.^8 - 25740*p.^6 + 1980*p.^4 - 55*p.^2).* sin(2*A);
% FFit(:,214)=42^(1/2)*(125970*p.^20 - 604656*p.^18 + 1225224*p.^16 - 1361360*p.^14 + 900900*p.^12 - 360360*p.^10 + 84084*p.^8 - 10296*p.^6 + 495*p.^4).* cos(4*A);
% FFit(:,215)=42^(1/2)*(125970*p.^20 - 604656*p.^18 + 1225224*p.^16 - 1361360*p.^14 + 900900*p.^12 - 360360*p.^10 + 84084*p.^8 - 10296*p.^6 + 495*p.^4).* sin(4*A);
% FFit(:,216)=42^(1/2)*(77520*p.^20 - 352716*p.^18 + 668304*p.^16 - 680680*p.^14 + 400400*p.^12 - 135135*p.^10 + 24024*p.^8 - 1716*p.^6).* cos(6*A);
% FFit(:,217)=42^(1/2)*(77520*p.^20 - 352716*p.^18 + 668304*p.^16 - 680680*p.^14 + 400400*p.^12 - 135135*p.^10 + 24024*p.^8 - 1716*p.^6).* sin(6*A);
% FFit(:,218)=42^(1/2)*(38760*p.^20 - 162792*p.^18 + 278460*p.^16 - 247520*p.^14 + 120120*p.^12 - 30030*p.^10 + 3003*p.^8).* cos(8*A);
% FFit(:,219)=42^(1/2)*(38760*p.^20 - 162792*p.^18 + 278460*p.^16 - 247520*p.^14 + 120120*p.^12 - 30030*p.^10 + 3003*p.^8).* sin(8*A);
% FFit(:,220)=42^(1/2)*(15504*p.^20 - 58140*p.^18 + 85680*p.^16 - 61880*p.^14 + 21840*p.^12 - 3003*p.^10).* cos(10*A);
% FFit(:,221)=42^(1/2)*(15504*p.^20 - 58140*p.^18 + 85680*p.^16 - 61880*p.^14 + 21840*p.^12 - 3003*p.^10).* sin(10*A);
% % FFit(:,222)=42^(1/2)*(4845*p.^20 - 15504*p.^18 + 18360*p.^16 - 9520*p.^14 + 1820*p.^12).* cos(12*A);
% % FFit(:,223)=42^(1/2)*(4845*p.^20 - 15504*p.^18 + 18360*p.^16 - 9520*p.^14 + 1820*p.^12).* sin(12*A);
% % FFit(:,224)=42^(1/2)*(1140*p.^20 - 2907*p.^18 + 2448*p.^16 - 680*p.^14).* cos(14*A);
% % FFit(:,225)=42^(1/2)*(1140*p.^20 - 2907*p.^18 + 2448*p.^16 - 680*p.^14).* sin(14*A);
% % FFit(:,226)=42^(1/2)*(190*p.^20 - 342*p.^18 + 153*p.^16).* cos(16*A);
% % FFit(:,227)=42^(1/2)*(190*p.^20 - 342*p.^18 + 153*p.^16).* sin(16*A);
% % FFit(:,228)=42^(1/2)*(20*p.^20 - 19*p.^18).* cos(18*A);
% % FFit(:,229)=42^(1/2)*(20*p.^20 - 19*p.^18).* sin(18*A);
% % FFit(:,230)=42^(1/2)*(p.^20).* cos(20*A);
% % FFit(:,231)=42^(1/2)*(p.^20).* sin(20*A);
% % 

output(:,1:term)=FFit(:,1:term);
clear FFit;
FFit=output;

return