(%i1)
f(x,u,S,k) := gamma((1+k)/2) /
( gamma(1/2) * %pi^(k/2) * determinant(S)^(1/2) *
( 1 + transpose( x - u) . invert(S) . ( x - u ))^((1+k)/2));
(%i2)
q:f(matrix([x,y]),matrix([x0,y0]),matrix([s00,s01],[s10,s11]),2);
(%i3)
fortran(q);
Created with
wxMaxima
.