function [a,b,c,d,e] = fill_x(JM,KM,A,dx,dt,jmax,kmax,epse)

a = zeros(jmax,kmax,3,3);
b = a;
c = a;
d = a;
e = a;

hx = 0.5*dt/dx;
he = dt/dx*epse;
for n = 1:3
  for m = 1:3
    b(JM,KM,n,m) = -A(n,m)*hx;
    d(JM,KM,n,m) = A(n,m)*hx;
    c(JM,KM,n,m) = 0.0;
  end
  b(JM,KM,n,n) = b(JM,KM,n,n) - he;
  d(JM,KM,n,n) = d(JM,KM,n,n) - he;
  c(JM,KM,n,n) = 1.0 + c(JM,KM,n,n) + 2*he;
end