Shear stress in a rectangular beam follows a parabolic distribution, being zero at the top and bottom surfaces (free surfaces) and maximum at the neutral axis. This is derived from the shear stress formula τ = VQ/(Ib), where Q varies quadratically across the section.