Dividing both sides by $z_{n-1}$ gives
$$3b_n=1+\frac{1}{b_{n-1}}+z_{n-2}$$
Since $z_n\to 0$, $z_n>0$ for $n\geq 1$, and the solutions to
$$3x=1+\frac{1}{x}$$
are negative and positive respectively, it is sufficient to prove that $b_n$ converges to prove that
$$\lim_{n\to\infty} b_n= \frac{1+\sqrt{13}}{6}$$
Now, it is easy to prove from induction that $z_n<.9^n$ for all $n\geq 2$. We can manually check that this is the case for $n\in\{2,3,...,11\}$. Then for $n>11$ we have
$$z_n=\frac{1}{3}\left(z_{n-1}+z_{n-2}+z_{n-1}z_{n-2}\right)$$
$$<\frac{1}{3}\left(.9^{n-1}+.9^{n-2}+.9^{n-1}.9^{n-2}\right)$$
Since $n>11$ we have that $.9^n<\frac{1}{3}$. This then gives us
$$=\frac{1}{3}\left(\frac{.9^{n}}{.9}+\frac{.9^{n}}{.9^2}+\frac{.9^{n}.9^{n}}{.9^3}\right)$$
$$<\frac{1}{3}\left(\frac{.9^{n}}{.9}+\frac{.9^{n}}{.9^2}+\frac{.9^{n}}{3\cdot.9^3}\right)$$
$$=9^n\left(\frac{6130}{6561}\right)<.9^n$$
(this also provides a proof that $z_n\to 0$). Next, we will prove by induction that $\frac{7}{10}<b_n<\frac{10}{11}$ for $n\geq 6$. We can manually check that this is the case for $n\in\{6,7,...,14\}$. For $n>14$ we have
$$b_n=\frac{1}{3}+\frac{1}{3b_{n-1}}+\frac{z_{n-2}}{3}<\frac{1}{3}+\frac{1}{3\cdot \frac{7}{10}}+\frac{.9^{n-2}}{3}$$
$$=\frac{17}{21}+\frac{100}{243}.9^n$$
Since $n>14$ this implies $.9^n<\frac{23}{100}$ which gives us
$$<\frac{17}{21}+\frac{100}{243}\cdot \frac{23}{100}=\frac{1538}{1701}<\frac{10}{11}$$
For the lower bound we have
$$b_n=\frac{1}{3}+\frac{1}{3b_{n-1}}+\frac{z_{n-2}}{3}$$
$$>\frac{1}{3}+\frac{1}{3\cdot \frac{10}{11}}=\frac{7}{10}$$
as desired. In fact, we really will only need this lower bound for the rest of the proof. Define $c_n=|b_n-b_{n-1}|$. Then
$$c_n=b_n-b_{n-1}=\left|\frac{1}{3b_{n-1}}-\frac{1}{3b_{n-2}}+\frac{z_{n-2}-z_{n-3}}{3}\right|$$
$$\leq \left(\frac{1}{3b_{n-1}b_{n-2}}\right)|b_{n-1}-b_{n-2}|+\frac{z_{n-2}+z_{n-3}}{3}$$
$$=\left(\frac{1}{3b_{n-1}b_{n-2}}\right)c_{n-1}+\frac{z_{n-2}+z_{n-3}}{3}$$
Importantly, using our lower bound from above we have
$$\frac{1}{3b_{n-1}b_{n-2}}<\frac{1}{3\left(\frac{7}{10}\right)^2}=\frac{100}{147}<.7$$
Since $z_n<.9^n$ for all $n\geq 2$, we know that $n\geq 5$ we have
$$\frac{z_{n-2}+z_{n-3}}{3}<.91^{n-5}$$
Then for these $n\geq 5$ we have
$$c_n<.7 c_{n-1}+.91^{n-5}$$
Now, we will prove by induction that $c_n<.99^n$ for all $n\geq 5$. We can manually check that this is true for $n\in\{5,6,...,21\}$. For $n>21$ we have
$$c_n<.7 .99^{n-1}+.91^{n-5}=.99^n\left[\frac{70}{99}+\left(\frac{100}{91}\right)^5\cdot \left(\frac{91}{99}\right)^n\right]$$
However, since $n>21$ we have
$$\left(\frac{91}{99}\right)^n<\frac{180969322079}{990000000000}$$
which gives us
$$<.99^n\left[\frac{70}{99}+\left(\frac{100}{91}\right)^5\cdot \frac{180969322079}{990000000000}\right]=.99^n$$
as desired. We conclude that $|b_n-b_{n-1}|<.99^n$ (at least for $n\geq 5$). But at this point we are basically done! It is well known that this property implies the sequence in question is Cauchy (see here for example just replacing $.5$ with $.99$). Thus, $b_n$ converges and as shown above this implies
$$\lim_{n\to\infty}b_n=\frac{1+\sqrt{13}}{6}$$