%% Dirac with v = ao*t

function [current] = CurrentDiracLin(q,d,f,ao)

current = q*(exp(-2*1i*pi*f*sqrt(2*d/ao)).*((2*1i*pi*f*sqrt(2*d*ao)+ao)./((2*pi*f).^2))-ao./(2*pi*f).^2);