def funfit( (xvec, yvec), basisfuns ):
    return Mat([ map(form,xvec) for form in basisfuns ]).tr().solve(Vec(yvec))
Solves design matrix for approximating to basis functions
dataset/ETHPy150Open sunlightlabs/clearspending/completeness/statlib/matfunc.py/funfit