summaryrefslogtreecommitdiff
path: root/1535/CH1/EX1.3/Ch01Ex3.sci
blob: 8337ec8da0ddbd8a66afa084a79f1da1f4d4b42a (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
// Scilab Code Ex1.3 : Page-24 (2010)
A1 = 40;    // First amplitude of oscillation, cm
An_plus_1 = 4;    // Amplitude after 100 oscillations, cm
n = 100;    // Number of oscillations
T = 2.5;    // Time period of oscillations, s
t = T/4;    // Time taken to reach the first amplitude from the mean position, s
// Now A1 = x0*exp(-lambda*t) and An_plus_1 = x0*exp(-lambda*(t+nT))
// A1/An_plus_1 = exp(n*lambda*T), solving for lambda
lambda = log(A1/An_plus_1)/(n*T);    // Damping constant. per sec
printf("\nDamping constant = %3.2e per sec", lambda);

// Result 
// Damping constant = 9.21e-003 per sec