друге
> eq1:=f1 (p1) * c1 * c2 + f2 (p1) * c1 + f3 (p1) * c2 + f4 (p1) - (f1 (p2) * c1 * c2 + f2 (p2) * c1 + f3 ( p2) * c2 + f4 (p2));
збираємо щодо с1 і с2
> eq2:=collect (eq1, {c1, c2}, distributed);
висловлюємо с1
> c1:=solve (eq2, c1);
підставляємо с1 в перше рівняння
> eq3:=f1 (p1) * c1 * c2 + f2 (p1) * c1 + f3 (p1) * c2 + f4 (p1);
> eq4:=simplify (eq3);
беремо чисельник
> eq5:=f1(p1)*c2^2*f3(p2)+f1(p1)*c2*f4(p2)+f2(p1)*f3(p2)*c2+f2(p1)*f4(p2)-f3(p1)*c2^2*f1(p2)-f3(p1)*c2*f2(p2)-f4(p1)*f1(p2)*c2-f4(p1)*f2(p2);
збираємо щодо с2, отримуємо квадратне рівняння
> eq6:=collect (eq5, {c2}, distributed);
> W:=2; g:=9.8; d1:=1.5; d2:=0.2; J1:=0.25; J:=0.5; l1:=0.1; l2:=0.2; l:=0.3; p1:=0.01130833612; p2:=2.211662868;
> f1 (p1):=(- 2*l1*l^2*g*W*p1^4*l2*J1+l1*l^3*g*W*p1^2*l2^2*d2+l1*l^3*g^2*d1*p1^2*J1-l1^2*l^2*g*p1^4*J1*W+l1*l^3*g^2*p1^2*J1*d2+l1^3*l^3*g*d1*p1^2*W- 2*l1*l^4*g^2*d1*d2*l2-l1^2*l^4*g^2*d1*d2-g^2*l^4*l2^2*d1*d2-g*l^2*l2^2*W*p1^4*J1+g*l^3*l2^3*W*p1^2*d2+g^2*l^3*l2*d1*p1^2*J1+g^2*l^3*l2*p1^2*J1*d2+g*l^3*l2*d1*l1^2*p1^2*W)/l^4/g^2;
> f1 (p2):=(- 2*l1*l^2*g*W*p2^4*l2*J1+l1*l^3*g*W*p2^2*l2^2*d2+l1*l^3*g^2*d1*p2^2*J1-l1^2*l^2*g*p2^4*J1*W+l1*l^3*g^2*p2^2*J1*d2+l1^3*l^3*g*d1*p2^2*W- 2*l1*l^4*g^2*d1*d2*l2-l1^2*l^4*g^2*d1*d2-g^2*l^4*l2^2*d1*d2-g*l^2*l2^2*W*p2^4*J1+g*l^3*l2^3*W*p2^2*d2+g^2*l^3*l2*d1*p2^2*J1+g^2*l^3*l2*p2^2*J1*d2+g*l^3*l2*d1*l1^2*p2^2*W)/l^4/g^2;
> f2(p1):=(l1^2*l*p1^6*W^2*l2*J1-l1^2*l^2*p1^4*W^2*l2^2*d2-2*p1^4*J1*g*l^2*d1*l1^2*W-J^2*W^3*p1^4*g*l*l1+l1^3*l*p1^6*W^2*J1+J^2*W^2*p1^2*g^2*l^2*d1-J^2*W^3*p1^4*g*l*l2+J^2*W^2*p1^2*g^2*l^2*d2+l1^3*l^3*p1^2*W*d1*d2*g+p1^6*J1^2*g*l*l1*W-p1^4*J1^2*g^2*l^2*d2+p1^2*J1*g^2*l^3*d1*l1*d2+p1^6*J1^2*g*l*W*l2-l1^4*l^2*p1^4*W^2*d1+l1^2*l^3*p1^2*W*d1*g*d2*l2-l1^2*l^2*p1^4*W*J1*d2*g+p1^2*J1*g^2*l^3*d1*d2*l2-p1^4*J1^2*g^2*l^2*d1-p1^4*J1*g*l^2*W*l2^2*d2)/l^4/g^2;
> f2(p2):=(l1^2*l*p2^6*W^2*l2*J1-l1^2*l^2*p2^4*W^2*l2^2*d2-2*p2^4*J1*g*l^2*d1*l1^2*W-J^2*W^3*p2^4*g*l*l1+l1^3*l*p2^6*W^2*J1+J^2*W^2*p2^2*g^2*l^2*d1-J^2*W^3*p2^4*g*l*l2+J^2*W^2*p2^2*g^2*l^2*d2+l1^3*l^3*p2^2*W*d1*d2*g+p2^6*J1^2*g*l*l1*W-p2^4*J1^2*g^2*l^2*d2+p2^2*J1*g^2*l^3*d1*l1*d2+p2^6*J1^2*g*l*W*l2-l1^4*l^2*p2^4*W^2*d1+l1^2*l^3*p2^2*W*d1*g*d2*l2-l1^2*l^2*p2^4*W*J1*d2*g+p2^2*J1*g^2*l^3*d1*d2*l2-p2^4*J1^2*g^2*l^2*d1-p2^4*J1*g*l^2*W*l2^2*d2)/l^4/g^2;
> f3(p1):=(W^2*p1^6*l2^2*l*J1*l1-p1^4*J1^2*g^2*l^2*d2-J^2*W^3*p1^4*g*l*l2-p1^4*J1^2*g^2*l^2*d1+W^2*p1^6*l2^3*l*J1-W^2*p1^4*l2^2*l^2*d1*l1^2-2*p1^4*J1*g*l^2*W*l2^2*d2+J^2*W^2*p1^2*g^2*l^2*d1-J^2*W^3*p1^4*g*l*l1-W^2*p1^4*l2^4*l^2*d2+p1^6*J1^2*g*l*l1*W+p1^2*J1*g^2*l^3*d1*d2*l2+p1^6*J1^2*g*l*W*l2+J^2*W^2*p1^2*g^2*l^2*d2-W*p1^4*l2^2*l^2*d1*g*J1+W*p1^2*l2^3*l^3*d1*g*d2+W*p1^2*l2^2*l^3*d1*l1*d2*g+p1^2*J1*g^2*l^3*d1*l1*d2-p1^4*J1*g*l^2*d1*l1^2*W)/l^4/g^2;
> f3(p2):=(W^2*p2^6*l2^2*l*J1*l1-p2^4*J1^2*g^2*l^2*d2-J^2*W^3*p2^4*g*l*l2-p2^4*J1^2*g^2*l^2*d1+W^2*p2^6*l2^3*l*J1-W^2*p2^4*l2^2*l^2*d1*l1^2-2*p2^4*J1*g*l^2*W*l2^2*d2+J^2*W^2*p2^2*g^2*l^2*d1-J^2*W^3*p2^4*g*l*l1-W^2*p2^4*l2^4*l^2*d2+p2^6*J1^2*g*l*l1*W+p2^2*J1*g^2*l^3*d1*d2*l2+p2^6*J1^2*g*l*W*l2+J^2*W^2*p2^2*g^2*l^2*d2-W*p2^4*l2^2*l^2*d1*g*J1+W*p2^2*l2^3*l^3*d1*g*d2+W*p2^2*l2^2*l^3*d1*l1*d2*g+p2^2*J1*g^2*l^3*d1*l1*d2-p2^4*J1*g*l^2*d1*l1^2*W)/l^4/g^2;
> f4(p1):=(-J^2*W^3*p1^4*l1*d1*g*l-J^2*W^3*p1^4*l1*d2*g*l+W^4*p1^6*l2^2*J^2-W^2*p1^8*l2^2*J1^2-2*p1^4*J1*l1*W*d1...
