Daniel@0: % PLOT_ELLIPSE
Daniel@0: % h=plot_ellipse(x,y,theta,a,b)
Daniel@0: %
Daniel@0: % This routine plots an ellipse with centre (x,y), axis lengths a,b
Daniel@0: % with major axis at an angle of theta radians from the horizontal.
Daniel@0: 
Daniel@0: %
Daniel@0: % Author: P. Fieguth
Daniel@0: %         Jan. 98
Daniel@0: %
Daniel@0: %http://ocho.uwaterloo.ca/~pfieguth/Teaching/372/plot_ellipse.m
Daniel@0: 
Daniel@0: function h=plot_ellipse(x,y,theta,a,b)
Daniel@0: 
Daniel@0: np = 100;
Daniel@0: ang = [0:np]*2*pi/np;
Daniel@0: R = [cos(theta) -sin(theta); sin(theta) cos(theta)];
Daniel@0: pts = [x;y]*ones(size(ang)) + R*[cos(ang)*a; sin(ang)*b];
Daniel@0: h=plot( pts(1,:), pts(2,:) );