Biyokimyasal reaksiyon sistemleri farkli reaksiyonlar araciligiyla etkilesime giren bircok farkli turu icerir. Sistem icerisinde yer alan turlerin sayilari ve miktarlari cok yuksek oldugunda, diferansiyel denklemlere dayanan saf modelleme yaklasimlari cok boyutluluktan muzdarip olurlar. Eger bir sistem korunumlu donguler icerirse, bazi turlerin miktarlari cebirsel baglantilar yoluyla elde edilebilir bu da sistemin dinamiklerini temsil eden diferansiyel denklemlerin boyutunu dusurur. Bu calismada, biyokimyasal reaksiyon sistemlerinde yer alan korunumlu donguleri elde etmek icin Gauss-Jordan metodunu kullanan bir numerik algoritma oneriyoruz. Algoritmayi stokastik modelleme yaklasiminda konum vektorunun tam realizasyonlarini elde eden Direk Metod (DM), Ilk Reaksiyon Metodu (FRM) ve Sonraki Reaksiyon Metodu (FRM) icerisinde verdik. Bu uc algoritmayi korunum bagintilarini icerecek/icermecek sekilde farkli boyutlardaki biyokimyasal sistemlere uyguladik ve her tam lagoritmanin farkli iki versiyonunun hesaplama miktarlari kiyasladik.