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

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

return