3523/CH4/EX4.10.3/Ex4_3.sce


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

//Example 3// Ch 4
clc;
clear;
close;
// given data
a=4*1e4;
b=15*1e5;
lb=0;
ub=0.0005;
i=integrate('(a-b*sqrt(x))','x',lb,ub)
as=exp(i);//Avalanche size
printf('Avalache size %f',as)