Add C-part to local matrix for linear element

Solver.cpp

@@ -16,10 +16,10 @@ void Solver::Execute_Linear(double val1, double val2) {
     // Локальный вектор нагрузки
     VectorXd local_load(2);
 
-    local(0, 0) = -A / L - B / 2.;// +C * L / 3.;
-    local(0, 1) = A / L + B / 2.; // +C * L / 6.;
-    local(1, 0) = A / L - B / 2.; // +C * L / 6.;
-    local(1, 1) = -A / L + B / 2.;// +C * L / 3.;
+    local(0, 0) = -A / L - B / 2. + C * L / 3.;
+    local(0, 1) = A / L + B / 2. + C * L / 6.;
+    local(1, 0) = A / L - B / 2. + C * L / 6.;
+    local(1, 1) = -A / L + B / 2. + C * L / 3.;
 
     local_load(0) = -D * L / 2.;
     local_load(1) = -D * L / 2.;