例
s = %s; p = [s ; s*(s+1)^2 ; 2*s^2+s^3] [Y, rk, ac] = colcompr(p*p'); Y clean(ac) |  |  |
--> p = [s ; s*(s+1)^2 ; 2*s^2+s^3]
p =
s
s +2s² +s³
2s² +s³
--> [Y, rk, ac] = colcompr(p*p');
--> Y
Y =
1 +2s +s² 2s +s² 1
-1 0 0
0 -1 0
--> clean(ac)
ans =
0 0 s²
0 0 s² +2s³ +s⁴
0 0 2s³ +s⁴