float(integrate(x*sin(x),x,-%pi,%pi));
float(integrate(x*sin(2*x),x,-%pi,%pi));

load("simpsn");
f(x):=x*sin(3*x);
simpson(f,-float(%pi),float(%pi),100);
float(integrate(x*sin(3*x),x,-%pi,%pi));