  In earlier papers [3,4,5,6] Gursey et al. showed development of a bilocal
baryon-meson field from two quark-antiquark fields. The Hamiltonian in the case
of vanishing quark masses was shown to have a very good agreement with
experiments [5]. The theory for vanishing mass was solved using Confluent
Hypergeometric functions [6].
  In this paper I construct the normalized wave function for the spin-free
Hamiltonian with light quark masses (only up to the first order of the mass of
quark). I develop the new kind of special function theory in mathematics that
generalize all existing theories of Confluent Hypergeometric types. I call it
the 'Grand Confluent Hypergeometric (GCH) Function.' My solution produces
previously unknown extra "hidden" radial quantum numbers relevant for
description of supersymmetry and for generating new mass formulas.
  This paper is 1st out of 10 in series "Special functions and three term
recurrence formula (3TRF)". See section 6 for all the papers in the series. The
next paper in the series describes generalization of three term recurrence
relation in linear ordinary differential equations and its applications [8].
