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

rovnice rovnováhy styčníku

RR:[
-N_1*sin(α)+N_3*sin(β)=0,       /*směr x*/
N_1*cos(α)+N_2+N_3*cos(β)-F=0   /*směr y*/
];

fyzikální rovnice pro jednotlivé deformovatelné pruty

Δ_1:N_1*(h/cos(α))/E/S;
Δ_2:N_2*h/E/S;
Δ_3:N_3*(h/cos(β))/E/S;

deformační podmínky

DP:[
'(Δ_1=u*sin(α)-v*cos(α)),
'(Δ_2=0-v),
'(Δ_3=-u*sin(β)-v*cos(β))
];

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

DP:ratsimp(''DP);

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

linsolve(append(RR,DP),[u,v,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;