Как численно решить сложный интеграл, используя MATLAB

Question

Я пытаюсь решить это уравнение:

K = SQRT ((R * T)/(4 * пи * лямбда)) * интеграл от -inf до inf (- ((lambda + F * neta)/R * Tx)^2 * R * T/4 * lambda)/exp (x) +1 по отношению к x

где, neta является интервал от 0 до 1, а остальные символы (R, T, F, лямбда и pi) имеют постоянные значения.

Я пытался использовать эти коды:

код 1

clear all; 
close all; 
clc; 
F = 96485.34; 
R = 8.3145; 
T = 298.15; 
lambda = 0.2; 
neta=0:0.1:1; 
pi=3.1415; 
f=@(x) exp(-((lambda+F*neta)/R*T-x).^2*R*T/4*lambda)/(exp(x)+1); 
Q=integral(f,-inf,inf); 
k= sqrt((R*T)/(4*pi*lambda)).*Q 

код 2

clear all; 
close all; 
clc; 
F = 96485.34; 
R = 8.3145; 
T = 298.15; 
lambda = 0.2; 
neta=0:0.1:1; 
pi=3.1415; 
x= 0:100; 
f(x)=exp(-((lambda+F*neta)/R*T-x).^2*R*T/4*lambda)/(exp(x)+1); 
q=quadl('f', 0, 100); 
k= sqrt((R*T)/(4*pi*lambda)).*q 

но эти коды возврата ошибок, которые я не знаю, чтобы решить. Может кто-нибудь мне помочь, пожалуйста?

спасибо

Answer 1

Проблема с neta. Если это скалярное значение, то код в порядке. Но как вы предлагаете интегрироваться по интервалу neta по отношению к другой переменной? Что это значит? Я попытался взять отдельные значения для neta от 0 до 1 и вычислить значение Q для каждого, но каждый Q возвращает 0.