Update implicite.cpp

src/implicite.cpp

@@ -60,9 +60,9 @@ void settab(Materiaux mat, int u0, double f)
     double coef= (mat.getlambda()*dt)/(mat.getc()*mat.getkho()*(dx*dx));
     double F2=16*f*f;
     vector<double> F(101,0);
-    F[10]=F2;
+    F[10]=F2/(mat.getc()*mat.getkho());
     for(int j=50;j<61;j++){
-        F[j]=F2*0.75;
+        F[j]=F2*0.75/(mat.getc()*mat.getkho());
     }
 
     vector<vector<double> > temp;