Path: news.mathworks.com!not-for-mail
From: "nouha " <nouha@voila.fr>
Newsgroups: comp.soft-sys.matlab
Subject: probleme: integrate a vecter
Date: Sat, 27 Nov 2010 10:25:06 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 28
Message-ID: <icqma2$3q5$1@fred.mathworks.com>
Reply-To: "nouha " <nouha@voila.fr>
NNTP-Posting-Host: webapp-02-blr.mathworks.com
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: fred.mathworks.com 1290853506 3909 172.30.248.37 (27 Nov 2010 10:25:06 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Sat, 27 Nov 2010 10:25:06 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 2624892
Xref: news.mathworks.com comp.soft-sys.matlab:690412

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.