%
%   Debugs application of SEM to bar elements
%


phi_undam = sin(2*pi * Nodes(:,2) / max(Nodes(:,2)) );

phi_dam = phi_undam;
phi_dam(23) = phi_dam(23) * 0.90;

interpopt = 'spline';
numdiv = 2;
mode_anal = 1;

DImat = dibars(Nodes,Bars,phi_undam,phi_dam,interpopt,numdiv,mode_anal);
