L = 10; % length of the domain
d = .0005; % diffusion coefficient
h=L/N;
dt = 2*h;
    alf = 0.05;
    epsilon = 0.01;
    ga = 8.0;
  gs = 8.0;
  mu1 = 0.05;
  mu2 = 0.3;
  eps1 = 0.03;
  eps2 = 0.0001;
  Istim = 0.02;
 
  
