PCE Coefficients calculation at multiple frequency points

UQ with UQLab Community Q&A and How To
1

Hi, I have a governing equation:
f = 1:5000;
modelopts.mString = '(sqrt(1jP.wP.mju_0*(X+1jP.wP.epsju_rP.epsju0)))’;
modelopts.Parameters.w = 2pi*f;
using these above lines of code, its taking too much time, is there any fast way to calculate PCEs at all the points.? Thanks