function integralsin2x2
h=quad(@sin2x2,0,2*pi);
fprintf('Integral z sin2x2 od %.2f do %.2f je %.2f\n',0,2*pi,h);
h=quad(@sin2x2,3,7);
fprintf('Integral z sin2x2 od %.2f do %.2f je %.2f\n',3,7,h);
x=0:0.01:2*pi;
for k=1:length(x)
    h(k)=quad(@sin2x2,0,x(k));
end
plot(x,sin2x2(x),x,h)
legend('sin^2x^2','integral z sin^x^2')
end
function y=sin2x2(x)
y=sin(x.^2).^2;
end

