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

