Create a PID tunning function and also wrap it to matlab compatibility.

It must heritage the tf class as answer and print in the format Kp + Ki * 1/s + Kd * s instead num / den format.

Originally posted by @kingaza in #139 (comment)