Bessoin d'aide intégration sous matlab

Bonjour,



mon problème consiste à calculer moyennant l'outil matlab:

deltaG=(1/Mprime)*intégrale de 0.000026 à M de (R×T×Ln(a) dM)

avec:

a=[0.01, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95];


M=[0.000026, 0.000404, 0.000496, 0.000564, 0.000621, 0.000674, 0.000724, 0.000773, 0.000822, 0.000872, 0.000925, 0.000981, 0.001040, 0.001108, 0.001185, 0.001255, 0.001385, 0.001530, 0.001788, 0.002106]


Mprime=[0, 0.000144, 0.000236, 0.000304, 0.000361, 0.000414, 0.000464, 0.000513, 0.000562, 0.000612, 0.000665, 0.000721, 0.000780, 0.000848, 0.000925, 0.000995, 0.001125, 0.001270, 0.001528, 0.001846]

R=8.314;

T=304.15;

puisque je suis entrain de vérifier le travail de quelqu'un d'autre, je dois trouvé comme résultats:

deltaG=[inf, -9.513, -8.395, -7.708, -7.183, -6.748, -6.373, -6.035, -5.725, -5.434, -5.151, -4.877, -4.614, -4.339, -4.060, -3.772, -3.466, -3.125, -2.734, -2.225]

est ce qu'il y a quelqu'un qui peut m'aider et merci d'avance.