tried using 2011 newton basis, there might be a numerical problem (in colloc_1dbvp: increase N, value matrix has NaN entries)
- condition # decreases quickly for odes?
- very slightly better in condition # for pde