| (%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); |