Monte Carlo Random Walk

%Random yürüme işlemi buradan başlıyor
        lmax = imax;
        gmax = jmax;
        p = zeros(lmax,gmax); %rastgele yürüme matrisi

        l = i;
        g = j; %l satır nodu, g sütun nodu
        S = 0;
        while g >= 1 || g <= gmax; %g=1 ve g=gmax sınır şartları bu değerlere
            %ulaştığı zaman döngü durup ilgili nodun sıcaklığı kaydedilecek
            r = rand(1);
            if r<= 0.25
                g = g+1;
            elseif r<=0.5
               if l == 1
                   l = l+0;
               else
                   l = l+1;
               end
            elseif r<=0.75
                g = g-1;
            else
                if l == lmax
                    l = l-0;
                else
                    l = l-1;
                end
            end

            if g==1 %Sınır yüzeyine ulaşıldı, sıcaklığı kaydet.
                S = S+Ts1;
            elseif g== gmax
                S = S+Ts2;
            end