diff --git a/Comparaison_of_methods.rmd b/Comparaison_of_methods.rmd
index b4f38ad..dc66d70 100644
--- a/Comparaison_of_methods.rmd
+++ b/Comparaison_of_methods.rmd
@@ -285,23 +285,25 @@ ScoreDistribElisa <- function(lambda0, lambda1, T){
 ```
 
 ```{r}
-distrib_score_mc=ScoreDistribEmpiric(2,3,10000,T)
+distrib_score_mc = ScoreDistribEmpiric(2,3,10000,T)
+distrib_score_theo = ScoreDistribElisa(2,3,T)
 
-distrib_score_theo=ScoreDistribElisa(2,3,T)
+plot_graph_distrib_score <- function(distrib_score_theo, distrib_score_mc){
+    # length(distrib_score_mc[,2])
+    # length(distrib_score_theo[,2])
 
+    #diff_distrib_score=abs(distrib_score_mc[,2]-distrib_score_theo[,2])
 
-length(distrib_score_mc[,2])
-length(distrib_score_theo[,2])
+    #par(mfrow = c(1,2))
+    barplot(distrib_score_mc[,2],col="blue",axes=F)
+    mtext("Distribution of scores via Monte Carlo",side=1,line=2.5,col="blue")
+    axis(2, ylim=c(0,10))
+    par(new = T)
+    barplot(distrib_score_theo[,2],col="red",axes=F)
+    mtext("Distribution of scores using the theoretical method",side=1,line=4,col="red") 
+}
 
-#diff_distrib_score=abs(distrib_score_mc[,2]-distrib_score_theo[,2])
-
-#par(mfrow = c(1,2))
-barplot(distrib_score_mc[,2],col="blue",axes=F)
-mtext("Distribution des scores via Monte_Carlo",side=1,line=2.5,col="blue")
-axis(2, ylim=c(0,10))
-par(new = T)
-barplot(distrib_score_theo[,2],col="red",axes=F)
-mtext("Distribution des scores via la méthode théorique",side=1,line=4,col="red") 
+plot_graph_distrib_score(distrib_score_theo, distrib_score_mc)
 ```