Update README

README.md

@@ -69,4 +69,18 @@ Pour passer à la suite, le point courant devient le nouveau point et on recalcu
 
 ## Résultats
 
+Dans cette partie, nous allons détailler comment nous sommes arrivées à créer notre simulation de lignes de champ et pourquoi nous pensons que notre résultat est bonne représentation du phénomène physique. Ci-dessous, vous pouvez observer le résultat de notre simulation avec deux particules. La première particule est chargée positivement et placée à la position x=0.25 et y=0.5, la deuxième particule est chargée négativement et placé à la position x=0.75 et y=0.5. Il est important de noter que notre univers est de taille [0, 1] x [0, 1], le point (0, 0) étant en haut à gauche et le point (1, 1) en bas à droite de l'écran.
+
+![](assets/results_simulation_result.png)
+
+<center>Figure 2 - Premier résultat de la simulation</center>
+
+Si vous comparer l'image ci-dessus et l'image initial que nous avons mis dans l'introduction, la simulation ce comporte comme prédite, cependant nous allons créer un cas plus complexe afin de pouvoir être plus confiant et pouvoir dire que notre simulation  est correcte par rapport à la théorie. 
+
+Afin de pouvoir valider notre simulation, nous avons utilisé ce site internet https://academo.org/demos/electric-field-line-simulator qui nous permet de créer des prédictions et de les comparer à notre simulation. Le site academo.org est un site gratuit qui met à disposition des resources sur de sujets comme les mathématiques, la physique, etc. Nous avons estimé qu'il était suffisament de confiance afin de pouvoir se baser dessus.
+
+|                          Prédiction                          |                       Notre simulation                       |
+| :--------------------------: | :--------------------------: |
+|Particule 1 (Q = 1, inconnue), Particule 2 (Q = -5, inconnue)![](assets/results_simulation_pred_1.png)| Particule 1 (Q = 1, (0.25, 0.8)), Particule 2 (Q = -5, (0.85, 0.2))![](assets/results_simulation_res_1.png) |
+
 ## Conclusion

src/main.c

@@ -22,12 +22,12 @@ int main(int argc, char *argv[]) {
     double x1 = 1;
     double y1 = 1;
 
-    int num_charges = 4;
+    int num_charges = 2;
     charge_t *charges = (charge_t *)malloc(sizeof(charge_t) * num_charges);
-    charges[0] = charge_create(ELEMENTARY_CHARGE * 5, vector2_create(.25, .5));
-    charges[1] = charge_create(-ELEMENTARY_CHARGE, vector2_create(.75, .5));
-    charges[2] = charge_create(-ELEMENTARY_CHARGE * 4, vector2_create(.5, .25));
-    charges[3] = charge_create(ELEMENTARY_CHARGE * 5, vector2_create(.7, .75));
+    charges[0] = charge_create(ELEMENTARY_CHARGE, vector2_create(.25, .8));
+    charges[1] = charge_create(-ELEMENTARY_CHARGE * 5, vector2_create(.85, .2));
+    // charges[2] = charge_create(-ELEMENTARY_CHARGE * 4, vector2_create(.5, .25));
+    // charges[3] = charge_create(ELEMENTARY_CHARGE * 5, vector2_create(.7, .75));
 
     double dx = compute_delta_x(SCREEN_WIDTH, SCREEN_HEIGHT);