We say that a system {zmF(z)}m=0∞ is a Beurling system if F is an outer function. Beurling’s approximation theorem asserts that if F is an outer function from H2(D) then the system {zmF(z)}m=0∞ is complete in the space H2(D). We prove that a Beurling system with F∈Hp(D),1≤p<∞ is an M−bases in Hp(D) with an explicit dual system. Any function f∈Hp(D),1≤p<∞ can be expanded as a series by the system {zmF(z)}m=0∞. For different methods of summation we characterize outer functions F for which the expansion converges to f. Related results for weighted Hardy spaces in the unit disc are studied. Particularly we prove Rosenblum’s hypothesis.