1D Lorentz
(%i1) | f1(x,x0,g) := ( g^2 / ( ( x - x0)^2 + g^2) ); |
2D symmetric Lorentz
(%i2) | f2(x,y,x0,y0,g) := ( g^2 / ( ( x- x0)^2 + (y-y0)^2 + g^2)); |
ND Lorentz
(%i3) | fn(x,u,S) := 1/( 1 + transpose(x-u) . invert(S) . (x-u)); |
Specialize 2D
(%i4) | fn(matrix([x,y]),matrix([x0,y0]),matrix([s00,s01],[s10,s11])); |
(%i5) | fortran(%); |