function integral1
h=quad(@cos,0,2*pi);
fprintf('Integral z cos od 0 do 2*pi: %.2f\n',h);
h=quad(@cos,0,pi/2);
fprintf('Integral z cos od 0 do pi/2: %.2f\n',h);
h=quad(@cos,0,pi);
fprintf('Integral z cos od 0 do pi: %.2f\n',h);
h=quad(@cos,0,3);
fprintf('Integral z cos od 0 do 3: %.2f\n',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(x)','integral z cos(x)')
end

