У меня есть функция,
F = exp (-B (i, j). * Cos (x))
где B (i, j) - постоянная матрица размерности i-j
Теперь я хочу численно интегрировать вышеуказанную функцию w.r.t 'x' в интервал (0, pi / 2) в Matlab.
Заранее спасибо !!!
function Z = Z_normalization(E0,f)
%Evaluating the Normalization partition function for the calculation
% Detailed explanation goes here
% E0 = electric field amplitude;
% f = frequency
% B(i,j) = product of electric dipole moment and field
%f = 10:100:1000000;
%E0 = [650, 750, 850, 950, 1000];
for j = 1:length(E0)
for i = 1:length(f)
%B(i,j) = E0(i).*f(j);
B = Scalar_electric_potential(E0,f);
F = @(x) exp(-B(j,i).*cos(x)).*cos(2.*x);
Z(j,i) = integral(F, 0, pi/2)
end
end
end
Здесь B (i, j) является функцией от f, E0. Таким образом, интегрирование прекрасно, когда у меня нет exp (F). При добавлении exp (F) результат интегрирования одинаков для всех значений.