%% well_1

Pwf=150;
Pr=300;

mo=461997;
mg=56429;
mw=144599;

%GOR D
disp('Darcy')
kg=mg/(Pwf^2-Pr^2)
ko=mo/(Pwf-Pr)
kw=mw/(Pwf-Pr)

%GOR 1
disp('GOR 1')
kg=mg/((Pwf-Pr)*(Pwf^2-Pr^2))

%GOR 2
disp('GOR 2')
kg=mg/(Pwf-Pr)^2

%GOR 3
disp('GOR 3')
kg=mg/(Pwf-Pr)^3


