Question: problem in fsolve to solve some complicated systems of nonlinear algebraic equations

Hello Dear Maple Users and Experts

I am running this code for N=3, but fsoolve can not work and did not give me any output. Could you help me how can I get the result?

Actually, I got the result for N=2. Exact solution is a[0]=0, a[1]=1, a[2]=1, b[1]^2=sqrt(2) and b[1]^2+b[2]^2=sqrt(3). But, for N=3, I can not receive any results from fsolve

Here is my code

restart;
Digits := 20;
L := 1;
N := 3;
alpha := 1;
xexact := t -> t^sqrt(2) + t^sqrt(3);
f := simplify(fracdiff(t^sqrt(2), t, alpha)) + simplify(fracdiff(t^sqrt(3), t, alpha));
f := unapply(f, t);
xapp := a[0] + sum(a[j]*t^sum(b[i]^2, i = 1 .. j), j = 1 .. N);
xapp := unapply(xapp, t);
xfrac := sum(a[jj]*simplify(GAMMA(sum(b[ii]^2, ii = 1 .. jj) + 1)/GAMMA(sum(b[ii]^2, ii = 1 .. jj) + 1 - alpha))*t^(sum(b[ii]^2, ii = 1 .. jj) - alpha), jj = 1 .. N);
xfrac := unapply(xfrac, t);
xfrac1 := sum(a[jj]*simplify(sum(b[ii]^2, ii = 1 .. jj)^(alpha + 1)/(sum(b[ii]^2, ii = 1 .. jj) - alpha))*t^(sum(b[ii]^2, ii = 1 .. jj) - alpha), jj = 1 .. N);
xfrac1 := unapply(xfrac1, t);
S1 := {seq(evalf(xfrac(k/(2*N)*L)) - evalf(f(k/(2*N)*L)) = 0, k = 1 .. 2*N)};
S2 := {xapp(0) = 0};
S := S1 union S2;
sol := fsolve(S);