scilab/modules/differential_equations/tests/unit_tests/intg.dia.ref

   1 // to check that intg works
   2 function y=f(x),y=x*sin(30*x)/sqrt(1-((x/(2*%pi))^2)),endfunction
   3 exact=-2.5432596188;
   4 I=intg(0,2*%pi,f);
   5 if abs(exact-I) > 1 then bugmes();quit;end