diff --git a/src/field.c b/src/field.c
index 18a18c2b7d7adccc68129336ea17e4879798085c..71cf15b27f91a7b75e616d5a1b0ae3f902a4fbc3 100644
--- a/src/field.c
+++ b/src/field.c
@@ -1,12 +1,22 @@
-#include "field.h"
+#include "physics.h"
 #include <math.h>
+#include <stdbool.h>
 #include <stdlib.h>
 
+#include "physics.h"
+#include "../utils/utils.h"
+
 
 // Compute E*qP/norm(qP)
 // Return false if norm(qP) < eps
-bool compute_e(charge_t c, vec2 p, double eps, vec2 *e){
- return EXIT_SUCCESS;
+/// qP = vectoriel(P-q)
+/// Compute the vector value of the field generated by a charge at a given point in space.
+bool compute_e(charge_t c, vec2 p, double eps, vec2 *e) {
+	vec2 relative_position = vec2_sub(p, c.pos);
+	if (vec2_norm(relative_position) < eps) return false;
+	double field = K * c.q / vec2_norm_sqr(relative_position);
+	*e = vec2_mul(field, vec2_normalize(relative_position));
+	return true;
 }
 
 // Compute the normalized sum of Ei*qiP/norm(qiP)
@@ -19,13 +29,12 @@ bool compute_total_normalized_e(charge_t *charges, int num_charges, vec2 p, doub
 // starting from pos0.
 // Returns false if pos0 is not a valid position
 // (for example if pos0 is too close to a charge).
-static bool draw_field_line(struct gfx_context_t *ctxt, charge_t *charges, int num_charges, double dx, vec2 pos0, double x0, double x1, double y0, double y1){
- return EXIT_SUCCESS;
+static bool draw_field_line(struct gfx_context_t *ctxt, charge_t *charges, int num_charges, double dx, vec2 pos0, double x0, double x1, double y0, double y1) {
+	return EXIT_SUCCESS;
 }
 
 // Draw all the charges
 // A circle with minus sign for negative charges
 // A circle with a plus sign for positive charges
-static void draw_charges(struct gfx_context_t *context, charge_t *charges, int num_charges, double x0, double x1, double y0, double y1){
-
+static void draw_charges(struct gfx_context_t *context, charge_t *charges, int num_charges, double x0, double x1, double y0, double y1) {
 }
diff --git a/src/field.h b/src/field.h
index 8c1a1626c9dd9d50e6256c40d00a4721d997d2a5..eac1d4eb869076510343e9f4f39c0a5c7d8ce54f 100644
--- a/src/field.h
+++ b/src/field.h
@@ -5,11 +5,10 @@
 #include "../utils/utils.h"
 #include <stdio.h>
 
-typedef struct charge_t
-{
-    double q;
-    vec2 pos;
-} charge;
+typedef struct _charge_t {
+	double q;
+	vec2 pos;
+} charge_t;
 
 // Compute E*qP/norm(qP)
 // Return false if norm(qP) < eps