1D Lorentz

(%i1) f1(x,x0,g) := ( g^2 / ( ( x - x0)^2 + g^2) );

Result

2D symmetric Lorentz

(%i2) f2(x,y,x0,y0,g) := ( g^2 / ( ( x- x0)^2 + (y-y0)^2 + g^2));

Result

ND Lorentz

(%i3) fn(x,u,S) := 1/( 1 + transpose(x-u) . invert(S) . (x-u));

Result

Specialize 2D

(%i4) fn(matrix([x,y]),matrix([x0,y0]),matrix([s00,s01],[s10,s11]));

Result

(%i5) fortran(%);

Result


Created with wxMaxima.