Rutina con matlab para graficar

clear all
clc
disp('**********Opciones programa proyecto final**********')
disp('Opcion 1 Grafica 3D')
disp('Opcion 2 Grafica de Nivel')
disp('Opcion 3 Derivadas Cruzadas')disp('Opcion 4 Grafica de la Gradiente')
disp('Opcion 5 Funcion de la Divergencia')
disp('Opcion 6 Rotacional de la funcion')
disp('Opcion 7 Continuar o Salir')
K=input('Digite 1 parainiciar ');
while K==1
A=input('Seleccione la Opcion Deseada ');

switch A

case 1
[x,y]=meshgrid(-150:0.2:150,-150:0.2:150);
z=x.^4+(y-50).^4;mesh(x,y,z)
xlabel ('Eje X')
ylabel ('Eje Y')
title('Grafica 3D')
pause
close all
case 2
[x,y]=meshgrid(-150:0.2:150,-150:0.2:150);z=x.^4+(y-50).^4;
contour3 (x,y,z,16);
xlabel ('Eje X')
ylabel ('Eje Y')
title('Grafica de nivel')
pause
close all
case3
syms x y
z=x.^4+(y-50).^4;
diff(z,x)
diff(z,y)

case 4
[x,y]=meshgrid(-100:10:100);
z=x.^4+(y-50).^4;[px,py]=gradient(z,0.1,0.1);
contour(z);
quiver(px,py)
xlabel ('Eje X')
ylabel ('Eje Y')
title('Grafica de la Gradiente')
pauseclose all
case 5
syms x y z
m=(x/(x^2+y^2)^(3/2));
n=(y/(x^2+y^2)^(3/2));
r=1;
f=diff(m,x)+diff(n,y)+diff(r,z);
pretty(f)case 6
syms x y z
m=(x/(x^2+y^2)^(3/2));
n=(y/(x^2+y^2)^(3/2));
r=1;
F=[diff(r,y)-diff(n,z) diff(m,z)-diff(r,x) diff(n,x)-diff(m,y)];pretty(F)
case 7
disp('Digite 1 Para volver al menu u otra tecla para salir ');
K=input('Digite su opcion');
end
end
disp('Gracias por visitarnos')