newton basis stuff seems to work, but I still use the K matrix to compute the needed values. This might be bad for ill-conditioned K
made requested corrections to NewtonBasisThoughts, still not sure why newton basis solns don't converge on right bdry