function D = disparity(f,G,b)

D = ones(size(G));
%G = G.^(-1);
D = (f*b)./G;