Define mathematica functions by simply sending the definition to the interpreter.
sage: m = mathematica sage: _ = mathematica('f[p_] = p^2'); sage: m('f[9]') 81