% Detta är en skriptfil
% Integralexempel

clear all, clf
disp(['Integralen av funktionen: '...
    ' 2.*sin(x)-cos(x)./(1+x.^2)']);
a = input('Ange undre integrationsgräns:  ');
b = input('Ange övre integrationsgräns:  ');

int = quadl(@fkn,a,b);
disp(['Integralvärdet blir:  '   num2str(int)]);

fplot(@fkn,[a,b])
grid, 
xlabel('x-axel'), ylabel('y-axel')
title('Integralexempel')