function [c, ceq] = conFunEq(gp,gc,msEq,opt)
% global gp msEq
% s = tf('s');
% k = t(1);
% I = t(2);
% D = t(3);

%gc = k*((1-I-D)+I/s+D*s);
c   = [];
ceq = - msEq + ms(gp,gc,opt);

return