clear all;

Nt=1000;
T = linspace(0,3,Nt);
F = zeros(1,Nt);

delta=T(2)-T(1);

for it=1:Nt
    t=T(it);
    F(it)  = coeff(t);
    
end

%figure(1)
%plot(T,F,'*-b')

%Ith = 3^3/3+3;

I1=0;
for k=1:Nt-1
    I1 = I1 + delta*(F(k)+F(k+1))/2;
end

display(['to olokliroma tis coeff gia ' num2str(T(1)) '<t<' num2str(T(end)) ])
I1