function y=f(c) n=numel(c);y=[[0*(3:n);eye(n-2)] -c(1:n-1)'/c(n)]; end % neede in online interpreter only f([23., 1., 92., 8., -45., 88., 88.]) % call function
[23., 1., 92., 8., -45., 88., 88.]
ans = 0.00000 0.00000 0.00000 0.00000 0.00000 -0.26136 1.00000 0.00000 0.00000 0.00000 0.00000 -0.01136 0.00000 1.00000 0.00000 0.00000 0.00000 -1.04545 0.00000 0.00000 1.00000 0.00000 0.00000 -0.09091 0.00000 0.00000 0.00000 1.00000 0.00000 0.51136 0.00000 0.00000 0.00000 0.00000 1.00000 -1.00000