I'd like knowing any possible ways to speedup the execution of the following function:

function beta = beta_c(k,c,gamma) beta = zeros(size(k)); E = @(x) (1.453*x.^4)./((1 + x.^2).^(17/6)); for ii = 1:size(k,1) for jj = 1:size(k,2) E_int = integral(E,k(ii,jj),10000); beta(ii,jj) = c*gamma/(k(ii,jj)*sqrt(E_int)); end end end