function [c, ceq] = beamConstraintsWithDisc(x)

x = beamMapVariables(x);


[c, ceq] = beamConstraints(x);


