Generally, PID controller is not suitable to control multi variable system because it is very difficult to tune the PID gains. However, this paper shows that it is not hard to tune the PID gains if we can find a Jacobian matrix of the system. The Jacobian matrix expresses the ratio of output variations according to input variations. It is possible to adjust the input values in order to reduce the output error using the Jacobian. When the colt function is composed of error related terms, the gradient approach can tune the PID gains to minimize the function. In simulation, a hydrofoil catamaran with two inputs and two outputs is applied as a multi variable system. We can easily get the multi variable PID controller by the proposed teaming method. When the controller is compared with LQR controller, the performance is as good as that of LQR controller with a modeling equation.