solving a system of eqns

I need help seeing the solutions to this problem. It is giving me a matrix answer or a huge fraction. Any help would be appreciated

%12.16
% For one of the binary systems listed in the preceding table, based on
% eq 10.5 and the Wilson eqn, make the following calculations for
% (1)Propanol/(2)water
%a. Bubl P: t=60 deg C, x1=.3
%b. Dew P: t=60 deg C, y1=.3
%c. P,T-flash t=60 deg C, P=.5(Pbubble+Pdew), z1=.3
%d. If an azeotrope exists at t=60 deg C, find Paz and x1_az and y1+az

clear
clc

%Propanol
P_sat1=exp(16.1154-(3483.67/(60+205.807)));
V1=74.05 ;

%Water
P_sat2=exp(16.3872-(3358.79/(60+212.041)));
V2=18.07;

%Given
A12=291.27;
A21=1448.01;
R=8.314;
T=60+273;

lambda12=V2/V1*exp(-A12/(R*T))
lambda21=V1/V2*exp(-A21/(R*T))

%Part A
x1=.3
x2=1-x1

gamma1=exp((x2*(lambda12/(x1+x2*lambda12)-(lambda21/(x2+x1*lambda21)))))/(x1+x2*lambda12)
gamma2=exp((-x1*(lambda12/(x1+x2*lambda12)-(lambda21/(x2+x1*lambda21)))))/(x2+x1*lambda21)

syms y1 y2 P_b
y=solve(y1*P_b==x1*gamma1*P_sat1, y2*P_b==x2*gamma2*P_sat2, y1+y2==1)

**I also tried

[y1 y2 P_b]=solve(y1*P_b==x1*gamma1*P_sat1, y2*P_b==x2*gamma2*P_sat2, y1+y2==1, y1,y2,P_b)