A = [0.0, 0.0];
B = [4.0, 0.0];
C = [2.0, 3.0];
P = 100.0;  % kN, C'de dusey asagi yuk

R_By = P * 2 / 4;
R_Ay = P - R_By;

uCA = (A - C) / norm(A - C);
uCB = (B - C) / norm(B - C);
A_mat = [uCA', uCB'];
b_vec = [0; P];
F = A_mat \ b_vec;
F_AC = F(1); F_BC = F(2);

uAC = (C - A) / norm(C - A);
uAB = (B - A) / norm(B - A);
F_AB = -(F_AC * uAC(1)) / uAB(1);

fprintf('RAy = %.1f kN, RBy = %.1f kN\n', R_Ay, R_By);
fprintf('F_AC = %.2f kN, F_BC = %.2f kN\n', F_AC, F_BC);
fprintf('F_AB = %.2f kN\n', F_AB);

% Bagimsiz dogrulama: B dugumu dengesi
uBC = (C - B) / norm(C - B);
uBA = (A - B) / norm(A - B);
checkX = F_BC*uBC(1) + F_AB*uBA(1);
checkY = F_BC*uBC(2) + F_AB*uBA(2) + R_By;
fprintf('B dugumu kontrolu: Fx=%.6f, Fy=%.6f\n', checkX, checkY);
