(17:03) gp > N=26; x='x+O('x^N); Vec(x/(1-2^(1/2)*x)^2 * ((1+2^(1/2)*x)/(1-2^(1/2)*x))^(3/(2^(3/2))-1))
%1 = [1, 3.0000000000000000000000000000000000000, 6.5000000000000000000000000000000000000, 12.500000000000000000000000000000000000, 22.375000000000000000000000000000000000, 38.425000000000000000000000000000000000, 63.962500000000000000000000000000000001, 104.26250000000000000000000000000000000, 167.02343750000000000000000000000000000, 264.19947916666666666666666666666666667, 413.30671875000000000000000000000000000, 641.11897253787878787878787878787878788, 986.89318063446969696969696969696969697, 1509.9825252221736596736596736596736597, 2297.3540452451194638694638694638694638, 3479.4358594933712121212121212121212120, 5247.1023141452460300116550116550116547, 7884.8309508947270177738927738927738924, 11808.343120106279896318958818958818958, 17634.137131279919282334989571831677093, 26261.806809904547684988166073692389480, 39019.960949689059662525431439905124111, 57844.517385675785323957072798280932247, 85584.858949683656541219176723064282863, 126387.14214006202771556654268694489985]
(17:04) gp >