changed calculate_newton_basis, now v*v'-KM ~ 0 but v != tril(v)
I guess this doesn't really matter and is an artifact of the adaptive point selection
Merge branch 'master' of https://github.com/palmerlao/iit-pde-project
merging tim's changes