clear,clc;
x=0:pi/10:2*pi;
y=x;
[X,Y]=meshgrid(x,y);%转换为空间点坐标
Z=sin(X).sin(2X+3);
surf(X,Y,Z)
box on
xlabel(‘x轴’),ylabel(‘y轴’),zlabel(‘z轴’)
title(‘三维曲面图’)
grid on

绘制同心圆
fcs.m
a=input(‘输入横坐标:’);
b=input(‘输入纵坐标:’);
r=input(‘输入半径:’);
step=3;
fcircle(a,b,r);

function fcircle=(a,b,r)
for rx=r:3:r+9
theta=0:pi/100:2pi;
x=a+rxcos(theta);
y=b+rx*sin(theta);
plot(x,y,‘b-’)
hold on
axis equal
end
end