%% Dirac with vo

function [current] = CurrentDiracCst(Q,d,f,vo)

current = Q*sinc(pi*f*d/vo).*exp(-1i*2*pi*f*d/vo);