function V =velvar(t,f)
if t<=50
    V=0;
elseif t>50 & t<=80
    %V=12*(0.01*2*pi*f*cos((2*pi*f.*t)));
    V=-0.125;
else
    V=25*(0.02*2*pi*f*cos((2*pi*f.*t)));
end

