load_file("nrngui.hoc")
taudt=1
func C() { // n
return (1 - exp(-$1*taudt))/(1 - exp(-taudt))
}
func T() {
return C($1)^4
}
func F() { return T($1) - 1 }
i=0
proc table1() {
printf("DT/tau\tF1\tF2\tF3\tF1000\n")
for case(&taudt, 1e-9, .5, .73, .75, 1, 1e9) {
printf("%-5.2g ", taudt)
for case(&i, 1, 2, 3, 1000) {
printf(" %-7.1f", F(i))
}
printf("\n")
}
}
objref Avec, Bvec, Nvec, g
{
Nvec = new Vector(11)
Nvec.indgen(1,1)
Avec = Nvec.c
Bvec = Nvec.c
{taudt = 0.83 Avec.apply("F")}
{taudt = 0.50 Bvec.apply("F")}
load_file("start.ses")
g = Graph[0]
Avec.mark(g, Nvec.c.sub(1), "+", 10, 1)
Bvec.mark(g, Nvec.c.sub(1), "+", 10, 2)
}
table1()