Добавлено решение СЛАУ

+ вывод локальной матрицы жёсткости

MemMAPR-MKE.cpp

@@ -4,9 +4,9 @@ using namespace Eigen;
 
 int main() {
     std::ifstream file("matrix.txt");
-    int N = 20;
-    int upper = 6, lower = 1;
-    double A = 1, B = 1, C = 1, L = upper - lower, du_dx = 1.;
+    int N = 3;
+    int upper = 1, lower = 0;
+    double A = 1, B = 0, C = 1, L = upper - lower;
 
     double dx = L / N;
 
@@ -15,10 +15,11 @@ int main() {
     local(0, 1) = -A / dx + B / 2;
     local(1, 0) = -A / dx - B / 2;
     local(1, 1) = A / dx + B / 2;
+    std::cout << "Local matrix:\n"  << local << std::endl;
 
     VectorXd local_load(2);
-    local_load(0) = C * dx / 2 - A * du_dx; 
-    local_load(1) = C * dx / 2 + A * du_dx;
+    local_load(0) = C * dx / 2; 
+    local_load(1) = C * dx / 2;
 
     // Ансамблированные матрица и вектор
     MatrixXd ansamb = MatrixXd::Zero(N + 1, N + 1);
@@ -72,5 +73,9 @@ int main() {
     std::cout << "Modified matrix:\n" << ansamb << std::endl;
     std::cout << "Modified load vector:\n" << global_load << std::endl;
 
+    VectorXd solution = ansamb.fullPivLu().solve(global_load);
+    std::cout << "\nSolution:" << std::endl;
+    std::cout << solution << std::endl;
+
     return 0;
 }