D=2.9;
s=2*pi*rand;
M=10;
L=1500;
Li=240;
G=2;
r=1.5;
x=1:10:1000;
y=1:10:1000;
m=1;
n=0;
W=L*[(G/L)^(D-2)]*sqrt(log(r)/M);
ni=19;
Z=0;
[X,Y]=meshgrid(x,y);
while(m<11)
while(n<ni)
S=[r^((D-3)*n)]*[cos(s)-cos((2*pi*(r^n))*sqrt(X.^2+Y.^2)*cos(atan(Y./X)-pi*m/M)/L+cos(s))];
Z=W*S+Z;
n=n+1;
end
m=m+1;
end
surf(X,Y,Z) |