Cointegration function for many variables

Is there a function or script similar to corrcoef that can return cointegration results in a matrix form, without the need for additional toolboxes? Input will have thousands of variables and observations.

Thanks