function s = roehre(s, gain, k1)
N = size(s);

k1 = 0;

for n=1:N(1)
    for m=1:N(2)

        if (gain == 0)
         gain = 000000000000001;
        end;

        if (k1 == 0)
        k1 = 0.000000000000001;
        end;


  if (s(n,m) ~= k1)

       s(n,m) = ((s(n,m)-k1)/(1-2.7182818284^(-gain*s(n,m)-k1))) + (k1/(1-2.7182818284^(gain*k1)));
  else

      s(n,m) = (1/gain) + (k1/(1-2.7182818284^(gain*k1)));
  end



        %2.7182818284


    end;
end;