tuning transfer function for given dynamic performances

I want to find a second order transfer function with a non minimum phase zero z=36.6 which has 2% overshooting and a 2% settling time of 0.2s. I created a tunable transfer function but I don’t know how to find the values for the tunable parameters w and xi that allows the performances I want. Is there an automatic way to find them ? I tried to sample the variables and then chose the best function, but I think that this is not the best way to do.

 w = realp('w',25);
 xi = realp('xi',0.8);
 z = 36.6; 

 G = tf(w^2*[-1/z 1],[1 2*xi*w w^2]);