clc; clear all;

x=-2;y=6;z=3;
P = [x y z] %cartesian

rho = sqrt(x.^2 + y.^2);
phi = atan2d(y,x);
z = z;

P1 = [rho phi z] %cylindrical point

x1 = rho.*cosd(phi); %-1.95
y1 = rho.*sind(phi); %5.99
z1 = z;

rho1 = (y1.*cosd(phi)) + ((x1+z1).*sind(phi));
phi1 = (-(y1).*sind(phi)) + ((x1+z1).*cosd(phi));
z1 = 0;

A_p = [rho1 phi1 z1] %vector