It is because you’re trying to use the MarkerFaceAlpha property in Plot function, but this property is the part of scatter chart rather than plot function.
Just change the line
p1=plot(x,y,'.')
to
p1=scatter(x,y,'.')
And this should not produce any error