4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

w.0 26.2202 :

σ.E 3.6177 :

N.p0 0.6689 :

N.E 0.3111 kN * :

u 4.8549 mm * - :

dl 5.244 m * :

E.y 210 kN mm 2^/ * :

A.0 86 mm 2^* :

i.0 6.9282 :

N.i Fon u 12 N.i N.E / + (1 N.i N.E / + (2^/ w.0 2 i.0 * / (2^* + (N.i dl * E.y A.0 * / * - :

?

x# Fon x 1 N * roots #

Chord method (Ioan prgm):

a 1 N * :

b 1000103^* N * - :

fa a Fon :

fb b Fon :

fa 0.0049 m * -

fb 0.2852 m *

iter 0 :

niter 300 :

iter niter ≤ (a b - abs 1 N * > (& fa fb * abs 0 mm 2^* > (& iter iter 1 + : fa a Fon : fb b Fon : c a b a - (fa fa fb - / * + : fc c Fon : fa fc * 0 mm 2^* ≥ a c : b c : if 61 line while

a kN 17.8336 -

b kN 17.8336 -

c kN 17.8336 -

c Fon 6.55 1019 -^* m *

iter 5

N.p c :

N.p kN 17.83 -

Chord method as function (Ioan prgm):

x Fct a b Sol fa a Fct : fb b Fct : niter 300 : iter 0 : iter niter ≤ (a b - abs 1 N * > (& fa fb * abs 0 mm 2^* > (& iter iter 1 + : fa a Fct : fb b Fct : c a b a - (fa fa fb - / * + : fc c Fct : fa fc * 0 mm 2^* ≥ a c : b c : if 61 line while c iter 1 * kN * 21 sys 61 line :

a 1 N * :

b 1000 kN * - :

N.pp N.i Fon a b Sol :

N.pp 1 el Fon 7.3751 1017 -^* m *

N.pp kN 17.8336 - 621 sys

Broyden's method - solving nonlinear system of eqiations,

Arguments:
f(1) - vector function returning vector of size n (f(x)≡0)
x0  - initial estimates, vector of size n
ε - tolerance (small number)
maxiter - naximal number of iterations (if exceeded - method not converging) 1 f x0 ε maxiter broyden next 1 : iter 1 : n x0 length : B n identity : next h B x0 f * - eval : h min abs x0 min ε + abs / ε < (h max abs x0 min ε + abs / ε < (& next 0 : x x0 h + : z x f x0 f - eval : B B B z * h - (h transpose * B * h transpose B * z * (1 el / - eval : iter iter 1 + : iter maxiter > next 0 : next 1 : x0 x : 21 line if 51 line if 21 line while konv iter maxiter ≥ not converging iter if : x konv 21 mat 71 line :

a 1 N * :

b 1000 kN * - :

ε 1020 -^:

maxiter 100 :

N.ppp x# Fon b ε maxiter broyden :

N.ppp #

1 f a b ε bisect b f 0 < a f 0 < res no zero in interval : b1 a : a1 b : 21 line if 11 line a f 0 > res no zero in interval : a1 a : b1 b : 21 line if 11 line if res IsDefined n 0 : res b1 a1 + (2 / : n 0 : res f abs ε > res f a1 f * 0 ≤ b1 res : a1 res : if res b1 a1 + (2 / : n n 1 + : 31 line while 31 line if res n 21 mat 31 line :

ε 108 -^:

a 1 N * :

b 1000 kN * - :

res x# Fon a b 108 -^bisect :

res 5105^* kg m * s 2^/ * - 021 mat

res 1 el Fon 0.14 m *