A class of exact axisymmetric tokamak equilibria with sheared flows, parallel to the magnetic field is constructed. The solutions form nested magnetic surfaces and contrary to the standard procedure are non-separable functions with respect to the radial and vertical variables (R, z). Also a particular solution satisfies a recently established sufficient condition for linear MHD stability and is constructed to agree with the geometrical data of actual experiments.