The Maple notebook is chapt2_5_1.mw but you are warned that it may download as a text file.

Explorations of d'Alembert's Series

The command SquareRoot[x,n] computes the first  n terms of the binomial series expansion of a and compares it to the actual value (10-digit accuracy).

> SquareRoot := (x,n) -> [1+evalf(sum((-1)^(k-1)*(2*k-2)!/2^(2*k-1)/k!/(k-1)!*x^k,k=1..n),10),evalf(sqrt(1+x),10)];

> SquareRoot(200/199,100);

The next command computes  this function at ten values of n: m/10, 2m/10, ..., m. If we chooose a value of m that is not a multiple of 10, it computes these at the floors. This is useful so that the number of terms is not always even or always odd.

> SquareRootList := (x,m) -> [seq([floor(j*m/10),SquareRoot(x,floor(j*m/10))[1]],j=1..10)];

> SquareRootList(200/199,1007);

The last command plots these values.

> plots[listplot](SquareRootList(200/199,1007));