function rij(a,b,c,d,R0,J0,N)



M=[a, b; c,d];

pom=[R0;J0];
R=R0;
J=J0;
for i=1:N
    pom=M*pom;
    R(i+1)=pom(1);
    J(i+1)=pom(2);
end

plot(0:N,R,0:N,J)
legend('Romeo', 'Julija')

[v,l]=eig(M);
A=v\[R0;J0];
JRn=A(1)*l(1,1)^N*v(:,1)+A(2)*l(2,2)^N*v(:,2)
R(end)
J(end)
    