function integral_sin2x2
h=quad(@sin2x2,0,pi/2);
fprintf('Integral z sin2x2 od %.2f do %.2f je %.2f\n',0,pi/2,h);
h=quad(@sin2x2,6,18);
fprintf('Integral z sin2x2 od %.2f do %.2f je %.2f\n',6,18,h);
x=0:0.1:2*pi;
for n=1:length(x)
    h(n)=quad(@sin2x2,0,x(n));
end
plot(x,sin2x2(x),x,h)
legend('sin^2x^2','integral z sin^2x^2')
end
function y=sin2x2(x)
y=sin(x.^2).^2;
end

