Contents

Abb1

clear
im = 0:150;
jm = 0:40;
X = zeros(length(im),length(jm));
Y = X;
Z = X;
%
for il=im
    for jl=jm
        X(il+1,jl+1) = (0.8 - 0.6*sin(sqrt(jl/80*pi))) .* cos(il/50*pi)  +0.08*sin(jl/40*pi);
        Y(il+1,jl+1) = (0.8 - 0.6*sin(sqrt(jl/80*pi))) .* sin(il/50*pi)  +0.08*sin(jl/40*pi);
        Z(il+1,jl+1) = 1.1*jl/40 + 0.7*(1-sin((150-il)/300*pi)).^2;
    end
end
%
Z = -Z;
figure,surf(X,Y,Z), shg
% >> axis off
% >> [a,e] = view
%
% a =
%
%   -21.1000
%
%
% e =
%
%    38.8000

Abb2

clear
im = 0:240;
jm = 0:60;
X = zeros(length(im),length(jm));
Y = X;
Z = X;
A = X;
B = X;
for il=im
    for jl=im
        A(il+1,jl+1) = 0.45 + 0.01*cos(il/3*pi) + jl/300*cos(il/240*pi).*cos((il/120*pi).^20);
        B(il+1,jl+1) = jl/300*cos(il/240*pi).*sin(il/120*pi)+0.125*(jl/60).^6 .* sin(il/120*pi);
        X(il+1,jl+1) = (0.4*cos(il/120*pi)+A(il+1,jl+1)) .* cos(jl/60*pi) + 0.48.*(jl/60).^6 .* sin(((il+60)/120*pi).^3);
        Y(il+1,jl+1) = 0.5 * sin(il/120*pi) + 0.01*sin(il/3*pi) + B(il+1,jl+1);
        Z(il+1,jl+1) = (0.45 + 0.4*cos(il/120*pi)) .* sin(jl/60*pi);
    end
end
%
figure, surf(X,Y,Z), shg

Abb3

clear
im = 0:800;
jm = 0:5;
X = zeros(length(im),length(jm));
Y = X;
Z = X;
A = X;
B = X;
C = X;
for il=im
    for jl=im
        X(il+1,jl+1) = cos(il/80*pi) + 3*jl/50*sin(21*il/800*pi);
        Y(il+1,jl+1) = sin((il/3200*pi).^0.1) * sin(il/80*pi);
        Z(il+1,jl+1) = 3*jl/50 + 0.5*sin(il/200*pi);
    end
end
figure, surf(X,Y,Z), shg

Abb4

clear
im = 0:150;
jm = 0:40;
X = zeros(length(im),length(jm));
Y = X;
Z = X;
%
for il=im
    for jl=jm
        X(il+1,jl+1) = (0.5 - 0.6*sin(jl/80*pi)) .* cos(il/50*pi)  +0.08*cos(jl/40*pi);
        Y(il+1,jl+1) = (0.8 - 0.6*sin(jl/80*pi)) .* sin(il/50*pi)  +0.08*cos(jl/40*pi);
        Z(il+1,jl+1) = 1.9*jl/40 + 2.*(1-sin((150-il)/300*pi));
    end
end
%
Z = -Z;
figure,surf(X,Y,Z), shg
% >> axis off
% >> [a,e] = view
%
% a =
%
%   -21.1000
%
%
% e =
%
%    38.8000
%
view(-21.2,38.8)
colormap jet
shading interp
axis off

Abb4spiel

clear
im = 0:150;
jm = 0:40;
X = zeros(length(im),length(jm));
Y = X;
Z = X;
%
for il=im
    for jl=jm
        X(il+1,jl+1) = (0.5 - 0.6*sin(jl/80*pi)) .* cos(il/50*pi)  +0.08*cos(jl/160*pi);
        Y(il+1,jl+1) = (0.8 - 0.6*sin(jl/80*pi)) .* sin(il/50*pi)  +0.08*cos(jl/160*pi);
        Z(il+1,jl+1) = 1.9*jl/40 + 2.*(1-sin((150-il)/300*pi));
    end
end
%
Z = -Z;
Y = Y.*X;
figure,surf(X,Y,Z), shg
% >> axis off
% >> [a,e] = view
%
% a =
%
%   -21.1000
%
%
% e =
%
%    38.8000
%
view(-21.2,38.8)
colormap jet
shading interp
axis off