function I=intpos(a,b,n)
% Calcola
%   int_a^b sin(x^2) dx
% usando il metodo dei trapezi con N partizioni.
% Esempio:
%   >> intpos(300, 400, 100000)

h = (b-a)/n;
x = a+h:h:b-h;
y = sin(x.^2);
J = mean([sin(a^2), sin(b^2)]);
I = (J + sum(y))*h;