* NPN with BF=100: force 10uA base current, recover beta = Ic/Ib.
* Ic is slightly > BF because of the Early effect (VAF=80) at Vce=5.
.model MYNPN NPN(IS=1e-16 BF=100 VAF=80)
IB 0 b DC 10u
Q1 c b 0 MYNPN
VC c 0 DC 5
.op
.control
  run
  let beta = -i(vc)/10u      $ negate: source current is into VC
  print beta
.endc
.end