function PredictionMatrix input Real[:, :, :] S "step response matrix"; input Integer Nl "lower prediction horizon"; input Integer Np "upper prediction horizon"; input Integer Nu "control horizon"; input Integer Nm "model horizon"; input Integer m "number of control systems inputs"; input Integer p "number of control systems outputs"; output Real P[(Np - Nl + 1)*p, Nu*m] "prediction matrix"; end PredictionMatrix;