% Ca³kowanie funkcji. Wyznaczanie zmiennoœci funkcji ca³kuj¹cej
clear;
clc;
F=inline('sin(x).*x');
Calk=quad(F,0,20);

F1=@(x) sin(x);
Calk1=quad(F1,0,20);




% Ca³kowanie w przedziale

x=[0:0.1:20];
dl=length(x);
Q=zeros(dl-1,1);
suma=Q;
for i=1:dl-1
    Q(i)=quad(F1,x(i),x(i+1));
    suma(i+1)=suma(i)+Q(i);
end;

plot(x,F1(x),'r-',x,suma,'b*');
