function [answer] = obs()

format compact;

n00=96
n01=104
n10=109
n11=666
p = 0.5

answer = [];
z = [];

ntot = n00+n01+n10+n11;
omg = (n10+n11)./ntot
pi0 = n01./(n00+n01)
pi1 = n11./(n10+n11)
R = pi1./pi0

for irf = 0:10;

zrow = [];

for ith = 0:10;

rf=1+irf.*0.1;
theta=1+ith.*0.1;

denom = 1-omg + omg.*theta;
p0 = p./denom
p1 = theta.*p
fac = ((1-p1)+p1.*rf)./((1-p0)+p0.*rf)
Rstar = R./fac

answer = vertcat(answer,[rf theta Rstar]);
zrow = [zrow Rstar];

end

z = vertcat(z,zrow);

end

x = 1 + (0:10)'./10;
y = 1 + (0:10)'./10;
surf(x,y,z)