function Rp = parallel_R(R1,R2)
% Compute parallel R value

if R1 < 0 | R2 < 0
    display('oops!  negative R...')
    return
end

Rp = 1 ./ (1/R1 + 1/R2);
return

end

