  Non-local interactions are assumed for the deuteron with form factors
  $g_C(k)= j_0(b_1k)\,j_0(b_2k)$ for the central part responsible for the
S-state and
  $g_T(k)= j_1(b_1k)\,j_1(b_2k)$ for the tensor part responsible for the
D-state, where $b_1$ and $b_2$ are range parameters for the proton and neutron,
respectively. The analytically obtained wavefunctions in coordinate space have
different forms in three different regions. The inner most region is between
$r=0$ and $r=b_2-b_1$ (assuming $b_2>b_1$), followed by a region between
$r=b_2-b_1$ and $b_1+b_2$ and finally the region $r>b_1+b_2$. The resulting
wavefunctions and their derivatives are found to be continuous at the
boundaries. The ensuing calculations are simplified by setting $b_1=b_2=b$.
Good agreement is obtained for the quadrupole moment. Neutron-proton scattering
calculations will follow in the next communication.
