#FIT_MAXITER = 2
set angles radians
s(x) = a * sin(2*pi*f*x - p ) + o + 1
f = 200000
a = 1
o = -1
p = 1

fit s(x) 'mcm05_ch1.dat' using 1:2 via p,a, o
fit s(x) 'mcm05_ch1.dat' using 1:2 via p,a,f, o
plot 'mcm05_ch1.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'mcm05_ch2.dat' using 1:2 via p,a, o
fit s(x) 'mcm05_ch2.dat' using 1:2 via p,a,f, o
plot 'mcm05_ch2.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'mcm04_ch1.dat' using 1:2 via p,a, o
fit s(x) 'mcm04_ch1.dat' using 1:2 via p,a,f, o
plot 'mcm04_ch1.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'mcm04_ch2.dat' using 1:2 via p,a, o
fit s(x) 'mcm04_ch2.dat' using 1:2 via p,a,f, o
plot 'mcm04_ch2.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'mcm03_ch1.dat' using 1:2 via p,a, o
fit s(x) 'mcm03_ch1.dat' using 1:2 via p,a,f, o
plot 'mcm03_ch1.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'mcm03_ch2.dat' using 1:2 via p,a
fit s(x) 'mcm03_ch2.dat' using 1:2 via p,a,f
plot 'mcm03_ch2.dat' using 1:2 with lines, s(x)
pause -1


fit s(x) 'waf02_ch1.dat' using 1:2 via p,a, o
fit s(x) 'waf02_ch1.dat' using 1:2 via p,a,f, o
plot 'waf02_ch1.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'waf02_ch2.dat' using 1:2 via p,a, o
fit s(x) 'waf02_ch2.dat' using 1:2 via p,a,f, o
plot 'waf02_ch2.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'mcm02_ch1.dat' using 1:2 via p,a, o
fit s(x) 'mcm02_ch1.dat' using 1:2 via p,a,f, o
plot 'mcm02_ch1.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) 'mcm02_ch2.dat' using 1:2 via p,a
fit s(x) 'mcm02_ch2.dat' using 1:2 via p,a,f
plot 'mcm02_ch2.dat' using 1:2 with lines, s(x)
pause -1

fit s(x) '1_left.dat' using 1:2 via p,a
#FIT_MAXITER = 10
fit s(x) '1_left.dat' using 1:2 via p,a,f
plot "1_left.dat" using 1:2 with linespoints, s(x)
pause -1

fit s(x) '2_right.dat' using 1:2 via p,a
#FIT_MAXITER = 10
fit s(x) '2_right.dat' using 1:2 via p,a,f
plot "2_right.dat" using 1:2 with linespoints, s(x)
pause -1

fit s(x) '2_left.dat' using 1:2 via p,a
#FIT_MAXITER = 10
fit s(x) '2_left.dat' using 1:2 via p,a,f
plot "2_left.dat" using 1:2 with linespoints, s(x)
pause -1

# DDS amplitude 4000
# wafer tester 2, left  a=0.658119 V
# wafer tester 2, right a=0.604414 V
# 3038.9641
# wafer tester 2, left a=0.482893 V
# 3308.99
# wafer tester 2, right a=0.480627 V
# DDS amplitude 3000
# mcm tester 2, ch1 a=0.865451 V
# mcm tester 2, ch2 a=0.825308 V
# DDS amplitude 3466
# mcm tester 3, ch1 a=0.997 V
# DDS amplitude 3635
# mcm tester 3, ch2 a=1.03905 V