The Code for Single Channel model

real u(0:10000)
iseed = 11208244
q1 = 1
q2 = 2
N = 10
Do j = 1, N
         t2 = 0.0
         Do i = 1, 15
                 t1 = t2
                 t2 = t2-log(1-ran(iseed))/q1
                 x = 1
                 Do k = int(t1*100), int(t2*100)-1
                             u(k) = u(k)+x
                 End Do
                 t1 = t2
                 t2 = t2-log(1-ran(iseed))/q1
                 x = 0
                 Do k = int(t1*100), int(t2*100)-1
                             u(k) = u(k)+x
                 End Do
        End do
End do

Do k = 0, 600
              write (10,*) k*0.01, u(k)/N
End Do
End