代码如下:
a=10;
clf;
x=-a:0.2:a;y=-a:0.2:a;
for i=1:length(y)
for j=1:length(x)
if 0.85 u(i,j)=10-10*x(j)/y(i); else 1<=x(j)/y(i)<1.15 z(i,j)=10*x(j)/y(i)-10; end end end surf(x,y,u);
代码如下:
a=10;
clf;
x=-a:0.2:a;y=-a:0.2:a;
for i=1:length(y)
for j=1:length(x)
if 0.85 u(i,j)=10-10*x(j)/y(i); else 1<=x(j)/y(i)<1.15 z(i,j)=10*x(j)/y(i)-10; end end end surf(x,y,u);