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

rovnice rovnováhy pro 1xSN stav

RR:N_1+N_2+F=0;

fyzikální rovnice

Δ_1:N_1*l/E/S_1;
Δ_2:N_2*l/E/S_2;

mezní síla při změně stavu ze staticky určitého na 1xSN stav

linsolve(ev(Δ_1,N_1:-F_m)=-δ,F_m);

deformační podmínka pro 1xSN stav

DP:Δ_1+δ=Δ_2;

řešení soustavy rovnic RR a DP

linsolve([RR,DP],[N_1,N_2]);

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

N_1:if F<F_m then -F else N_1$
N_2:if F<F_m then 0 else N_2$

wxplot2d(ev([N_1,N_2],E:1,S_1:100,S_2:200,δ:1,l:100),[F,0,3],[legend,"N1","N2"])$