base=getresult('CHARGE','base');
Ib = base.I/1e4;  # Amp/um
Vb = base.V_base;
Vc = base.V_collector;

collector=getresult('CHARGE','collector');
Ic = collector.I/1e4;  # Amp/um

plot(Ic,Ic/Ib,'Ic [A/um]','Beta','','log10x');
legend('');

plot(Vb,Ib,Ic,'Vbase [V]','Current [A/um]','','log10y');
legend('Ib','Ic');