foo : 'diff (1/(1 + sqrt(1 - a*x^2)), x, 2);

bar : ev (foo, nouns);

integrate (x^k, x);

integrate (x^k, x);

plot2d (subst (a = 1, bar), [x, -2, 2], [y, -10, 10], [svg_file, "maxplot.svg"],[plot_format, gnuplot])$