We analyse inclusive photoproduction of polarized P-3(1) quarkonium in
the framework of QCD. To separate non-perturbative and perturbative p
arts in the density matrix of the produced quarkonium we use a method,
which is equivalent to the diagrammatic expansion widely used in anal
yzing deeply inelastic scattering. A systematic expansion in a small v
elocity upsilon, with which a heavy quark moves inside the quarkonium
in its rest frame, is performed for the non-perturbative parts, and th
ey are expressed as matrix elements in non-relativistic QCD, At the le
ading order of upsilon there are four matrix elements representing non
-perturbative physics. The perturbative parts are calculated at the le
ading order of coupling constants. Some numerical results, especially
numerical results for HERA, are given.