[T1,T2]=solve(alpha*G-h1*(T1-Tinf)-sigma*epsilon*(T1^4-Te^4)-(k(T1-T2)/t),h2*(T2-Tinf)+sigma*epsilon*(T2^4-Te2^4)-(k(T1-T2)/t));
Or alternately,
[T1,T2]=solve(subs('alpha*G-h1*(T1-Tinf)-sigma*epsilon*(T1^4-Te^4)=(k(T1-T2)/t)'),subs('h2*(T2-Tinf)+sigma*epsilon*(T2^4-Te2^4)=(k(T1-T2)/t)'))
Best Answer