function IntNegNuBes=FpexColPosNu(y,nu,Beta);  
% sumS fr Fpex Collection wth Uniform IC 1/(b-a) 
% definte intgral from a to b to X (x^nu)*BesselJ(+nu,Beta*x) 
% . for matix dimensions D. Bakewell 19:00 Mo/16/8/1999
format long;
ynu=y.^nu;
IntNegNuBes=ynu.*BesselJ(nu,y*Beta);