0
|
3157 Views
|
2 Replies
|
1 Total Likes
View groups...
Share
GROUPS:

# Physics thermoelectricity related plot

Posted 9 years ago
 Why, I cannot see the plot when I type of the following code: In:= << PhysicalConstants (* tick mark used, next to number 1 \ key*) In:= BoltzmannConstant Out= (1.38065*10^-23 Joule)/Kelvin In:= Vt[T_] := BoltzmannConstant*T/ElectronCharge In:= Vt Out= (0.0000861734 Joule)/(Coulomb Kelvin) In:= Ies = 10^-9 (* Unita A Transistor reverse current *) Ic[Vbe_, T_] := Ies*(Exp[Vbe/(Vt[T])] - 1) In:= Plot[Ic[Vbe, 273], {Vbe, 0, 0.7}] 
2 Replies
Sort By:
Posted 9 years ago
 You might edit in a little note saying what version of Mathematica you are using. The way units were handled was completely changed starting with version 9 and examples for a version you don't have won't be helpful.If you are using version 9 or later then you might compare this with what you showed. And watch for the missing _ in your Ic definition. In:= boltz = UnitConvert[Quantity["BoltzmannConstant"], "Joules"/"Kelvins"] Out= Quantity[1.38065*10^-23, ("Joules")/("Kelvins")] In:= echarge = UnitConvert[Quantity["ElectronCharge"], "Coulombs"] Out= Quantity[1.6021766*10^-19, "Coulombs"] In:= Vt[T_] := boltz*T/echarge In:= Vt Out= Quantity[0.000086173, ("Joules")/("Coulombs" "Kelvins")] In:= Ies = 10^-9 (*Unita A Transistor reverse current*) Out= 1/1000000000 In:= Ic[Vbe_, T_] := Ies*(Exp[Vbe/QuantityMagnitude[Vt[T]]] - 1) In:= Ic[Vbe, 273] Out= (-1 + E^(42.5074 Vbe))/1000000000 In:= Plot[Ic[Vbe, 273], {Vbe, 0, 0.7}] ` If you insert PlotRange->All just before the final closing ] of that Plot then you will be able to see the current out to 0.7 volts.Please check this carefully and point out any mistakes that you see. Thank you
Posted 9 years ago