newton BDF function TP6
fixed diffdiv function
function to construct the matrix A
cholesky related function
function to compute z minimising ||Az - y||^2
function to compute z minimising ||Az - y||^2 with QR factorisation
plotting function