function dx = de_2(t,x)
% Indatafil för att lösa en andra ordningens differentialekvation:
% y" + ay' + by = h(t)
% där a och b och är lämpligt valda (konstanter).
% Vi skriver om ekvationen som ett system av första ordningens
% differentialekvationer: Sätt x1 = y och x2 = y. Då är 
%
% x1' = x2
% x2' = -a*x2 - b*x1 + h(t)
%
% Vi löser detta system genom att anropa Matlabfunktionen ode45
% (i kommandofönstret);
% 
% > [t,Y] = ode45(@de_2,[0 10], [1 0]);plot(t,Y(:,1))
% Utdata från ode45 kommer att innehålla en kolonnvektor som ger tidpunkter
% och en matris Y (2 kolonner) där första kolonnen innehåller beräknade
% värden för den sökta lösningen y till den ursprungliga diffekvationen.
% Den andra kolonnen innehåller beräknade värden för derivatan av lösningen y.
a=1.1; b = 3;

dx = zeros(2,1); % Vi initierar utdata som en kolonnvektor av längd 2.

dx(1) = x(2);
dx(2) = -a*x(2) - b*x(1) + t*exp(-t); % Här är funktionen  h(t) = t*exp(t).