function lambda = compute_lambda_msg_nor(bnode, bnet, p)

% Pearl p190 top eqn

n = bnode.self;
l0 = bnet.pearl.msg{n}.lambda(1);
l1 = bnet.pearl.msg{n}.lambda(2);

pi = prod_CPT_and_pi_msgs_nor(bnet, n, p);
j = bnet.pearl.parent_index{n}(p); % p is my j'th parent
q = bnode.inhibit(j+1); % skip the leak
lambda = zeros(2,1);
for u=0:1
  lambda(u+1) = l1 - (q^u)*(l1 - l0)*pi;
end
