clear;
m=0.166; //hmotnost telesa [kg]
LAB=0.5; //vzdalenost |AB| [m]
tAB=0.382; //cas za kterou teleso na naklonene rovine urazilo drahu LAB
g=9.81; //tihove zrychleni [m.s-2]
a=20/180*%pi; //uhel naklonene roviny (20°) [rad]
f=0.18; //soucinitel smykoveho treni [-]
//odpor vzduchu - vychazi z kvadratickeho odporu (1/2*rho*S*v^2), kde pro rychlosti okolo 1m/s lze uvazovat linearni zavislost na rychlosti.
c=1/2*1.27*0.03*0.04*1.05;

// substituce dle ulohy ze cviceni (gg=g*(sin(a)+f*cos(a)), kde jsme uvazovali pohyb "do kopce") pro pohyb smerem "z kopce" uvazujeme uhel a<0, takze sin(-a)=-sin(a)
gg=g*(-sin(a)+f*cos(a));

//vztahy pro rychlosti vA a vB odvozene na cviceni
vA=gg*((tAB+c/(m*gg)*LAB)/(1-exp(-c/m*tAB))-m/c);
vB=vA-gg*tAB-c/m*LAB;

disp('vA= '+ string(vA)+ ' m.s-2')
disp('vB= '+ string(vB)+ ' m.s-2')

fps=960;
//cas v ms na 1cm
disp('cas na 1cm drahy v bode A: '+string(0.01/vA*fps)+' ms')
disp('cas na 1cm drahy v bode B: '+string(0.01/vB*fps)+' ms')