D.3.2.17 hessenberg

Procedure from library linalg.lib (see linalg_lib).

Usage:

hessenberg(M); matrix M

Assume:

M constant square matrix

Return:

matrix H; Hessenberg form of M

Example:

LIB "linalg.lib";
ring R=0,x,dp;
matrix M[3][3]=3,2,1,0,2,1,0,0,3;
print(M);
print(hessenberg(M));