function [r, theta] = polarcoordinates(x, y)
% polarcoordinates(x, y) returns the polar coordinates 
% of the point with rectangular coordinates (x, y).
r = sqrt(x^2 + y^2);
theta = atan2(y, x);