diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp
index 7b0ffa5e5c1d548e73bf410484edc88dbfe30a54..0fb6058a12e75ce7055b8a3d3cc92b33651bfafe 100644
--- a/src/scheme/FiniteVolumesDiffusion.hpp
+++ b/src/scheme/FiniteVolumesDiffusion.hpp
@@ -216,7 +216,7 @@ private:
 
     const Kokkos::View<const double*>& Vl = m_mesh_data.Vl();
     const Kokkos::View<const double*>& Vj = m_mesh_data.Vj();
-
+    
     Kokkos::parallel_for(m_mesh.numberOfFaces(), KOKKOS_LAMBDA(const int& l) {
         Rd sum = zero;
 	double sum2 = 0.;
@@ -431,7 +431,7 @@ public:
 	
 	// test k non cst pour diff pure
 	uj[j] += std::exp(-t)*(dt*inv_mj[j])*Vj(j)*(std::sin(pi*xj[j][0])*(pi*pi*xj[j][0]-1.) - std::cos(xj[j][0]*pi)*pi); 
-	Ej[j] -= ((pi*pi*0.5*(std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0])-std::cos(pi*xj[j][0])*std::cos(pi*xj[j][0])) + pi*pi*pi*std::cos(xj[j][0])*std::sin(xj[j][0]))*(std::exp(-2.*t)-1.) - pi*0.5*std::exp(-2.*t)*std::cos(xj[j][0])*std::sin(xj[j][0]) + (1.+ xj[j][0])*( (2*pi*pi*pi*std::cos(xj[j][0])*std::sin(xj[j][0])-4*std::pow(pi,4)*(std::sin(xj[j][0])*std::sin(xj[j][0])-std::cos(xj[j][0])*std::cos(xj[j][0])))*(std::exp(-2.*t)-1.) + pi*pi*0.5*std::exp(-2.*t)*(std::sin(xj[j][0])*std::sin(xj[j][0])-std::cos(xj[j][0])*std::cos(xj[j][0]))))*(dt*inv_mj[j])*Vj(j);
+	Ej[j] -= ((pi*pi*0.5*(std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0])-std::cos(pi*xj[j][0])*std::cos(pi*xj[j][0])) + xj[j][0]*pi*pi*pi*std::cos(pi*xj[j][0])*std::sin(pi*xj[j][0]))*(std::exp(-2.*t)-1.) - pi*0.5*std::exp(-2.*t)*std::cos(pi*xj[j][0])*std::sin(pi*xj[j][0]) + (1.+xj[j][0])*((3.*pi*pi*pi*std::cos(pi*xj[j][0])*std::sin(pi*xj[j][0])-xj[j][0]*pi*pi*pi*pi*(std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0])-std::cos(pi*xj[j][0])*std::cos(pi*xj[j][0])))*(std::exp(-2.*t)-1.)+0.5*pi*pi*std::exp(-2.*t)*(std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0])-std::cos(pi*xj[j][0])*std::cos(pi*xj[j][0]))))*(dt*inv_mj[j])*Vj(j);
 
 	// ajout second membre pour kidder (k = 0.5)
 	//Ej[j] -= (dt*inv_mj[j])*Vj(j)*((0.5*t*t)/(((50./9.)-t*t)*((50./9.)-t*t)));
diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp
index 334450d066c67756d6d033f36785f05ae56cf88e..b1750d2d5f39dd0a6b961a3402dc260c57275578 100644
--- a/src/scheme/FiniteVolumesEulerUnknowns.hpp
+++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp
@@ -459,7 +459,7 @@ public:
 	m_Tj[j] = 2 - 0.5*std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0]); // k = x
       });
 
-    // Conditions aux bords de Dirichlet sur u et k
+    // Conditions aux bords de Dirichlet sur T et nu
     
     m_TL[0] = 2.;
     m_TR[0] = 2.;