is there an any alternate solution for this "P=freqcp(​numpl,denp​l,w);" command?

clc; clear all; close all; tic k = 1:0.5:10; a = 1:0.5:5; n1 = numel(k); n2 = numel(a); numpl = repmat(k, 1, n2)'; denpl = repmat([ones(n2,1), a.', zeros(n2,1)] ,n1, 1); w=[.1,.5,.8,1,2,8,15,50,100]; P=freqcp(numpl,denpl,w);

% plot templates plottmpl(w,w,P);

time=toc