clear all; clc; close all
format long
syms a b l W mu P NA NB
list={a, b, l, W, mu};
listn={5, 24.5, 25, 25, 0.3};
rA=[a, 0, 0];
rB=[0, b, 0];
rC=[a/2, b/2, 0];
G=[0, -W, 0];
FC=[P, 0, 0];
FA=[-mu*NA, NA, 0];
FB=[NB, 0, 0];
R = G+FA+FB+FC;
MA=cross(rB-rA, FB)+cross(rC-rA, FC+G);
sol=solve(R(1),R(2),MA(3),'P, NA, NB');
NAs=sol.NA;
NBs=sol.NB;
Ps=sol.P;
NAn=eval(subs(NAs,list,listn));
NBn=eval(subs(NBs,list,listn));
Pn=eval(subs(Ps,list,listn));
fprintf('NA = %s = %g (lb)\n',char(NAs),NAn);
fprintf('NB = %s = %g (lb)\n',char(NBs),NBn);
fprintf('P = %s = %g (lb)\n',char(Ps),Pn);
NA = W = 25 (lb)
NB = W*(a-b*mu)/b = -2.39796 (lb)
P = -W*(-2*b*mu+a)/b = 9.89796 (lb)