MATLAB Poynting Vector integral without symbolic pkg

IntegralPoyntingVector.m

% Define farfield radiation electric field (E) and magnetic field (H) functions
E = @(r, theta, phi) your_electric_field_function(r, theta, phi);
H = @(r, theta, phi) your_magnetic_field_function(r, theta, phi);

% Define spherical region of integration
rmin = your_min_r_value;
rmax = your_max_r_value;
thetamin = your_min_theta_value;
thetamax = your_max_theta_value;
phimin = your_min_phi_value;
phimax = your_max_phi_value;

% Create anonymous function for Poynting vector in spherical coordinates
PoyntingVectorSpherical = @(r, theta, phi) r^2 * sin(theta) * real(cross(E(r, theta, phi), conj(H(r, theta, phi))));

% Perform the integral over the specified spherical region
integral_result_spherical = integral3(PoyntingVectorSpherical, rmin, rmax, thetamin, thetamax, phimin, phimax);

disp(['Integral of Poynting Vector over the specified spherical region: ', num2str(integral_result_spherical)]);

• You do not need to do a Poynting Vector integral over r to getting radiation power, Only over Theta and Phi is needed!