Lab3_main

clc;
y = dlmread('y5.txt');

c1 = 0.1; c2 = 0.3; c3 = 0.2; c4 = 0.08;
m1 = 12; m2 = 21; m3 = 18;
b = [c1; c4; m2;];

eps = 1e-7;
t = 0.2;
ok = true;

disp ("The value of pi is:")

while ok
    A = [0 1 0 0 0 0;
        ((-b(1) - c2) / m1) 0 (c2 / m1) 0 0 0;
        0 0 0 1 0 0;
        (c2/b(3)) 0 ((-c2 -c3) / b(3)) 0 (c3 / b(3)) 0;
        0 0 0 0 0 1;
        0 0 (c3 / m3) 0 ((-c3 - b(2)) / m3) 0;];

    y1 = SolveDiff(A, y(:, 1));
    U = FindU(A, y1, b(1), c2, c3, b(2) , m1, b(3), m3);

    deltaB = CalcDelta(U, y1, y);

    b = b + deltaB;
    b
    dif = 0;
    for i = 1:251
        dif = dif + (y(:, i) - y1(:, i))' * (y(:, i) - y1(:, i));
    end;

    if(dif < eps)
        ok = false;
    end;


end;

b