clear all, clc
h=0.1;
fel=inf;
startX=0.3; slutX=10;
for i =1:2
    U=[0.1 ; 0.2];
    for x=startX:h:slutX-h
        U=U+h*fEXS710(x,U);
    end
    Uend(i)=U(1);
    h=h*0.5;
end
format long
disp(Uend)
diff(Uend)
    