sanweitu.m

来自「最近一段数学建模源程序汇编,适用于建模初学者练习」· M 代码 · 共 16 行

clear,clf,
a=-1;b=1;c=-15;d=15;n=20;eps1=0.01;
x=linspace(a,b,n);y=linspace(c,d,n);
[X,Y]=meshgrid(x,y);
for i=1:n
    for j=1:n
        if(1-X(i,j))<eps1|X(i,j)-Y(i,j)<eps1
            z(i,j)=NaN;
        else
            Z(i,j)=1000*sqrt(1-X(i,j))^-1.*log(X(i,j)-Y(i,j));
        end
    end
end
zz=-20*ones(1,n);plot3(x,x,zz),grid off,hold on
surf(X,Y,z)
xlabel('x'),ylabel('y'),zlabel('z'),box on