Интеграция IDL

Question

Я хочу интегрировать функцию, которую я создаю, но функция изменит каждую итерацию на основе заданного ввода. Например:

y=4e^(mx/4) 

Я хотел бы интегрировать по й с нижней и верхней границей, но значение т изменится. Я знаю все мои значения m.

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

Answer 1

QROMB (и другие интеграторы) хотят функцию от одной переменной, поэтому вам нужно получить m там через заднюю дверь. Один из способов заключается в общем блоке:

function integrand,x 
common int_common,int_m 
return,4*exp(int_m*x/4) 
end 

function integrator,m,xlow,xhigh 
common int_common,int_m 
int_m=m 
return,qromb('integrand',xlow,xhigh) 
end 

интегратор (m, xlow, xhigh) вернет необходимый вам интеграл.