//example 3 pagenumber 139
clear
v=4.3;//volt
q=4;//volt
dop=10^17;//per cubic centimetre
fi0=0.254*log(dop/(5.1*10^10));
fi01=0.407+q+0.55;
disp('fi0   =   '+string(fi01));