Question: How to extract data from implicit function with a parameter in MAPLE

I plotted number of functions y=y(x) with different parameter aa for given value of tau (see the code below) 

restart;
R0:=gamma+2*ln(2)+Re(Psi(1/2+(1+I)*tanh((1+I)*x)/((tau+0.25e-1*a)*y)))+ln(y);
tau := 9.975;
with(plots);
implicitplot({seq(R0, a = 1 .. 50)}, x = 0 .. 5, y = 0 .. 1, numpoints = 1000);

My question is how I can extract the values of y from these plots in the form of matrix 50×1000 (a×x) in data file?

Thanks in a advance.

Please Wait...