function [x,y] = circle(R,x0,y0)

%	<<<  COMPUTING FREQUENCY DOMAIN CIRCLES >>>
%
%	USAGE:	[RE,IM] = CIRCLE(R,X0,Y0)
%
%	The circles are computed from the equation:
%
%		z'*z = R^2
%
%	Input to the function is circle radius R and the center of the 
%	circle (x0,y0). Output is matched real and imaginary parts
%	of the circle at 50 linearily scaled points.

% S. Toffner-Clausen , Last revised:  95.08.17
% Copyright (c) by the authors
% All Rights Reserved.

t = 0:0.02:1;
z = [cos(2*pi*t)              % circle z'*z = 1
     sin(2*pi*t)];
x1 = R*z; 
x = x0 + x1(1,:);
y = y0 + x1(2,:);