% Laboration i FMM
%% Uppgift 1
clear all
load randwebs
load randinits
A = randwebscell{3};
x = randinitscell{3};
opts = odeset('NonNegative',1);
K = zeros(10);



tmax = 1001;
y= zeros(length(x), tmax);
y(:, 1) = x;
disp('hej')
for t = 1:tmax-1
    for i = 1:10
        K(i, :) = y(i, t).*A(i,:).*(y(i, t))';
    end
    for i = 1:10
        [tout,yout] = (ode45(@(t,y) sum(K(i,:)),[0 0.1],y(i, t),odeset('NonNegative',1))) ;
        y(i, t+1) = yout(end);
    end
end
plot(1:tmax, y)
% [r,t] = ode45(@(x, A) odefunc(x, A),[0,2],x(1))
% odefunc(x,A)
% function dxdt = odefunc(a,b)
%
% for i = 1:10
%      for j = 1:10
%          dxdt(i) = a(i).*a(j).*b(i,j);
%      end
% end
% end