%Response of photoreceptor
%Wilson page 35ff
%   kja 070504
aminit;
t=0:0.001:0.4;
k=4;tc=0.025;th=0.08;
A=[-1/tc -k/tc;1/th -1/th];
B=[1/tc;0];C=[1 0];D=0;
sys=ss(A,B,C,D);
y=step(sys,t);
subplot(331);
pl=plot(t,y,'b-');
set(pl,'Linewidth', AM_data_linewidth);
amaxis([0 0.4 0 0.601]);

xlabel('t');ylabel('y');

amprint('neural-step.eps');
