untitled2.asv

最近一段数学建模源程序汇编,适用于建模初学者练习

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=-20on