Rotatable Ellipsoid

// This creates an ellipsoid
// INPUTS:
//   wname - name of ywave having 361 points
//              (presumes xwave of same length is named wname + "x")
//   (xo,yo) - center position of ellipse
//   (a, b) - elliptical dimesions to above
//   alpha - angle of rotation (counterclockwise)
// OUTPUTS:
//   make all changes directly in ywave and xwave
 
Function Ellipsoid(wname,xo,yo,a,b,alpha)
    string wname
    variable xo,yo,a,b,alpha
 
    wave yw = $wname
    wave xw = $(wname +"x")
    
    make/FREE/n=361 urx,ury
    
    variable cosalpha = cos(alpha*Pi/180), sinalpha = sin(alpha*Pi/180)
    
    urx = a*(cos(p*Pi/180))
    ury = b*(sin(p*Pi/180))
    
    xw = xo + cosalpha* urx - sinalpha*ury
    yw = yo + sinalpha*urx + cosalpha*ury
        
    return (0)  
end