function x = frot2(w,th)
% Function to compute normalised radius of a rotating star.
% Equations from Collins, 1963.
% Note that theta ( = th here) must have dimensions of degrees.
% to use: x = frot2(w,th)
%
% by D. Holmgren, Jan 18 94. holmgren@phobos.astro.uwo.ca
b = zeros(3,11);
b(1,1)=1.0001956;
b(3,2)=-1.183952e-09;
b(2,3)=4.2098469e-05;
b(1,10)=-8.0843068e-02;
b(2,10)=2.5709925e-04;
b(3,10)=-5.598946e-08;
b(1,11)=8.0509197e-02;
b(2,11)=-2.545361e-04;
b(3,11)=5.9308193e-08;
y = zeros(11,1);
for j = 1:11
 y(j) = w^(j-1);
end
c = b*y;
v = [1;th^2;th^4];
x = c'*v;