function integral1
y=quad(@cos,0,2*pi);
fprintf('Integral z cos od 0 do 2*pi je: %.2f\n',y);
y=quad(@cos,0,pi);
fprintf('Integral z cos od 0 do pi je: %.2f\n',y);
y=quad(@cos,0,pi/2);
fprintf('Integral z cos od 0 do pi/2 je: %.2f\n',y);

x=0:0.1:2*pi;
for n=1:length(x)
    y(n)=quad(@cos,0,x(n));
end
plot(x,cos(x),x,y)
legend('cos(x)','integral z cos(x)')
grid    
end

