% rc1exact
t = 0:.1:10;
i0 = t.*exp(-t);
v = t.*i0/2;
[ax,h1,h2] = plotyy(t,i0,t,v);
xlabel('t (ms)','fontsize',14)
ylabel(ax(1),'i_0 (\mu A)','fontsize',14)
ylabel(ax(2),'x (mV)','fontsize',14)
set(h1,'linewidth',2)
set(h2,'linewidth',2)
