function bod_na_spirale(omega,v0,R,t)
%bod_na_spirale - kreasli drahu bodu pohybujiciho se po spirale
%omega..uhlova rychlost[rad/s]
%v0..pocatecni rychlodst bodu [m/s]
%R..polomer spiraly [m]
%t....casovy okamzik pro zvyrazneni [s], impl. t=30s
%priklad volani: bod_na_spirale(.3,.1,4,21),bod_na_spirale(.3,.5,4,5)


if nargin<4 || t<0
    t=30;
end
if R<=0
   error('Polomer musi byt kladny!') 
end
if v0<=0
   error('Pocatecni rychlost musi byt kladna!') 
end
if omega<0
    error('Uhlova rychlost musi kladna!')
end

R=real(R(1));    %kdyby R bylo komplexni a nebo vektor
v0=real(v0(1));   
omega=real(omega(1));   

%kreslime

tp=1:.5:max(t+10,150);  %vybere se vetsi z hodnot t+10 a 150, osetreni aby bod v case 30 s byl videt na spirale
x=R*sin(omega*tp);
y=R*cos(omega*tp);
z=v0*tp;
plot3(x,y,z)
hold on
plot3(R*sin(omega*t),R*cos(omega*t),v0*t, ...
    'color','r','marker','o','markersize',10);    %zvyrazneny bod v case t
title('draha bodu po spirale')
xlabel('x'),ylabel('y'),zlabel('z')
axis equal