function A = ret_mat_A( h ) 

A = [ 1 h; -h 1 - h^2 ];