function A = antimcrobial_solve_release(H,delta_z,D_a,ew,delta_ax,beta_a,Q,ex,Ab,phi,j,delta_t,t_bar,delta_a,delta_ab,eb)

    A = zeros(H/delta_z+1,H/delta_z+1);
    B = delta_z^2*beta_a/(ew*D_a)*phi(1:(H/delta_z+1))*heaviside(t_bar-j*delta_t);
    B(end) = B(end)+ 2*delta_z*Q/(ew*D_a)*Ab(end); 
    for i = 1:(H/delta_z+1)
        if i == 1
            A(i,i) = delta_z^2/(ew*D_a)*(delta_ax*ex(i)+delta_a*ew+delta_ab*eb) + 2;
            A(i,i+1) = -2;
        elseif i == H/delta_z+1
            A(i,i) = delta_z^2/(ew*D_a)*(delta_ax*ex(i)+delta_a*ew+delta_ab*eb) + 2 + 2*delta_z*Q/(ew*D_a); 
            A(i,i-1) = -2;
        else
            A(i,i) = delta_z^2/(ew*D_a)*(delta_ax*ex(i)+delta_a*ew+delta_ab*eb) + 2;
            A(i,i-1) = -1;
            A(i,i+1) = -1;
        end
    end
    A = linsolve(A,B);

end