function Quiz4_1 %Trapezoidal

%Input data
a = 0 ;
b = 1;
n = 224;
h = (b-a)/n;
x = a:h:b;
f = sin(x.^2);

%weight=(1,2,2,2,2,...,2,2,2,1)
w(1:n+1)=2;
w(1) = w(n+1) = 1;

%Sum
s = (h/2)*w*f'

end

%s =         0.310270096421712
%exact = 0.310268301723381
%|exact-s| = 1.79469833100354e-06 < e-5