it mthrfn works! need to test convergence