A double Jacobi series is introduced to approximate functions of two independent variables. It is then applied to analysing simultaneous linear distributed parameter systems. The solution for the coefficient matrices can be obtained directly from a Kronecker product formula. The method is algebraic and computer oriented. One illustrative example is given for demonstration. Very satisfactory results are obtained, owing to the rapid convergence of the shifted Jacobi series.