function []=uppg2()
% Uppgift 2a)

A=[2 -3 7 -13 5;...
   1 1 4 11 -1;...
  -2 -3 4 -8 -3;...
  -1 4 -1 4 -3;...
   3 1 1 1 3];
b=[1 2 3 4 5]';
x=A\b;
display(x)
pause

% Uppgift 2b)

[t,x] = ode45(@f,[0 2], [5 10]);
plot(x(:,1),x(:,2)),shg %fasplott
pause
figure(2)
plot(t,x)

function dx=f(t,x)
dx=ones(2,1);
dx(1)=x(1)+x(2).*sin(t);
ds(2)=-x(1).*x(2);
