pfaffian(M:MAT) where M is skew-symmetric. The resulting type depends on the entries of the matrix.
/**/ Use R ::= QQ[x,y]; /**/ pfaffian(mat([[0,y],[-y,0]])); y