softmatterlab · LauraPerezG · Jul 1, 2019 · Jul 10, 2019 · Jul 10, 2019 · Jul 10, 2019
diff --git a/.DS_Store b/.DS_Store
diff --git a/forma1d.m → 1d_trap/forma1d.m b/forma1d.m → 1d_trap/forma1d.m
diff --git a/forma1d.mat → 1d_trap/forma1d.mat b/forma1d.mat → 1d_trap/forma1d.mat
diff --git a/test_forma1d.m → 1d_trap/test_forma1d.m b/test_forma1d.m → 1d_trap/test_forma1d.m
@@ -20,7 +20,7 @@
 
 load forma1d.mat
 
-[k_forma1d, D_forma1d] = forma1d(s-s0, Dt, gamma);
+[k_forma1d, D_forma1d] = forma1d(x-x0, Dt, gamma);
 
 disp(['FORMA 1D results:'])
 disp(['k* = ' num2str( k_forma1d ) ' N/m'])

diff --git a/2d_forcefield/.DS_Store b/2d_forcefield/.DS_Store
diff --git a/BB_index_filter.m → 2d_forcefield/BB_index_filter.m b/BB_index_filter.m → 2d_forcefield/BB_index_filter.m
diff --git a/FORMA2D_F0_J0.m → 2d_forcefield/FORMA2D_F0_J0.m b/FORMA2D_F0_J0.m → 2d_forcefield/FORMA2D_F0_J0.m
diff --git a/doublebeam_dataxvsdx.mat.zip → 2d_forcefield/doublebeam_dataxvsdx.mat b/doublebeam_dataxvsdx.mat.zip → 2d_forcefield/doublebeam_dataxvsdx.mat
diff --git a/forma_2D_field.m → 2d_forcefield/forma_2D_field.m b/forma_2D_field.m → 2d_forcefield/forma_2D_field.m
diff --git a/2d_trap/.DS_Store b/2d_trap/.DS_Store
diff --git a/forma2d.m → 2d_trap/forma2d.m b/forma2d.m → 2d_trap/forma2d.m
diff --git a/forma2d.mat → 2d_trap/forma2d.mat b/forma2d.mat → 2d_trap/forma2d.mat
diff --git a/test_forma2d.m → 2d_trap/test_forma2d.m b/test_forma2d.m → 2d_trap/test_forma2d.m
@@ -19,8 +19,8 @@
 clc
 
 load forma2d.mat
-x=p-p0;
-y=q-q0;
+x=x-x0;
+y=y-y0;
 [k_forma2d, Omega_forma2d, theta_forma2d, D_forma2d] = forma2d([x,y], Dt, gamma);
 
 disp(['FORMA 2D results:'])

diff --git a/README.md b/README.md
@@ -9,12 +9,9 @@ Laura Pérez García, Jaime Donlucas Pérez, Giorgio Volpe, Alejandro V. Arzola,
 Nature Communications 9, 5166 (2018)
 https://doi.org/10.1038/s41467-018-07437-x
 
-# 1. Requirements to implement FORMA codes in Matlab
-
-The data points of your trajectories should be in in a way that the origin is in the center of the trap, i. e. where the $F=0$
-
 
 http://softmatterlab.org
 
 
+![double_well](https://user-images.githubusercontent.com/47359519/60442443-5246d500-9c19-11e9-9fcc-77b35dfe46bc.jpg)
 
diff --git a/test_forma2d.mat b/test_forma2d.mat
-Original file line number
+Diff line change
@@ Expand Up / @@ -19,8 +19,8 @@ @@
     clc
     load forma2d.mat
-    x=p-p0;
-    y=q-q0;
+    x=x-x0;
+    y=y-y0;
     [k_forma2d, Omega_forma2d, theta_forma2d, D_forma2d] = forma2d([x,y], Dt, gamma);
     disp(['FORMA 2D results:'])
@@ Expand Down @@