]
%% Mesh generation

[Mesh.elements, Mesh.nodes] = CreateMesh2(nelx, nely, nelz, sets);
%% Material properties

MP.E = 200e9;      	% Elastic modulus (homogeneous, linear, isotropic material)
MP.nu = 0.3;       	% Poisson ratio
%% General Settings

[sets.nel, sets.nxe]  = size(Mesh.elements);    % Number of elements in the mesh & Number of nodes per element
[sets.nnod, sets.dim] = size(Mesh.nodes);       % Number of nodes in the mesh & Space dimension
sets.dxn = 3;                                   % Numbe