function [Edof,Ex,Ey,nelm,ndof,bc]=geom37

L=1.5;      % Längd på stänger
nelm=14;    % Definera totalt antal element 'Number of ELeMents'
nen=2;      % Definera antal noder per element 'Number of Elements Nodes'
ndof=18;    % Definera totalt antal frihetsgrader 'Number of Degrees Of Freedom'


% Definiera koordinaterna till varje nod
Coord= [0 	0		
       	L 	0		
        2*L 	0		
	3*L	0		
	4*L	0		
	0	L
	L	L
	2*L	L
	3*L	L];


% Definiera vilka frihetsgrader som hör till varje nod
% Notera att vi har två frihetsgrader per nod (förskjuting i x-led/ förskjutning i y-led)
for i=1:length(Coord)		
	Dof(i,:)=[i*2-1 i*2];
end

% Definera topologin
% Rad 1:Element 1 sitter ihop med frihetsgraderna 1   2   3   4	
% Rad 2:Element 2 sitter ihop med frihetsgraderna 3   4   5   6
% Rad 3:Element 3 sitter ihop med frihetsgraderna 7   8   5   6   
% o s v
Edof=[1	 1   2   3   4		
      2  3   4   5   6   
      3  7   8   5   6   
      4  7   8   9   10
      5  11  12  13  14
      6  13  14  15  16
      7  15  16  17  18
      8  3   4   13  14
      9  5   6   15  16
      10 7   8   17  18
      11 11  12  3   4
      12 13  14  5   6
      13 15  16  7   8
      14 17  18  9   10];


% Kommandot coordxtr tar ut koordinaterna för varje element.      
[Ex,Ey]=coordxtr(Edof,Coord,Dof,nen);


% Randvillkor: Eftersom strukturen är helt låst i två
% noder (med frihetsgrader 1 2 11 12) skall bc (boundary conditions)
% su ut enligt nedan
bc=[1   0
    2   0
    11  0      
    12  0];  
