We present measurements of the elliptic (k to view the MathML source">v2), triangular (k to view the MathML source">v3) and quadrangular (k to view the MathML source">v4) anisotropic azimuthal flow over a wide range of pseudorapidities (k to view the MathML source">−3.5<η<5). The measurements are performed with Pb–Pb collisions at using the ALICE detector at the Large Hadron Collider (LHC). The flow harmonics are obtained using two- and four-particle correlations from nine different centrality intervals covering central to peripheral collisions. We find that the shape of k to view the MathML source">vn(η) is largely independent of centrality for the flow harmonics k to view the MathML source">n=2–4, however the higher harmonics fall off more steeply with increasing i8" class="mathmlsrc">i8.gif&_user=111111111&_pii=S0370269316303586&_rdoc=1&_issn=03702693&md5=8f072ec63cf8a21187020234d0a227ae" title="Click to view the MathML source">|η|. We assess the validity of extended longitudinal scaling of k to view the MathML source">v2 by comparing to lower energy measurements, and find that the higher harmonic flow coefficients are proportional to the charged particle densities at larger pseudorapidities. Finally, we compare our measurements to both hydrodynamical and transport models, and find they both have challenges when it comes to describing our data.