>> 
>> test
>> o
>> v
>> test
>> test
>> plot(v,a')
>> clf
>> whos
>> test
>> close all
>> whos
>> x = v-tha
>> whos
>> plot(x)
>> a = Ra  * x / (1 - exp(-x/qa));
>> whos
>> a = Ra  * x 
>> a = Ra  * (x ./ (1 - exp(-x/qa)));
>> whos
>> plot(v,a)
>> clf
>> tes
>> test
>> test
>> test
>> test
>> test
>> close all
>> test
>> test
>> close all
>> close all
>> test
>> test
>> isnan(a)
>> isnan(b)
>> 
>> x = -120:100;
>> whos
>> x = -120:0.9:100;
>> whos
>> x = -120:01.1:100;
>> whos
>> 290./199
>> 220./199
>> x = -120:01.1055:100;
>> whos
>> x = -120:01.1056:100;
>> whos
>> find(x<20)
>> x(126:129)
>> tan( 0.284082)
>> help deg
>> tan(2*pi)
>> tan(pi/4)
>> tan(2.393618)
>> tan(1.539648)
>> floor(-1.5)
>> plot(x,tan(x)
>> plot(x,tan(x))
>> close all
>> plot(x,tan(x))
>> x = -10*pi:0.01:10*pi;
>> plot(x,tan(x))
>> x = -10*pi:0.00001:10*pi;
>> x = -10*pi:0.0001:10*pi;
>> plot(x,tan(x))
>> x = -2*pi:0.001:2*pi;
>> plot(x,tan(x))
>> x = -pi:0.001:pi;
>> plot(x,tan(x))
>> x = -2:0.001:2;
>> plot(x,tan(x))
>> x = -pi:0.0001:pi;
>> plot(x,x)
>> plot(x,tan(x))
>> semilogy(x,tan(x))
>> 
>> x = -pi/2:0.0001:pi/2;
>> semilogy(x,tan(x))
>> 
>> 
>> 
>> pi/2
>> pi/2 - 1.567510
>> pi/2 - 1.561663
>> x = 0:0.0001:pi;
>> semilogy(2*x/pi,tan(x))
>> 
>> 
>> 
>> plot(x,tan(x))
>> semilogy(x,tan(x))
>> 
>> axis([0 pi -100 100])
>> axis([1.5 1.6 -100 100])
>> pi/2
>> pi//2 - 1.56
>> pi/2 - 1.56
>> clear all
>> close 
>> help random
>> x = random('unif',1,5000);
>> y = random('unif',1,5000);
>> z = random('unif',1,5000);
>> plot(z*pi,tan(z*pi))
>> close all
>> plot(z*pi,tan(z*pi))
>> plot(z*pi,tan(z*pi),'o')
>> whos
>> help random
>> help rand
>> x = rand(1,5000);
>> y = rand(1,5000);
>> z = rand(1,5000);
>> plot(z*pi,tan(z*pi))
>> figure
>> sum(find(x.^2 + y.^2 < 1.0) 
>> sum(find(x.^2 + y.^2 < 1.0) )
>> sum(find(x.^2 + y.^2 > 1.0) )
>> max(x)
>> min(x)
>> x = 2.0*x = 1;
>> x = 2.0*x - 1;
>> min(x)
>> max(x)
>> y = 2.0*y - 1;
>> sum(find(x.^2 + y.^2 > 1.0) )
>> sum(find(x.^2 + y.^2 1 1.0) )
>> sum(find(x.^2 + y.^2 < 1.0) )
>> (x.^2 + y.^2 < 1.0)
>> (x.^2 + y.^2 > 1.0)
>> plot(1:5000,(x.^2 + y.^2 < 1.0))
>> plot(1:5000,(x.^2 + y.^2 < 1.0),'o')
>> plot(1:5000,(x.^2 + y.^2 < 1.0),'.')
>> plot(1:5000,(x.^2 + y.^2 < 1.0),'o')
>> plot(1:5000,(x.^2 + y.^2 > 1.0),'o')
>> clf
>> sum(find(x.^2 + y.^2 < 1.0) )
>> j = (x.^2 + y.^2 > 1.0);
>> whos
>> sum(j)
>> j = (x.^2 + y.^2 < 1.0);
>> sum(j)
>> j = (x.^2 + y.^2 > 1.0);
>> x(j) = [];
>> y(j) = [];
>> whos
>> x./y;
>> plot(1:length(x),x./y)
>> semilog(1:length(x),x./y)
>> semilogy(1:length(x),x./y)
>> 
>> close
>> quit
>> testatn
>> what
>> testtan
>> z
>> plot(z,tan(z))
>> close 
>> 
>> clear all
>> testtan
>> plot(z,tan(z))
>> plot(z,tan(z),'o')
>> hold on
>> plot(tz,t,'o')
>> plot(tz,t,'ro')
>> close
>> semilogy(tz,t,'ro')
>> 
>> semilogy(z,tan(z),'bo')
>> 
>> semilogy(z,tan(z),'bo')
>> 
>> semilogy(tz,t,'ro')
>> 
>> semilogy(tz,tan(tz),'go')
>> 
>> semilogy(x,t,'ro')
>> 
>> 
>> close all
>> 2/0.05
>> 2/0.01
>> 1/0.01
>> 1/0.00011
>> 
>> tan(pi/2)
>> help gamma
>> gammaln(tan(pi/2))
>> 
>> whos
>> save
>> clear all
>> load a_kd3
>> losd a_sk
>> load a_sk
>> plotxy(a_sk)
>> plotxy(a_kd3)
>> plotxy(a_kd3,'r')
>> hold on
>> load b_sk
>> load b_kd3
>> plotxy(b_kd3)
>> plotxy(b_kd3,'r')
>> plotxy(b_sk)
>> plotxy(b_sk.*2)
>> clf
>> plot(b_sk(:,1),b_sk(:,2)*2)
>> hold on
>> plotxy(b_kd3,'r')
>> figure
>> plot(a_sk(:,1),a_sk(:,2)*2)
>> hold on
>> plotxy(a_kd3,'r')
>> figure
>> plotxy(a_kd3(:,1),a_kd3(:,2)-a_sk(:,2)*2,'r')
>> plot(a_kd3(:,1),a_kd3(:,2)-a_sk(:,2)*2,'r')
>> figure
>> plot(b_kd3(:,1),b_kd3(:,2)-b_sk(:,2)*2,'r')
>> 1/0.025
>> v = 40
>> rate(v)
>> rates(v)
>> rates_kd3(v)
>> close all
>> load b_sk1
>> plotxy(b_sk)
>> plotxy(b_sk1,'r')
>> close
>> figure
>> plotxy(a_sk)
>> axis([-0.1 60 -0.1 0.01])
>> close all
>> plot(b_kd3(:,1),b_kd3(:,2)-b_sk(:,2)*2.1576693,'r')
>> tan(0.0566722)
>> arctan(-7.61955e+231)
>> atan(-7.61955e+231)
>> pi/2
>> tan(1.58571)
>> 
>> whos
>> quit
>> load test_ran1.out
>> plot(test_ran1)
>> plot(test_ran1,'.')
>> 
>> psd(test_ran1)
>> spectrum(test_ran1)
>> close 
>> quit
>> load somas.dat
>> load soma.dat
>> plotxy(somas)
>> plotxy(soma,'r')
>> plotxy(hilloc,'r')
>> whos
>> load hilloc
>> load hilloc.dat
>> load hillocs.dat
>> plotxy(hilloc,'r')
>> plotxy(hillocs,'m')
>> load dend.dat
>> load dends.dat
>> plotxy(dends,'m')
>> plotxy(dend,'r')
>> close
>> figure
>> subplot(3,1,1)
>> load dends.dat
>> plotxy(dend,'r')
>> plotxy(dends,'m')
>> subplot(3,1,2)
>> plotxy(soma,'r')
>> plotxy(somas,'m')
>> subplot(3,1,3)
>> plotxy(hillocs,'m')
>> plotxy(hilloc,'r')
>> axis([20 40 -80 50]
>> axis([20 40 -80 50])
>> axis([20 40 -80 70])
>> 
>> quit
>> 
>> test
>> o
>> v
>> test
>> test
>> plot(v,a')
>> clf
>> whos
>> test
>> close all
>> whos
>> x = v-tha
>> whos
>> plot(x)
>> a = Ra  * x / (1 - exp(-x/qa));
>> whos
>> a = Ra  * x 
>> a = Ra  * (x ./ (1 - exp(-x/qa)));
>> whos
>> plot(v,a)
>> clf
>> tes
>> test
>> test
>> test
>> test
>> test
>> close all
>> test
>> test
>> close all
>> close all
>> test
>> test
>> isnan(a)
>> isnan(b)
>> 
>> x = -120:100;
>> whos
>> x = -120:0.9:100;
>> whos
>> x = -120:01.1:100;
>> whos
>> 290./199
>> 220./199
>> x = -120:01.1055:100;
>> whos
>> x = -120:01.1056:100;
>> whos
>> find(x<20)
>> x(126:129)
>> tan( 0.284082)
>> help deg
>> tan(2*pi)
>> tan(pi/4)
>> tan(2.393618)
>> tan(1.539648)
>> floor(-1.5)
>> plot(x,tan(x)
>> plot(x,tan(x))
>> close all
>> plot(x,tan(x))
>> x = -10*pi:0.01:10*pi;
>> plot(x,tan(x))
>> x = -10*pi:0.00001:10*pi;
>> x = -10*pi:0.0001:10*pi;
>> plot(x,tan(x))
>> x = -2*pi:0.001:2*pi;
>> plot(x,tan(x))
>> x = -pi:0.001:pi;
>> plot(x,tan(x))
>> x = -2:0.001:2;
>> plot(x,tan(x))
>> x = -pi:0.0001:pi;
>> plot(x,x)
>> plot(x,tan(x))
>> semilogy(x,tan(x))
>> 
>> x = -pi/2:0.0001:pi/2;
>> semilogy(x,tan(x))
>> 
>> 
>> 
>> pi/2
>> pi/2 - 1.567510
>> pi/2 - 1.561663
>> x = 0:0.0001:pi;
>> semilogy(2*x/pi,tan(x))
>> 
>> 
>> 
>> plot(x,tan(x))
>> semilogy(x,tan(x))
>> 
>> axis([0 pi -100 100])
>> axis([1.5 1.6 -100 100])
>> pi/2
>> pi//2 - 1.56
>> pi/2 - 1.56
>> clear all
>> close 
>> help random
>> x = random('unif',1,5000);
>> y = random('unif',1,5000);
>> z = random('unif',1,5000);
>> plot(z*pi,tan(z*pi))
>> close all
>> plot(z*pi,tan(z*pi))
>> plot(z*pi,tan(z*pi),'o')
>> whos
>> help random
>> help rand
>> x = rand(1,5000);
>> y = rand(1,5000);
>> z = rand(1,5000);
>> plot(z*pi,tan(z*pi))
>> figure
>> sum(find(x.^2 + y.^2 < 1.0) 
>> sum(find(x.^2 + y.^2 < 1.0) )
>> sum(find(x.^2 + y.^2 > 1.0) )
>> max(x)
>> min(x)
>> x = 2.0*x = 1;
>> x = 2.0*x - 1;
>> min(x)
>> max(x)
>> y = 2.0*y - 1;
>> sum(find(x.^2 + y.^2 > 1.0) )
>> sum(find(x.^2 + y.^2 1 1.0) )
>> sum(find(x.^2 + y.^2 < 1.0) )
>> (x.^2 + y.^2 < 1.0)
>> (x.^2 + y.^2 > 1.0)
>> plot(1:5000,(x.^2 + y.^2 < 1.0))
>> plot(1:5000,(x.^2 + y.^2 < 1.0),'o')
>> plot(1:5000,(x.^2 + y.^2 < 1.0),'.')
>> plot(1:5000,(x.^2 + y.^2 < 1.0),'o')
>> plot(1:5000,(x.^2 + y.^2 > 1.0),'o')
>> clf
>> sum(find(x.^2 + y.^2 < 1.0) )
>> j = (x.^2 + y.^2 > 1.0);
>> whos
>> sum(j)
>> j = (x.^2 + y.^2 < 1.0);
>> sum(j)
>> j = (x.^2 + y.^2 > 1.0);
>> x(j) = [];
>> y(j) = [];
>> whos
>> x./y;
>> plot(1:length(x),x./y)
>> semilog(1:length(x),x./y)
>> semilogy(1:length(x),x./y)
>> 
>> close
>> quit