data1 = load('data/data_0.2/part2_position.txt');
data2 = load('data/data_0.3/part2_position.txt');
data3 = load('data/data_0.4/part2_position.txt');
data4 = load('data/data_0.5/part2_position.txt');
data5 = load('data/data_0.6/part2_position.txt');

t1 = data1(:,1);
x1 = data1(:,2);
t2 = data2(:,1);
x2 = data2(:,2);
t3 = data3(:,1);
x3 = data3(:,2);
t4 = data4(:,1);
x4 = data4(:,2);
t5 = data5(:,1);
x5 = data5(:,2);

numPeaks = [2,2,2,3,4];
p1 = findPeaks(data1,numPeaks(1));
p2 = findPeaks(data2,numPeaks(2));
p3 = findPeaks(data3,numPeaks(3));
p4 = findPeaks(data4,numPeaks(4));
p5 = findPeaks(data5,numPeaks(5));

myFit1 = fit(p1(:,1),p1(:,2),'exp1');
myFit2 = fit(p2(:,1),p2(:,2),'exp1');
myFit3 = fit(p3(:,1),p3(:,2),'exp1');
myFit4 = fit(p4(:,1),p4(:,2),'exp1');
myFit5 = fit(p5(:,1),p5(:,2),'exp1');

a1 = 23.84;
b1 = -0.8044;
a2 = 7.332;
b2 = -0.4692;
a3 = 12.41;
b3 = -0.3557;
a4 = 11.08;
b4 = -0.2756;
a5 = 12.87;
b5 = -0.2195;

f1 = a1 * exp(b1 * t1);
f2 = a2 * exp(b2 * t2);
f3 = a3 * exp(b3 * t3);
f4 = a4 * exp(b4 * t4);
f5 = a5 * exp(b5 * t5);

figure();
hold on;
plot(t1,x1);
plot(t1,f1);
hold off;

figure();
hold on;
plot(t2,x2);
plot(t2,f2);
hold off;

figure();
hold on;
plot(t3,x3);
plot(t3,f3);
hold off;

figure();
hold on;
plot(t4,x4);
plot(t4,f4);
hold off;

figure();
hold on;
plot(t5,x5);
plot(t5,f5);
hold off;