Measurement of the electroweak production of Zγ and two jets in proton-proton collisions at s=13  TeV and constraints on anomalous quartic gauge couplings