Plotting the root locus for changing PID parameters?

I am trying to write some code, which for random system shows it's the root locus plot, and how it changes for a varying P,I Or D constant.

I am pretty lost on how i should make an Matlab script which does this.