fixed BC problem in Newton basis, but calculation of basis is definitely unstable. Need to figure out how to compute stably or find basis elem.s in closed form