kill(all)$globalsolve:true$ /*vyčistíme si paměť a nastavíme řešič*/

rovnice rovnováhy tuhého tělesa

RR:[
R_x+N_2=0,                          /*směr x*/
R_y+N_1-F-N_3=0,                    /*směr y*/
N_1*2*a+N_2*a-N_3*2*a-F*a=0         /*kolem z*/
];

fyzikální rovnice pro jednotlivé deformovatelné pruty

Δ_1:N_1*a/E/S+α*ΔT*a;
Δ_2:N_2*a/E/S;
Δ_3:N_3*a/E/S;

deformační podmínky

DP:[
'(Δ_1=-Δ_3),
'(Δ_1=2*Δ_2)
];

dosazením z FR do DP za jednotlivá Δ dostaneme:

DP:ratsimp(''DP);

řešení soustavy RR (3 rovnice) a DP (2 rovnice) pro neznámé reakce Rx a Ry a síly v prutech N1, N2 a N3:

linsolve(append(RR,DP),[R_x,R_y,N_1,N_2,N_3]);

vyjádření řešení v jednotlivých prutech

'(N_1)=N_1;             /*zápis z pohledu Maximy čuňárna, ale vypadá to hezky...*/
'(N_2)=N_2;
'(N_3)=N_3;