There are many functions having the property of orthogonality. Among them are Walsh functions, block-pulse functions, Chebyshev functions and Legendre functions. Recently, many scholars and researchers made new approaches to solving variational problems via orthogonal functions. This paper suggests a more simple and powerful orthogonal method—the Fourier series operational matrix for integration—which is derived and used to solve a variational problem.