function psi = mexican_hat(t,sigma)
    psi = 2/(sqrt(3*sigma)*pi^0.25) *(1-t.^2/(sigma^2)).*exp(-t.^2/(2*sigma^2));
    
    