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

