
% Generating the pattern of the tail of a mamal

k1=[1;0];
k2=[-0.5;1/sqrt(2)];
k3=[0.5;1/sqrt(2)];
phi1=0;
phi2=0;
phi3=-phi1-phi2;
phi3=phi3+0*pi;

%phi1=pi/3;
%phi2=pi/3;
%phi3=pi/3;

r=1;
s=10;
Npts=100;
z=linspace(0,s,Npts);
l=linspace(0,2*pi*r,Npts);
[Z,L]=meshgrid(z,l);

U=cos(n*L/r).*cos(Z*m*pi/s);

fs=16;

figure(1)
clf
set(gca,'FontSize',[fs]);
surfl(X,Y,U);
shading interp
colormap(gray)
axis tight

figure(2)
clf
%subplot(4,1,1)
set(gca,'FontSize',[fs]);
surf(X,Y,U);
shading interp
colormap('default')
colormap(jet)
colormap(gray.^(0.5))
%colormap([0 0 0 ; 1 1 1]);
colormap([gray;gray*0+1;gray*0+1])
colorbar
view([0 90])
axis tight

