read\format\errskip ?1 (*) t,x
x=x/100		! cm -> m

clear
defaults
label\x `t, min'
label\y `x(t), m'
set pchar -5
set yleadz 1
set font triumf.2
scales 0 240 6 0 .85 4
set yvmax 0.8
set nsxinc 4
set nsyinc 2

graph t,x

scalar\vary a,P,x0,tau,drift,t0
x0=0.73
a=-0.18
P=10
tau=30
drift=0

t0=45
fit\weight\e2 (t>t0)*(t<145) x=a*exp(-(t-t0)/tau)*cos(2*Pi/P*(t-t0))+x0+drift*(t-t0)

tt=[t0:145]
set pchar 0
set colour 7
graph\noaxes tt a*exp(-(tt-t0)/tau)*cos(2*Pi/P*(tt-t0))+x0+drift*(tt-t0)

set cursor -1
set %xloc 30
set %yloc 20
text `T='//rchar(60*P,`F5.1')//`<pm>'//rchar(60*FIT$E2[2],`F3.1')//` s'

x0=0.45
a=0.35
P=10
tau=30
drift=0

t0=145
fit\weight\e2 (t>t0) x=a*exp(-(t-t0)/tau)*cos(2*Pi/P*(t-t0))+x0+drift*(t-t0)

tt=[t0:220]
set colour 6
graph\noaxes tt a*exp(-(tt-t0)/tau)*cos(2*Pi/P*(tt-t0))+x0+drift*(tt-t0)

set cursor -1
set %xloc 65
text `T='//rchar(60*P,`F5.1')//`<pm>'//rchar(60*FIT$E2[2],`F3.1')//` s'