(bind-func nintegrate "integrate using Simpsons rule - N must be odd" (lambda (f:[double,double]* xmin xmax) (nintegrate:[double,[double,double]*,double,double,i64]* f xmin xmax 501)))