% Výpočet deformace mostovky zatížené ve středu bodovou silou

% uzly
nd = [0 0; 10 0; 20 0; 30 0; 40 0; 5 8; 15 8; 25 8; 35 8];

% elementy
el = [1 2; 2 3; 3 4; 4 5; ...
      6 7; 7 8; 8 9; 1 6; ...
      6 2; 2 7; 7 3; 3 8; ...
      8 4; 4 9; 9 5];

% materiálové parametry
mat = 500 * ones(1,15);

% působící síla - uzel 3, Fx = 0, Fy = -10
force = [3 0 -10];

% nulové posuvy na krajích
fix = [1 1; 1 2; 5 2];

% vlastní výpočet
u = mkp(el, nd, mat, force, fix)