theta=(pi/180)*[0 10 20 30 40 50 60 70 80 90];
ft=[0 17.3 34.2 50 64.3 76.6 86.6 93.9 98.4 100];
for m=1:length(theta)
    for n=1:length(theta)
        X(m,n)=theta(m).^(10-n);
    end 
end
y=ft';
A=inv(X)*y;
x=0:0.01:pi/2;
p=A(1).*x.^9 + A(2).*x.^8 + A(3).*x.^7 + A(4).*x.^6 + A(5).*x.^5 + ...
    A(6).*x.^4 + A(7).*x.^3 + A(8).*x.^2 + A(9).*x.^1 + A(10);
polar(x,p);