Intégration IDL

Question

Je cherche à intégrer une fonction que je construis, mais la fonction changerait chaque itération en fonction d'une entrée donnée. Par exemple:

y=4e^(mx/4) 

je voudrais intégrer par rapport à x avec une limite inférieure et supérieure, mais la valeur de m changerait. Je connais toutes mes valeurs de m.

Puis-je travailler avec ça? Mon hypothèse initiale serait d'utiliser QROMB mais cela semble limité et incapable de gérer mon problème.

Answer 1

QROMB (et d'autres intégrateurs) veulent une fonction d'une variable, vous devez donc faire passer le m par la porte arrière. Une façon est un bloc commun:

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 

intégrateur (m, Très bas, xhigh) retourne l'intégrale que vous voulez.