%Grønn tekst er kommentarer
clc % clear screen
%PS: I dag anter vi alltid samme steglengde

%Funksjonen vår:
f = @(x) sin(x).*exp(x)./((x+1).^3)
fplot(f, [0,10])

h=0.01
xverdier = 0:h:10;
fverdier = f(xverdier);

subplot(2,1,1)
plot(xverdier, fverdier)

%Df er den deriverte av f
Df = deriver(xverdier, fverdier);

subplot(2,1,2)
plot(xverdier(1:length(xverdier)-1), Df, 'r')