function [rho,theta]=tran(x,y)
rho = sqrt(x^2 + y^2);
theta=atan(y/x);