800/CH2/EX2.4/2_4.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14

//clear//
clc
clear
exec("2.4data.sci");


VCSTR = FA0*X(7)*(1/-rAat);
VPFR = FA0*inttrap(X,p)
disp("VCSTR =")
disp(VCSTR)
disp("dm^3")
disp("VPFR =")
disp(VPFR)
disp ("dm^3")