function [a,b]=m36hw8(a0,b0)

a=zeros(1,100);
b=zeros(1,100);

a(1)=a0;
b(1)=b0;

for i=1:99
a(i+1)=a(i)+3*a(i)+4*b(i)-2;
b(i+1)=b(i)+2*a(i)-3*b(i)+1;
end


plot(a(1:20),'b*')
hold on
plot(b(1:20),'rd')
set(gca,'yscale','log')