Hi,
Is there a way to speed up this?
maxN = 120;x = -30:0.1:30;xElements = numel(x);u_mn = zeros(xElements, xElements);for m = -maxN:2:maxN for i = 1 : xElements for j = 1 : xElements u_mn(i, j) = sqrt((n+1)/pi) * besselj(m+1, 2*sqrt(x(i)^2 + x(j)^2)) / sqrt(x(i)^2 + x(j)^2)^(m+1) * (x(i) + 1i*x(j))^m; end endend
Best regards, Alex
Best Answer