m=zeros(6,15); %This creates a matrix of 0's where rach polynomial is entered. If order is changed, so should M CT=zeros(100,7); %This is another matrix of 0's so ...

The function polyval should be generalized and vectorized in order to receive an array of values at which the polynomial is to be evaluated. The function polyval_der could also be created; this would ...