dim v(25) as double

dim v1(25) as double
dim v2(25) as double
dim v3(25) as double
dim v4(25) as double
dim v5(25) as double

dim alph as double
dim beta as double

dim weight as double

dim count as integer

dim check as integer

dim con as integer

dim n as integer

dim track as integer
dim e1 as double
dim e2 as double


weight = .01

open "weight0p01acceltest.csv" for output as #1



for n = 0 to 25
    
    v(n) = 1
    
next




track = 1

again:

for count = track to track+5
    
    for n = 1 to 24
        
        v(n) = (v(n-1) + v(n+1)+1)/2
        
    next
    
    v(25) = v(24)
    
    print #1, count;",";v(1);",";v(25)
    
    for n = 1 to 25
        
        v1(n) = v2(n)
        v2(n) = v3(n)
        v3(n) = v4(n)
        v4(n) = v5(n)
        v5(n) = v(n)
    next
    
next

track = track + 6

count = track

con = 0

do
    
    for n = 1 to 24
        v(n) = (v(n-1) + v(n+1)+1)/2
    next
    
    v(25) = v(24)
    
    for n = 1 to 25
        
        v1(n) = v2(n)
        v2(n) = v3(n)
        v3(n) = v4(n)
        v4(n) = v5(n)
        v5(n) = v(n)
    next
    
    check = 1
    
    
    for n = 1 to 25
        
        alph = sqr((v5(n)-v4(n))/(v3(n)-v2(n)))

        if alph >= 1 then check = 0
        
    next
    
    if check = 1 then con = con + 1
    

    print #1, count;",";v(1);",";v(25)
    
    count = count + 1
    track = track + 1

loop until con >= 25 or count > 2000

if count > 2000 then goto finish 
    
    
REM Extrapolation step

for n = 1 to 25
    
        
    alph = sqr((v5(n)-v4(n))/(v3(n)-v2(n)))
    
    beta = (v5(n)+v4(n)+v3(n))/3 - (v4(n)+v3(n)+v2(n))*alph/3
    
    e1 = beta/(1-alph)
    
        
    alph = sqr((v4(n)-v3(n))/(v2(n)-v1(n)))
    
    beta = (v4(n)+v3(n)+v2(n))/3 - (v3(n)+v2(n)+v1(n))*alph/3
    
    e2 = beta/(1-alph)
    
    v(n) = (2-weight)/(weight/v(n) + (1-weight)*(1/e1+1/e2))
    
next


    print #1, count;",";v(1);",";v(25)
    
    count = count + 1
    track = track + 1
    
    goto again
    
    finish:
    
    close #1