update issue #36 
checking the matrices with respect to various trial functions