function r=restrictions(vdoft) global A; ndof=max(size(vdoft)); r=[vdoft(1)-A]; end