function ell=ellrun(ucoef,vcoef,n);
%ell=ellrun(ucoef,vcoef,n);
% ucoef and vcoef are generated from tidef.m
% n refers to the index of the frequency (based on input into tidef)
% default: M2: 1, M4: 2, M6: 3
% ell is a matrix of ellipse parameters.
% see help for ellipse.m for details.
n=n*2;
ell=ellipse(ucoef(:,n)',ucoef(:,n+1)',vcoef(:,n)',vcoef(:,n+1)');
ell=ellfix(ell);
