Hi @Anal
The error in your MATLAB script likely arises from the way you are using the “whittakerW” function and the “quad2d” function.
Ensure that the function handles you are using are properly vectorized. The “quad2d” function expects the integrand to accept arrays as inputs and return an array of outputs.
Make sure that the “whittakerW” function is available in your MATLAB environment. This function is part of the Symbolic Math Toolbox, so ensure you have it installed and the toolbox is loaded.
Kindly refer to the following code snippet for understanding.
% Define the function Part_1
Part_1 = @(r, theta) 1 + (2 .* (1 + r .* cos(theta))).^2;
% Define the function W_i, ensuring vectorization
W_i = @(r, theta) Part_1(r, theta) .* whittakerW(50, 0.5, (2 .* r) ./ 50);
% Perform the double integration using quad2d
Matrix_element = quad2d(W_i, 0.5, 2000, 0, pi);
% Display the result
disp(Matrix_element);
Kindly refer following MathWorks documentation for more information:
I hope this will be helpful.
