1 2 3 4 5 6
f: exp(-r^2/sigma^2)/sigma; laplacef: ratsimp(diff(r * diff(f,r), r) / r); f: exp(-r^2/(2*sigma^2))/(sqrt(2)*sigma); laplacef: ratsimp(diff(r * diff(f,r), r) / r);