Thread Subject:

probleme: integrate a vecter

hey,



my problem is to calculate this integral tool matlab subject :

deltaG=(1/Mprime)*[int (R×T×Ln(a) dM)] "from 0.000026 to M"

with:

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;

since I'm currently checking the work of someone else, I have to found the following results:

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]

is there anyone who can help me and thank you in advance.

"nouha " <nouha@voila.fr> wrote in message <icqma2$3q5$1@fred.mathworks.com>...
> hey,
>
>
>
> my problem is to calculate this integral tool matlab subject :
>
> deltaG=(1/Mprime)*[int (R×T×Ln(a) dM)] "from 0.000026 to M"
>
> with:
>
> 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;
>
> since I'm currently checking the work of someone else, I have to found the following results:
>
> 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]
>
> is there anyone who can help me and thank you in advance.
- - - - - - - - - -
  This sounds like a problem for matlab's 'cumtrapz'. See the documentation for details.

Roger Stafford

think you so much for the solution