I am using the impedande class to compute the electrotonic distance of some dendrites. To achieve this I used a slightly modified code from https://www.neuron.yale.edu/phpBB2/viewtopic.php?t=771. The code I am using is the following:
Code: Select all
// prepare to use Impedance class
// always a good idea to finitialize before computing impedance
v_init=-70
finitialize(v_init)
// demonstrate use of impedance class
objref zz
zz = new Impedance()
FREQ = 0 // Hz
WHERE = 0.5 // location in the soma that is the reference point
soma distance(0, WHERE) // sets origin for distance calculations
proc calcZ() {
soma zz.loc(WHERE) // sets origin for impedance calculations
zz.compute(FREQ, 1) // takes the impedance contributions of
// gating state differential equations into account
// but requires mechanisms to be compatible with CVODE
print "distance(x) ratio(x)"
forall {
print secname()
for (x) if (x==1) print x3d(x), y3d(x), z3d(x), log(1/(zz.ratio(x)))
}
}
calcZ()
Any ideas on how I can print more decimal places??
Thank you all in advance