function [ theta, r ] = mycart2pol(x,y,scelta)

% mycart2pol calcola le coordinate polari theta (tra 0 e pi greco) e r
% dalle coordinate cartesiane (x,y); se scelta = 'gradi' allora
% theta va calcolato in gradi, altrimenti in radianti.

[theta, r] = cart2pol (x,y);
if theta < 0
  theta = theta + 2* pi;
endif;
if strcmp (scelta, 'gradi')
  theta = theta * 180/pi;
endif;

endfunction