> `assign`(N, 20); 1
 

20 (1)
 

> `assign`(W, mul(`+`(z, `-`(k)), k = 1 .. N)); 1
 

`*`(`+`(z, `-`(1)), `*`(`+`(z, `-`(2)), `*`(`+`(z, `-`(3)), `*`(`+`(z, `-`(4)), `*`(`+`(z, `-`(5)), `*`(`+`(z, `-`(6)), `*`(`+`(z, `-`(7)), `*`(`+`(z, `-`(8)), `*`(`+`(z, `-`(9)), `*`(`+`(z, `-`(10)),...
`*`(`+`(z, `-`(1)), `*`(`+`(z, `-`(2)), `*`(`+`(z, `-`(3)), `*`(`+`(z, `-`(4)), `*`(`+`(z, `-`(5)), `*`(`+`(z, `-`(6)), `*`(`+`(z, `-`(7)), `*`(`+`(z, `-`(8)), `*`(`+`(z, `-`(9)), `*`(`+`(z, `-`(10)),...
(2)
 

> `assign`(p, expand(W)); 1
 

`+`(`*`(8037811822645051776, `*`(`^`(z, 4))), `-`(`*`(8752948036761600000, `*`(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(z, 2))), `*`(40171771630, `*`(`^`(z, 14))), `-`(`*`(16722808...
`+`(`*`(8037811822645051776, `*`(`^`(z, 4))), `-`(`*`(8752948036761600000, `*`(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(z, 2))), `*`(40171771630, `*`(`^`(z, 14))), `-`(`*`(16722808...
`+`(`*`(8037811822645051776, `*`(`^`(z, 4))), `-`(`*`(8752948036761600000, `*`(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(z, 2))), `*`(40171771630, `*`(`^`(z, 14))), `-`(`*`(16722808...
`+`(`*`(8037811822645051776, `*`(`^`(z, 4))), `-`(`*`(8752948036761600000, `*`(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(z, 2))), `*`(40171771630, `*`(`^`(z, 14))), `-`(`*`(16722808...
`+`(`*`(8037811822645051776, `*`(`^`(z, 4))), `-`(`*`(8752948036761600000, `*`(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(z, 2))), `*`(40171771630, `*`(`^`(z, 14))), `-`(`*`(16722808...
(3)
 

> `assign`(B, map(abs, p)); 1
 

`+`(`*`(8037811822645051776, `*`(`^`(abs(z), 4))), `*`(8752948036761600000, `*`(abs(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(abs(z), 2))), `*`(40171771630, `*`(`^`(abs(z), 14))), `...
`+`(`*`(8037811822645051776, `*`(`^`(abs(z), 4))), `*`(8752948036761600000, `*`(abs(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(abs(z), 2))), `*`(40171771630, `*`(`^`(abs(z), 14))), `...
`+`(`*`(8037811822645051776, `*`(`^`(abs(z), 4))), `*`(8752948036761600000, `*`(abs(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(abs(z), 2))), `*`(40171771630, `*`(`^`(abs(z), 14))), `...
`+`(`*`(8037811822645051776, `*`(`^`(abs(z), 4))), `*`(8752948036761600000, `*`(abs(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(abs(z), 2))), `*`(40171771630, `*`(`^`(abs(z), 14))), `...
`+`(`*`(8037811822645051776, `*`(`^`(abs(z), 4))), `*`(8752948036761600000, `*`(abs(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(abs(z), 2))), `*`(40171771630, `*`(`^`(abs(z), 14))), `...
`+`(`*`(8037811822645051776, `*`(`^`(abs(z), 4))), `*`(8752948036761600000, `*`(abs(z))), 2432902008176640000, `*`(13803759753640704000, `*`(`^`(abs(z), 2))), `*`(40171771630, `*`(`^`(abs(z), 14))), `...
(4)
 

> `assign`(Digits, `+`(N, 5)); 1
 

25 (5)
 

> plots[logplot](`/`(`*`(B), `*`(abs(p))), z = 0 .. `+`(N, 1), colour = BLACK, numpoints = 601)
 

Plot_2d
 

> plots[logplot]([seq([r, `/`(`*`(eval(B, z = r)), `*`(abs(`*`(r, `*`(eval(diff(p, z), z = r))))))], r = 1 .. N)], style = POINT, colour = BLACK, symbol = CROSS); 1
plots[logplot]([seq([r, `/`(`*`(eval(B, z = r)), `*`(abs(`*`(r, `*`(eval(diff(p, z), z = r))))))], r = 1 .. N)], style = POINT, colour = BLACK, symbol = CROSS); 1
 

Plot_2d
 

> plots[logplot]([seq([r, `/`(`*`(eval(W, z = r)), `*`(abs(`*`(r, `*`(eval(diff(p, z), z = r))))))], r = 1 .. N)], style = POINT, colour = BLACK, symbol = CROSS)
plots[logplot]([seq([r, `/`(`*`(eval(W, z = r)), `*`(abs(`*`(r, `*`(eval(diff(p, z), z = r))))))], r = 1 .. N)], style = POINT, colour = BLACK, symbol = CROSS)
 

Plot_2d
 

> [seq(coeff(p, z, `+`(N, `-`(k))), k = 0 .. N)]; 1
 

[1, -210, 20615, -1256850, 53327946, -1672280820, 40171771630, -756111184500, 11310276995381, -135585182899530, 1307535010540395, -10142299865511450, 63030812099294896, -311333643161390640, 1206647803...
[1, -210, 20615, -1256850, 53327946, -1672280820, 40171771630, -756111184500, 11310276995381, -135585182899530, 1307535010540395, -10142299865511450, 63030812099294896, -311333643161390640, 1206647803...
[1, -210, 20615, -1256850, 53327946, -1672280820, 40171771630, -756111184500, 11310276995381, -135585182899530, 1307535010540395, -10142299865511450, 63030812099294896, -311333643161390640, 1206647803...
[1, -210, 20615, -1256850, 53327946, -1672280820, 40171771630, -756111184500, 11310276995381, -135585182899530, 1307535010540395, -10142299865511450, 63030812099294896, -311333643161390640, 1206647803...
(6)
 

> CodeGeneration[Matlab](%); 1
 

cg = [1 -210 20615 -1256850 53327946 -1672280820 40171771630 -756111184500 11310276995381 -135585182899530 1307535010540395 -10142299865511450 63030812099294896 -311333643161390640 1206647803780373360 -3599979517947607200 8037811822645051776 -12870931245150988800 13803759753640704000 -8752948036761600000 2432902008176640000];
 

>