Date: Jan 20, 2013 1:36 PM
Author: Francesco Perrone
Subject: Speedup this function

Hallo everybody,

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

I thank you all in advance.

Best regards,
fpe