clear all home disp('create a matrix: ') m = magic(6) disp('isolate the bottom left corner with size 2*2 of the previous matrix: ') M = bottom_left(magic(6) , 2) disp('call hulk function with a vector of your choice: ') v = [1 2 1] h = hulk(v) disp('compute and place the roots of each formatted polynomial in matrix h into a new matrix H: ') H = zeros( (size(h,1) - 1) , size(h,2) ); for i = 1:size(h,2) H(:,i) = roots(h(:,i)); end H