kill(all);

definisco lunghezza totale albero

l : a + b;

equazioni di equilibrio alla rotazione attorno ad A e attorno a B

eqrAz : +RA*0 -F*a -C -P*(a+xi*b) + RB*l = 0 ;
eqrBz : -RA*l +F*b -C +P*(1-xi)*b + RB*0 = 0 ; 

risolvo la lista di equazioni (aka. sistema) in una lista di incognite

lista : [ ele1 ,ele2, ele3 ];

lista[3] restituisce ele3

sols : linsolve( [ eqrAz , eqrBz ] , [ RA , RB ] );

assegno a RA e RB le rispettive espressioni secondo sols
(uso assegnazione multipla)

[RA ,RB] : ev(  [RA ,RB] , sols );

Tratto AC: momento flettente

MfAC : RA*sAC;

Tratto AC: energia potenziale elastica

UAC : integrate( MfAC^2/2/E/J , sAC , 0 , a );