function V =velvar(t,f)
if t<=50
    V=0;
elseif t>50 & t<=100
    V=-(0.012*pi*f*sin((2*pi*f.*t)));
else
    V=(0.04*sin((2*pi*f.*t)));
end