f=vikingplot()

function f=vikingplot()

global spunkter rpos botten nx ny
ztop=0.5;             % 		ztop initieras.

ny=10;		      % Antal indelningar p bredden.
nx=45;       	      % Antal indelningar p lngden.
n1=round(nx*0.45);    % Del av nx fr den frsta bzierkurvan.
n2=nx-n1-1;           % Del av nx fr den andra bzierkurvan.

% Start- och stopppunkter fr relingens tv bzierkurvor.
p1=[0 0 ztop+0.5]; p2=[2 0.8 ztop]; p3=[4 0 ztop+0.7];

% Styrpunktsavstnd och riktningar fr relingens tv bzierkurvor.
alfa1=0.7; alfa2=1.2; gamma1=1.5; gamma2=0.7;
k1=[3 2 -2]; k2=[1 0 0]; k3=[-3 2 -1];

% Styrpunkter fr relingens tv bzierkurvor.
b1=p1+alfa1*k1/norm(k1);
c1=p2+gamma1*(-k2);
b2=p2+alfa2*k2;
c2=p3+gamma2*k3/norm(k3);

% Uppbyggnad av bzierkurvorna p bda sidor om y-axeln.
rpos=[bezier(p1,b1,c1,p2,n1); 
      bezier(p2,b2,c2,p3,n2);
      p3];
rneg=[rpos(:,1) -rpos(:,2) rpos(:,3)];

% Klens koordinater.
botten=[rpos(:,1) 0*rpos(:,2) vbottfunc(rpos(:,1), ztop+0.5,ztop+.7)];

% Uppbyggnad av skrovets x-, y- och z-punkter till SURFL-kommandot.
% Skrovet plottas utmed x-axeln.
X=[]; Y=[]; Z=[];
for i=1:nx
  Q=rpos(i,:); P=botten(i,:); Q2=Q.*[1 -1 1]; % Start- och stoppunkter.

  alfa=0.4*norm(Q-P); gamma=0.6*norm(Q-P);    % Styrpunktsavstnd
  k1=[0 0 -1]; k2=[0 1 1];                    % Styrpunktsriktingar
  
  b1=Q+alfa*k1;                               % Styrpunkter
  c1=P+gamma*k2/norm(k2);
  b2=(c1-P).*[1 -1 1]+P;
  c2=b1-Q+Q2;
  
  r=[bezier(Q,b1,c1,P,ny);                    % Bzierkurvor
     bezier(P,b2,c2,Q2,ny);
     Q2];   
  X=[X r(:,1)];                               % Punkter till SURFL
  Y=[Y r(:,2)];
  Z=[Z r(:,3)];
end

% Skrovet ritas upp med SURFL-kommandot.
surfl(X,Y,Z), hold on

% Uppbyggnad och utritning av skldar.
% Utnyttjar funktionen RITADELTA.
n=20; r=0.16;col=[0.3 0.2 0.1];
t=0:2*pi/n:2*pi;
x=r*cos(t);y=-0.9*ones(1,size(t,2));z=r*sin(t)+0.5;
for i=1.25:0.5:2.75
  fill3(x+i,y,z,col)
  fill3(x+i,-y,z,col)
  ritadelta(i+.1,-1,.54,0.015)
end

% Uppbyggnad och utritning av mast.
n=10; r=0.03;
t=0:2*pi/n:2*pi;
xx=r*cos(t)+2;yy=r*sin(t);zz=0*ones(size(xx));
xm=[];ym=[];zm=[];
for i=0:0.1:2.7
  xm=[xm;xx];ym=[ym;yy];zm=[zm;zz+i];
end
surfl(xm,ym,zm)

% Uppbyggnad av seglet.
sailsurface=0.2;
sailtop=2.5; sailbottom=0.7; sailhalfwidth=1.05; sailstep=0.1;
xm=[]; ym=[]; zm=[];
for z=sailbottom:sailstep:sailtop
  xx=[]; yy=[]; zz=[];
  for y=-sailhalfwidth:sailstep:sailhalfwidth
    yy=[yy y];
    zz=[zz z];
    xx=[xx 2.3-(z-(sailtop-sailbottom)/2-0.7)^2/3-y^2/5+sailsurface];
  end
  xm=[xm;xx];ym=[ym;yy];zm=[zm;zz];
end

% Uppbyggnad av magisk frgskalningsmatris som ser till att
% endast en begrnsad del av vald 'CLUT' anvnds.
magiccolornumber=(46+30)/(30+46+30);
sailmin=min(min(zm)); sailmax=max(max(zm)); sailrange=sailmax-sailmin;
magicallyscaledmatrix=(1-magiccolornumber)/sailrange*zm;
scaledsailmin=min(min(magicallyscaledmatrix));
finalmagicmatrix=magicallyscaledmatrix-(scaledsailmin-magiccolornumber);

noofstripes=7;
stripewidth=2*sailhalfwidth/noofstripes/sailstep;
for i=1:2*stripewidth:noofstripes*stripewidth
  for j=i:1:i+stripewidth-1
    finalmagicmatrix(:,round(j))=finalmagicmatrix(:,round(j))-magiccolornumber;
  end
end
  
% Utritning av seglet med magiska frger.
sail=surf(xm,ym,zm,finalmagicmatrix);

% Utritningsinstllningar.
title('Roddbt DeLuxe av d94-fsa och d94-pax')
colormap(vcol)
axis([0 4 -2 2 0 4]), axis square, axis off, view ([47 15])