t=linspace(0,3e-4,200);
K=3.951e9
a=4.444e4
b=4.444e4
h=K/b*exp(-a*t).*sin(b*t);
plot(t,h)
grid
xlabel t
ylabel h(t)
title 'Impulse response'
