diff --git a/.DS_Store b/.DS_Store
index b55547b..9595c9b 100644
Binary files a/.DS_Store and b/.DS_Store differ
diff --git a/RunSimulationKF.m b/RunSimulationKF.m
index ad02681..9419569 100644
--- a/RunSimulationKF.m
+++ b/RunSimulationKF.m
@@ -4,7 +4,7 @@
 mWz = 1200;
 v0 = 12/3.6;
 s0 = 0;
-Ts = 0.015;
+Ts = 0.010;
 Pres = 1/1000*[5.7/771 0 1.6]; %Strahl formula for m/s velocity
 i = -5/1000; %Gradient, uphill positive
 ssigma = 1e-1; % Slip at 68% of mumax, e.g. dry: 1%, wet: 10%
diff --git a/TractionCurve/TractionCurve.m b/TractionCurve/TractionCurve.m
new file mode 100644
index 0000000..2d7d39d
--- /dev/null
+++ b/TractionCurve/TractionCurve.m
@@ -0,0 +1,34 @@
+clear all, close all, %clc
+
+I = 50; %A
+U = 50; %V
+
+mu = 0.4;
+m = 600;
+mW = 1200;
+Fmax = mu*m*10;
+
+P = U*I;
+
+v = linspace(0, 15/3.6);
+
+F = P./v;
+Fres = 10*(m+mW)*1/1000*(5.7/771*(3.6*v).^2 +1.6)+10*(m+mW)*[0.01;0.02];
+
+Ind = find(F > Fmax);
+F(end) = 0;
+
+F(Ind) = Fmax;
+
+L = 2;
+plot(3.6*v, F, 'LineWidth', L)
+hold on;
+plot(3.6*v, Fres, 'LineWidth', L)
+ts = ['Traction curve, $P =$' num2str(U*I) ' W, $\mu =$' num2str(mu)  ', $m_L =$' num2str(m) ' kg, $m_{Wz} =$' num2str(mW) ' kg'];
+title(ts ,'interpreter','latex')
+legend('F_{t/b}', 'F_{R, 1%}', 'F_{R, 2%}')
+xlim([0 16])
+xlabel('$v /\left(kmh^{-1}\right)$','interpreter','latex')
+ylim([0 1.1*Fmax])
+ylabel('$F$/N','interpreter','latex')
+grid on
diff --git a/TractionCurve/TractionCurve.m~ b/TractionCurve/TractionCurve.m~
new file mode 100644
index 0000000..acd6229
--- /dev/null
+++ b/TractionCurve/TractionCurve.m~
@@ -0,0 +1,33 @@
+clear all, close all, clc
+
+I = 60; %A
+U = 50; %V
+
+mu = 0.3;
+m = 650;
+mW = 1200;
+Fmax = mu*m*10;
+
+P = U*I;
+
+v = linspace(0, 15/3.6);
+
+F = P./v;
+Fres = 10*(m+mW)*1/1000*(5.7/771*(3.6*v).^2 +1.6)+10*(m+mW)*[0.01;0.02];
+
+I = find(F > Fmax);
+
+F(I) = Fmax;
+
+L = 2;
+plot(3.6*v, F, 'LineWidth', L)
+hold on;
+plot(3.6*v, Fres, 'LineWidth', L)
+ts = ['Traction curve, $\mu =$' num2str(mu)  ', $m_L =$' num2str(m)];
+title(ts ,'interpreter','latex')
+legend('F_{t/b}', 'F_{R
+xlim([0 1.1*15])
+xlabel('v/kmh^{-1}')
+ylim([0 1.1*Fmax])
+ylabel('F/N')
+grid on
diff --git a/TractionCurve/TractionCurve.png b/TractionCurve/TractionCurve.png
new file mode 100644
index 0000000..74d4351
Binary files /dev/null and b/TractionCurve/TractionCurve.png differ