changed 1d code to no longer enforce derivative constraints on boundary points. Now we have rank deficiency problems (???)