(%i1) list(n):=makelist([k/n,(1+1/n)^k],k,0,n);

(%i2) list(10);

(%i3) wxplot2d([discrete, list(10)], [style, points]);

(%i4) wxplot2d([discrete, list(10)]);

(%i5) wxplot2d([[discrete, list(10)],exp(x)],[x,0,1],[style,points,lines],
[legend,"methode d'Euler","y=exp(x)"],[xlabel, "n = 10"]);

(%i6) wxplot2d([[discrete, list(100)],exp(x)],[x,0,1],[style,points,lines],
[legend,"methode d'Euler","y=exp(x)"],[xlabel, "n = 100"]);

(%i7) wxplot2d([[discrete, list(100)],exp(x)],[x,0,1],[style,[points,1,1,1],lines],
[legend,"methode d'Euler","y=exp(x)"],[xlabel, "n = 100"]);

(%i8) float(second(last(list(100))));

(%i9)