e10p7.m

% Example 10.7
alphae=[1 10];
A=[0 1;0 0]; B=[0;1]; C=[1 0]; D=0;
a11=A(1,1); Ae1=A(2,1);
A1e=A(1,2); Aee=A(2,2);
b1=B(1); Be=B(2);
G=polyvalm(alphae,Aee)*inv(A1e)*1