From 49a4876bb823730cd58d6ba151a084d9df8e01ca Mon Sep 17 00:00:00 2001
From: gergaud <gergaud@enseeiht.fr>
Date: Mon, 9 Nov 2020 23:03:09 +0100
Subject: [PATCH] ajout du fichier test-TP-Projet-Optinum.jl

---
 src/test-TP-Projet-Optinum.jl | 38 +++++++++++++++++++++++++++++++++++
 1 file changed, 38 insertions(+)
 create mode 100644 src/test-TP-Projet-Optinum.jl

diff --git a/src/test-TP-Projet-Optinum.jl b/src/test-TP-Projet-Optinum.jl
new file mode 100644
index 0000000..29c61fe
--- /dev/null
+++ b/src/test-TP-Projet-Optinum.jl
@@ -0,0 +1,38 @@
+#using Pkg; Pkg.add("LinearAlgebra"); Pkg.add("Markdown")
+# using Documenter
+using LinearAlgebra
+using Markdown                             # Pour que les docstrings en début des fonctions ne posent
+                                           # pas de soucis. Ces docstrings sont utiles pour générer 
+                                           # la documentation sous GitHub
+include("Algorithme_De_Newton.jl")
+
+# Affichage les sorties de l'algorithme des Régions de confiance
+function my_afficher_resultats(algo,nom_fct,point_init,xmin,fxmin,flag,sol_exacte,nbiters)
+	println("-------------------------------------------------------------------------")
+	printstyled("Résultats de : ",algo, " appliqué à ",nom_fct, " au point initial ", point_init, ":\n",bold=true,color=:blue)
+	println("  * xsol = ",xmin)
+	println("  * f(xsol) = ",fxmin)
+	println("  * nb_iters = ",nbiters)
+	println("  * flag = ",flag)
+	println("  * sol_exacte : ", sol_exacte)
+end
+
+# Fonction f0
+# -----------
+f0(x) =  sin(x)
+# la gradient de la fonction f0
+grad_f0(x) = cos(x)
+# la hessienne de la fonction f0
+hess_f0(x) = -sin(x)
+sol_exacte = -pi/2
+options = []
+
+x0 = sol_exacte
+xmin,f_min,flag,nb_iters = Algorithme_De_Newton(f0,grad_f0,hess_f0,x0,options)
+my_afficher_resultats("Newton","f0",x0,xmin,f_min,flag,sol_exacte,nb_iters)
+x0 = -pi/2+0.5
+xmin,f_min,flag,nb_iters = Algorithme_De_Newton(f0,grad_f0,hess_f0,x0,options)
+my_afficher_resultats("Newton","f0",x0,xmin,f_min,flag,sol_exacte,nb_iters)
+x0 = pi/2
+xmin,f_min,flag,nb_iters = Algorithme_De_Newton(f0,grad_f0,hess_f0,x0,options)
+my_afficher_resultats("Newton","f0",x0,xmin,f_min,flag,sol_exacte,nb_iters)