a ( i ) = ((((( 2 k ( h- (((x ( i ) 2) + (L ^2)) A.5) + dely) ( cos (( atan ( L /x (i ))))))
)- (.5 c +jj) fr)- ( 2 wts .5 (sin ((( 2 atan ( L /x ( i )))))))) / ( m (.5 c +
JJ)));
av(i)= (a(i- 1)+a(i))/2;
t(i)= (-v(i- 1 )+((v(i 1))2 +(4 av(i) (x (i)- x(i 1))))A .5)/(2 av(
i ));
v(i ) =v(i- 1)+av(i)*t(i);
T(i)=T(i- 1)+t(i);
end


% Calculation of stiffness
mint = min (t (i ));
X =T(i)/mint;
mingt = GT / X;
forj = 3 :(((L-xd)/.01)+1)
S(j) t(j)/mint;
b(j)= S(j)* mingt;
vg(j)= (x(j)-x(j-1))/b(j);
ag(j) =(vg(j)-vg(j- 1))/b(j);
end


% Stiffness calculation of the spring depending on acceleration (ag)
AV ( 2 )= 0;
forq =3: (((L -xd)/.01)+ 1)
KKK (q)= 2 ((( h- ((x ( q) 2) + (L A2)) A.5)+ dely) *( cos (( atan (L/ x(q))))));
KKKK(q) = (( .5 c +jj ) m ag ( q )) + (( .5 c +jj) fr) + (2 wts ( cos((( atan (L
/ x (q )))))));
KKKKK ( q ) = KKKK ( q ) / KKK( q );
AV (q) =KKKKK(q) + AV(q- 1 )
end
SS= AV(q)/q