Constraint satisfaction programs where each constraint depends on a constant number of variables have the following property: The randomized algorithm that guesses an assignment uniformly at random satisfies an expected constant fraction of the constraints. Combining constructions from interactive proof systems with harmonic analysis over finite Abelian groups, Håstad [J. ACM, 48 (2001), pp. 798--859] showed that for several constraint satisfaction programs this naive algorithm is essentially the best possible unless $P = NP$. While most of the predicates analyzed by Håstad depend on a small number of variables, Samorodnitsky and Trevisan [{\it Proceedings of the 32nd Annual ACM Symposium on Theory of Computing, Portland, OR, 2000, pp. 191--199] recently extended Håstad's result to predicates depending on an arbitrarily large, but still constant, number of Boolean variables. We combine ideas from these two constructions and prove that there exists a large class of predicates on finite non-Boolean domains such that for predicates in the class, the naive randomized algorithm that guesses a solution uniformly is essentially the best possible unless P = NP$. As a corollary, we show that it is NP-hard to approximate the Maximum k-CSP problem over domains with size d within $d^{k-2k^{1/2}} - \epsilon$, for every constant $\epsilon > 0$, unless P = NP. This lower bound extends the previously known bound for the case d = 2 and matches well with the best known upper bound, dk-1, of Serna, Trevisan, and Xhafa [Proceedings of the 15th Annual Symposium on Theoretical Aspects of Computer Science, Lecture Notes in Comput. Sci. 1373, M. Morvan, C. Meinel, and D. Krob, eds., Springer-Verlag, Berlin, 1998, pp. 488--498].