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

rovnice rovnováhy tuhého tělesa

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

fyzikální rovnice pro jednotlivé deformovatelné pruty

Δ_1:N_1*l_1/E/S;
Δ_2:N_2*l_2/E/S+α*l_2*ΔT;

deformační podmínka

DP:'(-Δ_1/a=Δ_2/2/a);

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

DP:[ratsimp(''DP)];

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

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

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;