This note is concerned with a semi-analytical method for the solution of 2-D Helmholtz equation in unit square. The method uses orthogonal functions to project the problem down to finite dimensional space. After the projection, the problem simplifies to that of obtaining solutions for second order constant coefficient differential equations which can be done analytically. Numerical results indicate that the method is particularly useful for very high wave numbers.