in spite of the title this page contains nothing about parcellfun,
I have built the following test routine:
vector_in = 1:5;
vector_out = pararrayfun(5,parquad(procid), vector_x)
function out=parquad(procid)
switch procid
case{1}
out=quadgk(@(r) ntgrnd1(r),0,1);
case{2}
out=2*quadgk(@(r) ntgrnd1(r),0,1);
case{3}
out=3*quadgk(@(r) ntgrnd1(r),0,1);
case{4}
out=4*quadgk(@(r) ntgrnd1(r),0,1);
case{5}
out=5*quadgk(@(r) ntgrnd1(r),0,1);
endswitch
endfunction