Codigo en m'-file para conversion entre coordenasa

opc=0;
while opc~=7;

disp('----------------------------------------------------------');
disp('Conversiones entre coordenadas ');disp('----------------------------------------------------------');



disp('Opcion 1: Cartesiana a polar ');
disp('Opcion 2: Cartesiana a cilindrica ');
disp('Opcion 3: Cartesiana a esferica ');
disp('Opcion 4: Polar a cartesiana ');disp('Opcion 5: Cilindrica a cartesiana ');
disp('Opcion 6: Esferica a cartesiana ');
disp('Opcion 7: Salir ');

opc=input('\nIngrese opcion: ');



switch opc

case 1fprintf('\nIngrese valores \n');
x=input('\nIngrese x: ');
y=input('\nIngrese y: ');

r=sqrt(x^2+y^2);
theta= atand(y/x);


fprintf('El valor delradio r es = %f\n',r);
fprintf('El valor del angulo thetha(grados) es = %f\n',theta);

case 2
fprintf('\nIngrese valores \n');
x=input('\nIngrese x: ');y=input('\nIngrese y: ');
z=input('\nIngrese z: ');

r=sqrt(x^2+y^2);
theta= atand(y/x);
z=z;

fprintf('El valor del radio r es = %f\n',r);fprintf('El valor del angulo thetha(grados) es = %f\n',theta);
fprintf('El valor de la altura z es = %f\n',z);

case 3
fprintf('\nIngrese valores \n');
x=input('\nIngresex: ');
y=input('\nIngrese y: ');
z=input('\nIngrese z: ');

p=sqrt(x^2+y^2+z^2);
theta= atand(y/x);
phi=acosd(z/(sqrt(x^2+y^2+z^2)));fprintf('El valor del radio r es = %f\n',r);
fprintf('El valor del angulo thetha(grados) es = %f\n',theta);
fprintf('El valor del angulo phi es = %f\n',phi);

case 4fprintf('\nIngrese valores \n');
r=input('\nIngrese radio r: ');
theta=input('\nIngrese angulo theta(grados): ');

x=r*cosd(theta);
y=r*sind(theta);...