It is very important to analyze the polluted water that streams transfer around them by diffusion. Dirty water leaking 
from streams like this causes soil and water pollution. The protection of surface and underground water resources is 
very important given the decreasing water resources. How deep a pollutant has infiltrated can be found by solving 
the diffusion equation. In this study, the concentration amount of the pollutant passing from a long-polluted stream 
to the soil was investigated by the finite difference method. Instead of using expensive packet programs, a widespread 
program Matlab used by engineers and scientists to analyze data, develop algorithms, solve problems analytically 
and numerically, and create models is preferred in this study. The cross-section of the stream is rectangular and its 
length is considered to be quite long. Therefore, the problem has been examined in two dimensions. First, the 
necessary finite difference equations are derived and this problem is analyzed using appropriate boundary conditions. 
A numerical solution has been obtained for two cases in which the bottom of the stream is permeable and it is not. 
The mass flow rate of both of the cases was calculated and it is found that the mass flow rate in the first case in 
which the bottom of the stream is permeable is almost double of the second case. The flow velocities on the 
boundaries have been found to be less in the case in which the bottom of the stream is impermeable.
 
Akarsuların etraflarına difüzyon yolu ile aktardıkları kirli suyun analizi oldukça önemlidir. Bu akarsulardan sızan 
kirli su toprak ve su kirliliğine neden olmaktadır. Yüzey ve yeraltı su kaynaklarının korunması giderek azalan su 
kaynakları dikkate alındığında çok önemlidir.  Bir kirleticinin ne kadar derine sızdığı difüzyon denklemi çözülerek 
bulunabilir. Bu çalışmada atıksularla kirlenmiş uzun bir dereden toprağa geçen kirleticinin konsantrasyon miktarı 
sonlu farklar yöntemi ile incelenmiştir. Bu çalışmada pahalı paket programlar kullanmak yerine mühendisler ve bilim 
insanlarının verileri analiz etmek, algoritma geliştirmek, problemleri analitik ve sayısal olarak çözmek ve model 
oluşturmak için kullandıkları yaygın bir program olan Matlab tercih edilmiştir. Derenin kesiti diktörtgen alınmış ve 
boyu oldukça uzun kabul edilmiştir. Bundan dolayı problem iki boyutta incelenmiştir. Önce gerekli sonlu farklar 
denklemleri türetilmiş ve uygun sınır şartları kullanılarak problem analiz edilmiştir. Akarsuyun tabanının geçirimli 
olduğu ve olmadığı iki durum için sayısal çözüm elde edilmiştir. Her iki durumun kütlesel debisi hesaplanmış ve 
akarsuyun tabanının geçirimli olduğu birinci durumdaki kütlesel debinin, ikinci durumun neredeyse iki katı olduğu 
bulunmuştur. Akarsu tabanının geçirimsiz olması durumunda sınırlardaki akış hızlarının daha az olduğu tespit 
edilmiştir.