You are supposed to be solving for vo/vs, yet in numerous places you have the ratio vo/ve2; isn't ve2 the same as vo?
You seem to have done something different from what you show. For example, you have two KCL equations for vc1:
v_{C_{1}} = (\frac{v_{E_{2}}}{βr_{e}} - \frac{v_{S}(R_{Th}||βr_{e})}{(R_{S} + R_{Th}||βr_{e})βr_{e}})\frac{1}{\frac{1}{R_{C}||R_{Th}} + \frac{1}{βr_{e}}}
and
v_{C_{1}} = \frac{v_{E_{2}}(\frac{1}{βr_{e}} + \frac{1}{r_{e}} + \frac{1}{R_{E}||R_{L}})}{\frac{1}{βr_{e}} + \frac{1}{r_{e}}}
and then you say that you set these two together and solved for vo/ve2. Of course, you meant that you solved for vo/vs.
The problem I see is that if you do in fact use these two equations, you don't get what you say you got. Setting these two equations, as you have shown them, together and solving for vo/vs gets this result:
\frac{v_{O}}{v_{E_{2}}} = \frac{(R_{Th}||βr_{e})(R_{C}||R_{Th})(\frac{βr_{e}}{R_{C}||R_{Th}} + 1)(\frac{1}{β} + 1)}{(\frac{1}{β} + 1 - (\frac{βr_{e}}{R_{C}||R_{Th}} + 1)(\frac{1}{β} + 1 + \frac{r_{e}}{R_{E}||R_{L}}))(R_{S} + R_{Th}||βr_{e})(r_{e}β + R_{C}||R_{Th})}
You'll notice that there is a missing β in the numerator. One wonders, how did you get the result you posted? I can't have come from the two KCL equations you posted for vc1.
The answer is that the first KCL equation you show has an extra β. It should be:
v_{C_{1}} = (\frac{v_{E_{2}}}{βr_{e}} - \frac{v_{S}(R_{Th}||βr_{e})}{(R_{S} + R_{Th}||βr_{e})r_{e}})\frac{1}{\frac{1}{R_{C}||R_{Th}} + \frac{1}{r_{e}}}
Somehow in bulling through all of this, you got a good result, even though your first KCL equation, as you have shown it here, has an error.
If you used the first KCL equation as you show it in post #44, your result for vo/vs would be off by a factor of 150.
The advantage of setting up a nodal solution as I showed is that you let the computer (or calculator) do all the mistake-prone algebra. The elements of the 4 KCL equations are individually very simple. You don't solve them individually; you just set them up, which is almost trivially easy.
They can be solved in various ways, but formatting them as an array gives a compact representation. Solving them as a linear system of 4 simultaneous equations gives all 4 node voltages at once. There is no "bulling through" required. The computer or calculator does all the algebra for you, without any mistakes, like this:
The answer is an array of the 4 node voltages, (vs, vb1, vc1, vo). Since the input vs was set to 1, the output voltage at vo, -90.22, is the gain Av.
You can instantly see that the gain from vs to vb1 is .62, and from vs to vc1 is -91.83
This problem was of moderate complexity. As soon as you have a few more nodes, probably with feedback, there will be no possibility of getting a mistake free solution doing all the algebra by hand. This method is what simulators like pSpice use.