! PEP 2012 - testing the T ~ sqrt(l/g) for a simple pendulum

! small ball
stat {18.57;18.93;18.88} d_s\mean	! in mm
d_s=d_s*0.1

l_s={96.5;48.0;24.0;12.0;9.6}
dls=l_s*0+0.5
T_s={1.9792;1.3986;0.9812;0.7201;0.6147}
dTs={0.0008;0.0002;0.0001;0.0001;0.0005}

! large ball
stat {24.43;24.29;24.37} d_l\mean
d_l=d_l*0.1

l_l={94.1;47.0;22.4;10.1}
dll=l_l*0+0.5
T_l={1.9625;1.3923;0.9837;0.7020}
dTl={0.0012;0.0004;0.0003;0.0003}

clear
defaults
scales 0 0 0 0 2 4
label\y `T, s'
label\x `Length of pendulum'

set pchar -13
ls=l_s+0.5*d_s
graph ls,T_s,dTs,dls

g=980.67	! in cm/s^2
scalar\vary A
A=1
fit T_s=A*sqrt(ls/g)
generate l 0,,100 1000
set pchar 0
graph\noax l,A*sqrt(l/g)

A_s=A

colour 2 2
set pchar -15
ll=l_l+0.5*d_l
graph\noax ll,T_l,dTl,dll
fit T_l=A*sqrt(ll/g)
set pchar 0
graph\noax l,A*sqrt(l/g)

=`small ball: A = '//rchar(A_s)
=`large ball: A = '//rchar(A)

colour 6 6
graph\noax l,2*Pi*sqrt(l/g)