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

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

Originally posted by@kingaza in#139 (comment)