## Set file for koho3.ode on Mon Apr 18 11:04:19 2016
304   Number of equations and auxiliaries
6   Number of parameters
# Numerical stuff
8    nout
40    nullcline mesh
0   Discrete
80000  total
1  DeltaT
0  T0
40000  Transient
100  Bound
1e-12  DtMin
1  DtMax
0.001  Tolerance
0.001  Abs. Tolerance
0  Max Delay
100   Eigenvector iterates
0.001  Eigenvector tolerance
0.001  Newton tolerance
0  Poincare plane
1e-05  Boundary value tolerance
1e-05  Boundary value epsilon
20   Boundary value iterates
0   Poincare None
1   Poincare variable
1   Poincare sign
0   Stop on Section
0   Delay flag
80000  Current time
80000  Last Time
0   MyStart
1   INFLAG
# Delays
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
# Bndry conds
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Old ICs
0  WA0
0  WA1
0  WA2
0  WA3
0  WA4
0  WA5
0  WA6
0  WA7
0  WA8
0  WA9
0  WA10
0  WA11
0  WA12
0  WA13
0  WA14
0  WA15
0  WA16
0  WA17
0  WA18
0  WA19
0  WA20
0  WA21
0  WA22
0  WA23
0  WA24
0  WA25
0  WA26
0  WA27
0  WA28
0  WA29
0  WA30
0  WA31
0  WA32
0  WA33
0  WA34
0  WA35
0  WA36
0  WA37
0  WA38
0  WA39
0  WA40
0  WA41
0  WA42
0  WA43
0  WA44
0  WA45
0  WA46
0  WA47
0  WA48
0  WA49
0  WA50
0  WA51
0  WA52
0  WA53
0  WA54
0  WA55
0  WA56
0  WA57
0  WA58
0  WA59
0  WA60
0  WA61
0  WA62
0  WA63
0  WA64
0  WA65
0  WA66
0  WA67
0  WA68
0  WA69
0  WA70
0  WA71
0  WA72
0  WA73
0  WA74
0  WA75
0  WA76
0  WA77
0  WA78
0  WA79
0  WA80
0  WA81
0  WA82
0  WA83
0  WA84
0  WA85
0  WA86
0  WA87
0  WA88
0  WA89
0  WA90
0  WA91
0  WA92
0  WA93
0  WA94
0  WA95
0  WA96
0  WA97
0  WA98
0  WA99
0  WB0
0  WB1
0  WB2
0  WB3
0  WB4
0  WB5
0  WB6
0  WB7
0  WB8
0  WB9
0  WB10
0  WB11
0  WB12
0  WB13
0  WB14
0  WB15
0  WB16
0  WB17
0  WB18
0  WB19
0  WB20
0  WB21
0  WB22
0  WB23
0  WB24
0  WB25
0  WB26
0  WB27
0  WB28
0  WB29
0  WB30
0  WB31
0  WB32
0  WB33
0  WB34
0  WB35
0  WB36
0  WB37
0  WB38
0  WB39
0  WB40
0  WB41
0  WB42
0  WB43
0  WB44
0  WB45
0  WB46
0  WB47
0  WB48
0  WB49
0  WB50
0  WB51
0  WB52
0  WB53
0  WB54
0  WB55
0  WB56
0  WB57
0  WB58
0  WB59
0  WB60
0  WB61
0  WB62
0  WB63
0  WB64
0  WB65
0  WB66
0  WB67
0  WB68
0  WB69
0  WB70
0  WB71
0  WB72
0  WB73
0  WB74
0  WB75
0  WB76
0  WB77
0  WB78
0  WB79
0  WB80
0  WB81
0  WB82
0  WB83
0  WB84
0  WB85
0  WB86
0  WB87
0  WB88
0  WB89
0  WB90
0  WB91
0  WB92
0  WB93
0  WB94
0  WB95
0  WB96
0  WB97
0  WB98
0  WB99
0  WC0
0  WC1
0  WC2
0  WC3
0  WC4
0  WC5
0  WC6
0  WC7
0  WC8
0  WC9
0  WC10
0  WC11
0  WC12
0  WC13
0  WC14
0  WC15
0  WC16
0  WC17
0  WC18
0  WC19
0  WC20
0  WC21
0  WC22
0  WC23
0  WC24
0  WC25
0  WC26
0  WC27
0  WC28
0  WC29
0  WC30
0  WC31
0  WC32
0  WC33
0  WC34
0  WC35
0  WC36
0  WC37
0  WC38
0  WC39
0  WC40
0  WC41
0  WC42
0  WC43
0  WC44
0  WC45
0  WC46
0  WC47
0  WC48
0  WC49
0  WC50
0  WC51
0  WC52
0  WC53
0  WC54
0  WC55
0  WC56
0  WC57
0  WC58
0  WC59
0  WC60
0  WC61
0  WC62
0  WC63
0  WC64
0  WC65
0  WC66
0  WC67
0  WC68
0  WC69
0  WC70
0  WC71
0  WC72
0  WC73
0  WC74
0  WC75
0  WC76
0  WC77
0  WC78
0  WC79
0  WC80
0  WC81
0  WC82
0  WC83
0  WC84
0  WC85
0  WC86
0  WC87
0  WC88
0  WC89
0  WC90
0  WC91
0  WC92
0  WC93
0  WC94
0  WC95
0  WC96
0  WC97
0  WC98
0  WC99
0  IP
# Ending  ICs
-0.4325703050856919  WA0
-0.4256948931440818  WA1
-0.421042994792578  WA2
-0.4159843443777203  WA3
-0.4060355495681908  WA4
-0.3857712263916847  WA5
-0.3504946015250518  WA6
-0.2979820336721327  WA7
-0.2293997031268988  WA8
-0.1490684971433481  WA9
-0.06352210889450818  WA10
0.01959824085920886  WA11
0.0927456656612041  WA12
0.1498902680381344  WA13
0.1878786010093164  WA14
0.2070803509990505  WA15
0.2109412427257937  WA16
0.2046825839901381  WA17
0.1937437648619318  WA18
0.1825053924956924  WA19
0.1735902502059724  WA20
0.1678062404441446  WA21
0.1645910048580769  WA22
0.1626416681114393  WA23
0.1604081853602362  WA24
0.1563865123358104  WA25
0.149402075670783  WA26
0.1389725672722915  WA27
0.1255415452315191  WA28
0.1103940900129286  WA29
0.09540476998825283  WA30
0.08289032555329137  WA31
0.07558213704756564  WA32
0.0764938352592668  WA33
0.08850685848900271  WA34
0.1136822503928575  WA35
0.1524409320133773  WA36
0.2028758560613259  WA37
0.2605864996903977  WA38
0.3193521384997532  WA39
0.3725720866638056  WA40
0.4149585896559912  WA41
0.443798523600883  WA42
0.4592746161810321  WA43
0.463772661783527  WA44
0.4606010781954944  WA45
0.4527760270263666  WA46
0.4423123156376656  WA47
0.4300463357874271  WA48
0.4157773043183713  WA49
0.3985532666129314  WA50
0.3770526774335231  WA51
0.3500364983311299  WA52
0.3167904542313013  WA53
0.2774329617258691  WA54
0.2329692448274515  WA55
0.1850503382109572  WA56
0.1355389089245782  WA57
0.08610263353727139  WA58
0.03804739029074725  WA59
-0.007541611616897054  WA60
-0.04946824981994473  WA61
-0.08613075962883986  WA62
-0.1154568569253136  WA63
-0.1351421417779385  WA64
-0.1432081695144411  WA65
-0.138836560725872  WA66
-0.1231087259231449  WA67
-0.09894017727754501  WA68
-0.0699696039151028  WA69
-0.0392889132028407  WA70
-0.009123576443982852  WA71
0.01858257696899028  WA72
0.04170168187047558  WA73
0.05831729012482856  WA74
0.0677999018246369  WA75
0.0717535179402041  WA76
0.07401300272723263  WA77
0.07939092604433819  WA78
0.09161118578404161  WA79
0.1113677343179235  WA80
0.1354380220547149  WA81
0.1573357540697696  WA82
0.1692823036432784  WA83
0.1645838693912441  WA84
0.1393200761848496  WA85
0.09282790801553624  WA86
0.02731388955720741  WA87
-0.05272550877414443  WA88
-0.1409484672463054  WA89
-0.2297037192007949  WA90
-0.3110516399750881  WA91
-0.3782827364203012  WA92
-0.427369265084549  WA93
-0.4576160027352628  WA94
-0.4712335107254438  WA95
-0.4721976514545731  WA96
-0.4650300007818954  WA97
-0.4539332670509846  WA98
-0.4423257724610709  WA99
1.916920368917645  WB0
1.834415626427921  WB1
1.743138412675568  WB2
1.643353544064046  WB3
1.536094843731491  WB4
1.423189694540903  WB5
1.307300911062856  WB6
1.191770471365263  WB7
1.080089357923037  WB8
0.9752712445347156  WB9
0.8796936360838278  WB10
0.7955248108028462  WB11
0.725145909336  WB12
0.670863087004639  WB13
0.6338413877080077  WB14
0.6130018135615186  WB15
0.6048583848541356  WB16
0.6046117402555574  WB17
0.6078000858862695  WB18
0.6114360016301522  WB19
0.6141258648687503  WB20
0.6154651061409351  WB21
0.6153079467895485  WB22
0.6133373691357761  WB23
0.6090957779989367  WB24
0.602418557053774  WB25
0.593983623170802  WB26
0.5855992065484017  WB27
0.5800958800341101  WB28
0.581044472292268  WB29
0.5925600221962382  WB30
0.6191473624794969  WB31
0.665278355869895  WB32
0.7344821126207247  WB33
0.8281383968272806  WB34
0.9445228944086758  WB35
1.078631981896248  WB36
1.222997966810883  WB37
1.369345270328549  WB38
1.510545875063206  WB39
1.641982280080267  WB40
1.761583734057457  WB41
1.86866066345247  WB42
1.962514933896979  WB43
2.041815396609165  WB44
2.104979545059226  WB45
2.151052927171097  WB46
2.180398844064731  WB47
2.194893993323996  WB48
2.197804715499186  WB49
2.193611934125317  WB50
2.187787297664943  WB51
2.186347310898558  WB52
2.195158150229172  WB53
2.219207678891412  WB54
2.262110725565767  WB55
2.325942193457219  WB56
2.411249806821285  WB57
2.516951687055379  WB58
2.639944424546242  WB59
2.774664221567475  WB60
2.913229283004145  WB61
3.046629593229202  WB62
3.166627386598277  WB63
3.267330748381773  WB64
3.345696917756522  WB65
3.401261552248983  WB66
3.435882812943535  WB67
3.453615647732717  WB68
3.460050798358396  WB69
3.460868305067264  WB70
3.46039476389787  WB71
3.461046721863104  WB72
3.463670571727238  WB73
3.468179044776201  WB74
3.47400792454326  WB75
3.480283697170288  WB76
3.48577023656304  WB77
3.488681036701219  WB78
3.486489891485872  WB79
3.47595340060426  WB80
3.453523434829783  WB81
3.416122614245674  WB82
3.362013654509967  WB83
3.291388122593267  WB84
3.206386088018096  WB85
3.110491135915718  WB86
3.007563383002203  WB87
2.901003616338281  WB88
2.793461208478909  WB89
2.687125225792804  WB90
2.584227980676407  WB91
2.487226473703089  WB92
2.398366139128197  WB93
2.318826277816663  WB94
2.248010557344783  WB95
2.183488056499379  WB96
2.121644459089899  WB97
2.058641307168813  WB98
1.991190657145399  WB99
2.608903088537248  WC0
2.359781499820363  WC1
2.105776901026963  WC2
1.849675531363642  WC3
1.597909684363248  WC4
1.359815505730616  WC5
1.146346080930527  WC6
0.9685138098307668  WC7
0.8358225820509718  WC8
0.7551066358618733  WC9
0.7300047735758717  WC10
0.7607639068223405  WC11
0.8439490896431752  WC12
0.9722524697223981  WC13
1.135178068593525  WC14
1.321039867006788  WC15
1.519662123486478  WC16
1.724476077968234  WC17
1.932991660078413  WC18
2.145628027225935  WC19
2.363845846889336  WC20
2.588672077289897  WC21
2.819980011700354  WC22
3.05614808281905  WC23
3.293820005304695  WC24
3.528134254650394  WC25
3.753898699002911  WC26
3.967396163284021  WC27
4.167719607135888  WC28
4.356669198400379  WC29
4.537083201991951  WC30
4.710252514959864  WC31
4.873523791195254  WC32
5.019353221645567  WC33
5.136582770666358  WC34
5.213452968187411  WC35
5.240828113600762  WC36
5.214328539925052  WC37
5.135043175488429  WC38
5.009030433054614  WC39
4.845775256715777  WC40
4.65586670508864  WC41
4.448576013923089  WC42
4.230189415888533  WC43
4.003549122790949  WC44
3.768678781646991  WC45
3.524097024858857  WC46
3.268408810537842  WC47
3.001767453096023  WC48
2.726808717542647  WC49
2.448781443760008  WC50
2.174852677447123  WC51
1.912839596659832  WC52
1.66978819004037  WC53
1.450833847697853  WC54
1.258686599366221  WC55
1.093932842091478  WC56
0.9560962122669101  WC57
0.8450473883048795  WC58
0.7621048812278166  WC59
0.7103049939391708  WC60
0.6938486823989517  WC61
0.7172678161333741  WC62
0.7848199948927166  WC63
0.8998526639736095  WC64
1.06320108827885  WC65
1.27041715081464  WC66
1.5097125675742  WC67
1.76343224968571  WC68
2.01340880729448  WC69
2.246898535174412  WC70
2.459402139592374  WC71
2.653656101855899  WC72
2.836676216715255  WC73
3.016907164476116  WC74
3.202299045732849  WC75
3.39909753960286  WC76
3.610907931530243  WC77
3.837863970351192  WC78
4.07603034603105  WC79
4.317405333974717  WC80
4.551009893526912  WC81
4.76511660089986  WC82
4.949719510694422  WC83
5.09792846827258  WC84
5.205714656586651  WC85
5.270560242453387  WC86
5.29009934114714  WC87
5.261617180717227  WC88
5.182667411979827  WC89
5.052509198775389  WC90
4.873755461705197  WC91
4.65348452947308  WC92
4.403068310540985  WC93
4.136289472824895  WC94
3.866117549095778  WC95
3.601515706963153  WC96
3.345932391600368  WC97
3.098052362135915  WC98
2.853883427299406  WC99
59  IP
# Parameters
20  sig
1  k
0.05  eps
4  br
6  bc
0  t0
# Graphics
0.7071097251469876  rm
-0.4999999999826661  rm
0.4999958366025516  rm
0.7071038372138505  rm
0.5000041633974482  rm
-0.4999999999826661  rm
0  rm
0.7071038372138505  rm
0.7071097251469876  rm
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
-1000   
1000   
0    
0   3DFlag
1   Timeflag
0   Colorflag
0   Type
1  color scale
0   minscale
100   xmax
1   xmin
5.334949970245361   ymax
0.7048584818840027   ymin
12   zmax
-12   zmin
50.5   
0.0202020202020202   
3.019904226064682   
0.431956907337013   
0   
0.08333333333333333   
45   Theta
45   Phi
0    xshft
0    yshft
0    zshft
1   xlo
0.7048584818840027   ylo
0   
-1   
100   xhi
5.334949970245361   yhi
20   
1   
# Transpose variables etc
Wc0
100   n columns
1   n rows
1   row skip
1   col skip
1   row 0
# Coupling stuff for H funs
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Array plot stuff
Wc0
100   NCols
0   Row 1
201   NRows
24   RowSkip
0.5520036220550537  Zmin
5.452836513519287  Zmax
# Torus information 
0    Torus flag 1=ON
6.283185307179586  Torus period
# Range information
sig
-1   eq-range stab col
0   shoot flag 1=on
10   eq-range steps
0  eq_range low
1  eq_range high
WA0
WA0
20   Range steps
0   Cycle color 1=on
1   Reset data 1=on
1   Use old I.C.s 1=yes
0  Par1 low
0  Par2 low
1  Par1 high
1  Par2 high
sig
0   BVP side
0   color cycle flag 1=on
10   BVP range steps
0  BVP range low
1  BVP range high
RHS etc ...
WA0(n+1)=WA0+EPS*(ZA-WA0)*H(MX(3)-0)
WA1(n+1)=WA1+EPS*(ZA-WA1)*H(MX(3)-1)
WA2(n+1)=WA2+EPS*(ZA-WA2)*H(MX(3)-2)
WA3(n+1)=WA3+EPS*(ZA-WA3)*H(MX(3)-3)
WA4(n+1)=WA4+EPS*(ZA-WA4)*H(MX(3)-4)
WA5(n+1)=WA5+EPS*(ZA-WA5)*H(MX(3)-5)
WA6(n+1)=WA6+EPS*(ZA-WA6)*H(MX(3)-6)
WA7(n+1)=WA7+EPS*(ZA-WA7)*H(MX(3)-7)
WA8(n+1)=WA8+EPS*(ZA-WA8)*H(MX(3)-8)
WA9(n+1)=WA9+EPS*(ZA-WA9)*H(MX(3)-9)
WA10(n+1)=WA10+EPS*(ZA-WA10)*H(MX(3)-10)
WA11(n+1)=WA11+EPS*(ZA-WA11)*H(MX(3)-11)
WA12(n+1)=WA12+EPS*(ZA-WA12)*H(MX(3)-12)
WA13(n+1)=WA13+EPS*(ZA-WA13)*H(MX(3)-13)
WA14(n+1)=WA14+EPS*(ZA-WA14)*H(MX(3)-14)
WA15(n+1)=WA15+EPS*(ZA-WA15)*H(MX(3)-15)
WA16(n+1)=WA16+EPS*(ZA-WA16)*H(MX(3)-16)
WA17(n+1)=WA17+EPS*(ZA-WA17)*H(MX(3)-17)
WA18(n+1)=WA18+EPS*(ZA-WA18)*H(MX(3)-18)
WA19(n+1)=WA19+EPS*(ZA-WA19)*H(MX(3)-19)
WA20(n+1)=WA20+EPS*(ZA-WA20)*H(MX(3)-20)
WA21(n+1)=WA21+EPS*(ZA-WA21)*H(MX(3)-21)
WA22(n+1)=WA22+EPS*(ZA-WA22)*H(MX(3)-22)
WA23(n+1)=WA23+EPS*(ZA-WA23)*H(MX(3)-23)
WA24(n+1)=WA24+EPS*(ZA-WA24)*H(MX(3)-24)
WA25(n+1)=WA25+EPS*(ZA-WA25)*H(MX(3)-25)
WA26(n+1)=WA26+EPS*(ZA-WA26)*H(MX(3)-26)
WA27(n+1)=WA27+EPS*(ZA-WA27)*H(MX(3)-27)
WA28(n+1)=WA28+EPS*(ZA-WA28)*H(MX(3)-28)
WA29(n+1)=WA29+EPS*(ZA-WA29)*H(MX(3)-29)
WA30(n+1)=WA30+EPS*(ZA-WA30)*H(MX(3)-30)
WA31(n+1)=WA31+EPS*(ZA-WA31)*H(MX(3)-31)
WA32(n+1)=WA32+EPS*(ZA-WA32)*H(MX(3)-32)
WA33(n+1)=WA33+EPS*(ZA-WA33)*H(MX(3)-33)
WA34(n+1)=WA34+EPS*(ZA-WA34)*H(MX(3)-34)
WA35(n+1)=WA35+EPS*(ZA-WA35)*H(MX(3)-35)
WA36(n+1)=WA36+EPS*(ZA-WA36)*H(MX(3)-36)
WA37(n+1)=WA37+EPS*(ZA-WA37)*H(MX(3)-37)
WA38(n+1)=WA38+EPS*(ZA-WA38)*H(MX(3)-38)
WA39(n+1)=WA39+EPS*(ZA-WA39)*H(MX(3)-39)
WA40(n+1)=WA40+EPS*(ZA-WA40)*H(MX(3)-40)
WA41(n+1)=WA41+EPS*(ZA-WA41)*H(MX(3)-41)
WA42(n+1)=WA42+EPS*(ZA-WA42)*H(MX(3)-42)
WA43(n+1)=WA43+EPS*(ZA-WA43)*H(MX(3)-43)
WA44(n+1)=WA44+EPS*(ZA-WA44)*H(MX(3)-44)
WA45(n+1)=WA45+EPS*(ZA-WA45)*H(MX(3)-45)
WA46(n+1)=WA46+EPS*(ZA-WA46)*H(MX(3)-46)
WA47(n+1)=WA47+EPS*(ZA-WA47)*H(MX(3)-47)
WA48(n+1)=WA48+EPS*(ZA-WA48)*H(MX(3)-48)
WA49(n+1)=WA49+EPS*(ZA-WA49)*H(MX(3)-49)
WA50(n+1)=WA50+EPS*(ZA-WA50)*H(MX(3)-50)
WA51(n+1)=WA51+EPS*(ZA-WA51)*H(MX(3)-51)
WA52(n+1)=WA52+EPS*(ZA-WA52)*H(MX(3)-52)
WA53(n+1)=WA53+EPS*(ZA-WA53)*H(MX(3)-53)
WA54(n+1)=WA54+EPS*(ZA-WA54)*H(MX(3)-54)
WA55(n+1)=WA55+EPS*(ZA-WA55)*H(MX(3)-55)
WA56(n+1)=WA56+EPS*(ZA-WA56)*H(MX(3)-56)
WA57(n+1)=WA57+EPS*(ZA-WA57)*H(MX(3)-57)
WA58(n+1)=WA58+EPS*(ZA-WA58)*H(MX(3)-58)
WA59(n+1)=WA59+EPS*(ZA-WA59)*H(MX(3)-59)
WA60(n+1)=WA60+EPS*(ZA-WA60)*H(MX(3)-60)
WA61(n+1)=WA61+EPS*(ZA-WA61)*H(MX(3)-61)
WA62(n+1)=WA62+EPS*(ZA-WA62)*H(MX(3)-62)
WA63(n+1)=WA63+EPS*(ZA-WA63)*H(MX(3)-63)
WA64(n+1)=WA64+EPS*(ZA-WA64)*H(MX(3)-64)
WA65(n+1)=WA65+EPS*(ZA-WA65)*H(MX(3)-65)
WA66(n+1)=WA66+EPS*(ZA-WA66)*H(MX(3)-66)
WA67(n+1)=WA67+EPS*(ZA-WA67)*H(MX(3)-67)
WA68(n+1)=WA68+EPS*(ZA-WA68)*H(MX(3)-68)
WA69(n+1)=WA69+EPS*(ZA-WA69)*H(MX(3)-69)
WA70(n+1)=WA70+EPS*(ZA-WA70)*H(MX(3)-70)
WA71(n+1)=WA71+EPS*(ZA-WA71)*H(MX(3)-71)
WA72(n+1)=WA72+EPS*(ZA-WA72)*H(MX(3)-72)
WA73(n+1)=WA73+EPS*(ZA-WA73)*H(MX(3)-73)
WA74(n+1)=WA74+EPS*(ZA-WA74)*H(MX(3)-74)
WA75(n+1)=WA75+EPS*(ZA-WA75)*H(MX(3)-75)
WA76(n+1)=WA76+EPS*(ZA-WA76)*H(MX(3)-76)
WA77(n+1)=WA77+EPS*(ZA-WA77)*H(MX(3)-77)
WA78(n+1)=WA78+EPS*(ZA-WA78)*H(MX(3)-78)
WA79(n+1)=WA79+EPS*(ZA-WA79)*H(MX(3)-79)
WA80(n+1)=WA80+EPS*(ZA-WA80)*H(MX(3)-80)
WA81(n+1)=WA81+EPS*(ZA-WA81)*H(MX(3)-81)
WA82(n+1)=WA82+EPS*(ZA-WA82)*H(MX(3)-82)
WA83(n+1)=WA83+EPS*(ZA-WA83)*H(MX(3)-83)
WA84(n+1)=WA84+EPS*(ZA-WA84)*H(MX(3)-84)
WA85(n+1)=WA85+EPS*(ZA-WA85)*H(MX(3)-85)
WA86(n+1)=WA86+EPS*(ZA-WA86)*H(MX(3)-86)
WA87(n+1)=WA87+EPS*(ZA-WA87)*H(MX(3)-87)
WA88(n+1)=WA88+EPS*(ZA-WA88)*H(MX(3)-88)
WA89(n+1)=WA89+EPS*(ZA-WA89)*H(MX(3)-89)
WA90(n+1)=WA90+EPS*(ZA-WA90)*H(MX(3)-90)
WA91(n+1)=WA91+EPS*(ZA-WA91)*H(MX(3)-91)
WA92(n+1)=WA92+EPS*(ZA-WA92)*H(MX(3)-92)
WA93(n+1)=WA93+EPS*(ZA-WA93)*H(MX(3)-93)
WA94(n+1)=WA94+EPS*(ZA-WA94)*H(MX(3)-94)
WA95(n+1)=WA95+EPS*(ZA-WA95)*H(MX(3)-95)
WA96(n+1)=WA96+EPS*(ZA-WA96)*H(MX(3)-96)
WA97(n+1)=WA97+EPS*(ZA-WA97)*H(MX(3)-97)
WA98(n+1)=WA98+EPS*(ZA-WA98)*H(MX(3)-98)
WA99(n+1)=WA99+EPS*(ZA-WA99)*H(MX(3)-99)
WB0(n+1)=WB0+EPS*(ZB-WB0)*H(MX(3)-0)
WB1(n+1)=WB1+EPS*(ZB-WB1)*H(MX(3)-1)
WB2(n+1)=WB2+EPS*(ZB-WB2)*H(MX(3)-2)
WB3(n+1)=WB3+EPS*(ZB-WB3)*H(MX(3)-3)
WB4(n+1)=WB4+EPS*(ZB-WB4)*H(MX(3)-4)
WB5(n+1)=WB5+EPS*(ZB-WB5)*H(MX(3)-5)
WB6(n+1)=WB6+EPS*(ZB-WB6)*H(MX(3)-6)
WB7(n+1)=WB7+EPS*(ZB-WB7)*H(MX(3)-7)
WB8(n+1)=WB8+EPS*(ZB-WB8)*H(MX(3)-8)
WB9(n+1)=WB9+EPS*(ZB-WB9)*H(MX(3)-9)
WB10(n+1)=WB10+EPS*(ZB-WB10)*H(MX(3)-10)
WB11(n+1)=WB11+EPS*(ZB-WB11)*H(MX(3)-11)
WB12(n+1)=WB12+EPS*(ZB-WB12)*H(MX(3)-12)
WB13(n+1)=WB13+EPS*(ZB-WB13)*H(MX(3)-13)
WB14(n+1)=WB14+EPS*(ZB-WB14)*H(MX(3)-14)
WB15(n+1)=WB15+EPS*(ZB-WB15)*H(MX(3)-15)
WB16(n+1)=WB16+EPS*(ZB-WB16)*H(MX(3)-16)
WB17(n+1)=WB17+EPS*(ZB-WB17)*H(MX(3)-17)
WB18(n+1)=WB18+EPS*(ZB-WB18)*H(MX(3)-18)
WB19(n+1)=WB19+EPS*(ZB-WB19)*H(MX(3)-19)
WB20(n+1)=WB20+EPS*(ZB-WB20)*H(MX(3)-20)
WB21(n+1)=WB21+EPS*(ZB-WB21)*H(MX(3)-21)
WB22(n+1)=WB22+EPS*(ZB-WB22)*H(MX(3)-22)
WB23(n+1)=WB23+EPS*(ZB-WB23)*H(MX(3)-23)
WB24(n+1)=WB24+EPS*(ZB-WB24)*H(MX(3)-24)
WB25(n+1)=WB25+EPS*(ZB-WB25)*H(MX(3)-25)
WB26(n+1)=WB26+EPS*(ZB-WB26)*H(MX(3)-26)
WB27(n+1)=WB27+EPS*(ZB-WB27)*H(MX(3)-27)
WB28(n+1)=WB28+EPS*(ZB-WB28)*H(MX(3)-28)
WB29(n+1)=WB29+EPS*(ZB-WB29)*H(MX(3)-29)
WB30(n+1)=WB30+EPS*(ZB-WB30)*H(MX(3)-30)
WB31(n+1)=WB31+EPS*(ZB-WB31)*H(MX(3)-31)
WB32(n+1)=WB32+EPS*(ZB-WB32)*H(MX(3)-32)
WB33(n+1)=WB33+EPS*(ZB-WB33)*H(MX(3)-33)
WB34(n+1)=WB34+EPS*(ZB-WB34)*H(MX(3)-34)
WB35(n+1)=WB35+EPS*(ZB-WB35)*H(MX(3)-35)
WB36(n+1)=WB36+EPS*(ZB-WB36)*H(MX(3)-36)
WB37(n+1)=WB37+EPS*(ZB-WB37)*H(MX(3)-37)
WB38(n+1)=WB38+EPS*(ZB-WB38)*H(MX(3)-38)
WB39(n+1)=WB39+EPS*(ZB-WB39)*H(MX(3)-39)
WB40(n+1)=WB40+EPS*(ZB-WB40)*H(MX(3)-40)
WB41(n+1)=WB41+EPS*(ZB-WB41)*H(MX(3)-41)
WB42(n+1)=WB42+EPS*(ZB-WB42)*H(MX(3)-42)
WB43(n+1)=WB43+EPS*(ZB-WB43)*H(MX(3)-43)
WB44(n+1)=WB44+EPS*(ZB-WB44)*H(MX(3)-44)
WB45(n+1)=WB45+EPS*(ZB-WB45)*H(MX(3)-45)
WB46(n+1)=WB46+EPS*(ZB-WB46)*H(MX(3)-46)
WB47(n+1)=WB47+EPS*(ZB-WB47)*H(MX(3)-47)
WB48(n+1)=WB48+EPS*(ZB-WB48)*H(MX(3)-48)
WB49(n+1)=WB49+EPS*(ZB-WB49)*H(MX(3)-49)
WB50(n+1)=WB50+EPS*(ZB-WB50)*H(MX(3)-50)
WB51(n+1)=WB51+EPS*(ZB-WB51)*H(MX(3)-51)
WB52(n+1)=WB52+EPS*(ZB-WB52)*H(MX(3)-52)
WB53(n+1)=WB53+EPS*(ZB-WB53)*H(MX(3)-53)
WB54(n+1)=WB54+EPS*(ZB-WB54)*H(MX(3)-54)
WB55(n+1)=WB55+EPS*(ZB-WB55)*H(MX(3)-55)
WB56(n+1)=WB56+EPS*(ZB-WB56)*H(MX(3)-56)
WB57(n+1)=WB57+EPS*(ZB-WB57)*H(MX(3)-57)
WB58(n+1)=WB58+EPS*(ZB-WB58)*H(MX(3)-58)
WB59(n+1)=WB59+EPS*(ZB-WB59)*H(MX(3)-59)
WB60(n+1)=WB60+EPS*(ZB-WB60)*H(MX(3)-60)
WB61(n+1)=WB61+EPS*(ZB-WB61)*H(MX(3)-61)
WB62(n+1)=WB62+EPS*(ZB-WB62)*H(MX(3)-62)
WB63(n+1)=WB63+EPS*(ZB-WB63)*H(MX(3)-63)
WB64(n+1)=WB64+EPS*(ZB-WB64)*H(MX(3)-64)
WB65(n+1)=WB65+EPS*(ZB-WB65)*H(MX(3)-65)
WB66(n+1)=WB66+EPS*(ZB-WB66)*H(MX(3)-66)
WB67(n+1)=WB67+EPS*(ZB-WB67)*H(MX(3)-67)
WB68(n+1)=WB68+EPS*(ZB-WB68)*H(MX(3)-68)
WB69(n+1)=WB69+EPS*(ZB-WB69)*H(MX(3)-69)
WB70(n+1)=WB70+EPS*(ZB-WB70)*H(MX(3)-70)
WB71(n+1)=WB71+EPS*(ZB-WB71)*H(MX(3)-71)
WB72(n+1)=WB72+EPS*(ZB-WB72)*H(MX(3)-72)
WB73(n+1)=WB73+EPS*(ZB-WB73)*H(MX(3)-73)
WB74(n+1)=WB74+EPS*(ZB-WB74)*H(MX(3)-74)
WB75(n+1)=WB75+EPS*(ZB-WB75)*H(MX(3)-75)
WB76(n+1)=WB76+EPS*(ZB-WB76)*H(MX(3)-76)
WB77(n+1)=WB77+EPS*(ZB-WB77)*H(MX(3)-77)
WB78(n+1)=WB78+EPS*(ZB-WB78)*H(MX(3)-78)
WB79(n+1)=WB79+EPS*(ZB-WB79)*H(MX(3)-79)
WB80(n+1)=WB80+EPS*(ZB-WB80)*H(MX(3)-80)
WB81(n+1)=WB81+EPS*(ZB-WB81)*H(MX(3)-81)
WB82(n+1)=WB82+EPS*(ZB-WB82)*H(MX(3)-82)
WB83(n+1)=WB83+EPS*(ZB-WB83)*H(MX(3)-83)
WB84(n+1)=WB84+EPS*(ZB-WB84)*H(MX(3)-84)
WB85(n+1)=WB85+EPS*(ZB-WB85)*H(MX(3)-85)
WB86(n+1)=WB86+EPS*(ZB-WB86)*H(MX(3)-86)
WB87(n+1)=WB87+EPS*(ZB-WB87)*H(MX(3)-87)
WB88(n+1)=WB88+EPS*(ZB-WB88)*H(MX(3)-88)
WB89(n+1)=WB89+EPS*(ZB-WB89)*H(MX(3)-89)
WB90(n+1)=WB90+EPS*(ZB-WB90)*H(MX(3)-90)
WB91(n+1)=WB91+EPS*(ZB-WB91)*H(MX(3)-91)
WB92(n+1)=WB92+EPS*(ZB-WB92)*H(MX(3)-92)
WB93(n+1)=WB93+EPS*(ZB-WB93)*H(MX(3)-93)
WB94(n+1)=WB94+EPS*(ZB-WB94)*H(MX(3)-94)
WB95(n+1)=WB95+EPS*(ZB-WB95)*H(MX(3)-95)
WB96(n+1)=WB96+EPS*(ZB-WB96)*H(MX(3)-96)
WB97(n+1)=WB97+EPS*(ZB-WB97)*H(MX(3)-97)
WB98(n+1)=WB98+EPS*(ZB-WB98)*H(MX(3)-98)
WB99(n+1)=WB99+EPS*(ZB-WB99)*H(MX(3)-99)
WC0(n+1)=WC0+EPS*(ZC-WC0)*H(MX(3)-0)
WC1(n+1)=WC1+EPS*(ZC-WC1)*H(MX(3)-1)
WC2(n+1)=WC2+EPS*(ZC-WC2)*H(MX(3)-2)
WC3(n+1)=WC3+EPS*(ZC-WC3)*H(MX(3)-3)
WC4(n+1)=WC4+EPS*(ZC-WC4)*H(MX(3)-4)
WC5(n+1)=WC5+EPS*(ZC-WC5)*H(MX(3)-5)
WC6(n+1)=WC6+EPS*(ZC-WC6)*H(MX(3)-6)
WC7(n+1)=WC7+EPS*(ZC-WC7)*H(MX(3)-7)
WC8(n+1)=WC8+EPS*(ZC-WC8)*H(MX(3)-8)
WC9(n+1)=WC9+EPS*(ZC-WC9)*H(MX(3)-9)
WC10(n+1)=WC10+EPS*(ZC-WC10)*H(MX(3)-10)
WC11(n+1)=WC11+EPS*(ZC-WC11)*H(MX(3)-11)
WC12(n+1)=WC12+EPS*(ZC-WC12)*H(MX(3)-12)
WC13(n+1)=WC13+EPS*(ZC-WC13)*H(MX(3)-13)
WC14(n+1)=WC14+EPS*(ZC-WC14)*H(MX(3)-14)
WC15(n+1)=WC15+EPS*(ZC-WC15)*H(MX(3)-15)
WC16(n+1)=WC16+EPS*(ZC-WC16)*H(MX(3)-16)
WC17(n+1)=WC17+EPS*(ZC-WC17)*H(MX(3)-17)
WC18(n+1)=WC18+EPS*(ZC-WC18)*H(MX(3)-18)
WC19(n+1)=WC19+EPS*(ZC-WC19)*H(MX(3)-19)
WC20(n+1)=WC20+EPS*(ZC-WC20)*H(MX(3)-20)
WC21(n+1)=WC21+EPS*(ZC-WC21)*H(MX(3)-21)
WC22(n+1)=WC22+EPS*(ZC-WC22)*H(MX(3)-22)
WC23(n+1)=WC23+EPS*(ZC-WC23)*H(MX(3)-23)
WC24(n+1)=WC24+EPS*(ZC-WC24)*H(MX(3)-24)
WC25(n+1)=WC25+EPS*(ZC-WC25)*H(MX(3)-25)
WC26(n+1)=WC26+EPS*(ZC-WC26)*H(MX(3)-26)
WC27(n+1)=WC27+EPS*(ZC-WC27)*H(MX(3)-27)
WC28(n+1)=WC28+EPS*(ZC-WC28)*H(MX(3)-28)
WC29(n+1)=WC29+EPS*(ZC-WC29)*H(MX(3)-29)
WC30(n+1)=WC30+EPS*(ZC-WC30)*H(MX(3)-30)
WC31(n+1)=WC31+EPS*(ZC-WC31)*H(MX(3)-31)
WC32(n+1)=WC32+EPS*(ZC-WC32)*H(MX(3)-32)
WC33(n+1)=WC33+EPS*(ZC-WC33)*H(MX(3)-33)
WC34(n+1)=WC34+EPS*(ZC-WC34)*H(MX(3)-34)
WC35(n+1)=WC35+EPS*(ZC-WC35)*H(MX(3)-35)
WC36(n+1)=WC36+EPS*(ZC-WC36)*H(MX(3)-36)
WC37(n+1)=WC37+EPS*(ZC-WC37)*H(MX(3)-37)
WC38(n+1)=WC38+EPS*(ZC-WC38)*H(MX(3)-38)
WC39(n+1)=WC39+EPS*(ZC-WC39)*H(MX(3)-39)
WC40(n+1)=WC40+EPS*(ZC-WC40)*H(MX(3)-40)
WC41(n+1)=WC41+EPS*(ZC-WC41)*H(MX(3)-41)
WC42(n+1)=WC42+EPS*(ZC-WC42)*H(MX(3)-42)
WC43(n+1)=WC43+EPS*(ZC-WC43)*H(MX(3)-43)
WC44(n+1)=WC44+EPS*(ZC-WC44)*H(MX(3)-44)
WC45(n+1)=WC45+EPS*(ZC-WC45)*H(MX(3)-45)
WC46(n+1)=WC46+EPS*(ZC-WC46)*H(MX(3)-46)
WC47(n+1)=WC47+EPS*(ZC-WC47)*H(MX(3)-47)
WC48(n+1)=WC48+EPS*(ZC-WC48)*H(MX(3)-48)
WC49(n+1)=WC49+EPS*(ZC-WC49)*H(MX(3)-49)
WC50(n+1)=WC50+EPS*(ZC-WC50)*H(MX(3)-50)
WC51(n+1)=WC51+EPS*(ZC-WC51)*H(MX(3)-51)
WC52(n+1)=WC52+EPS*(ZC-WC52)*H(MX(3)-52)
WC53(n+1)=WC53+EPS*(ZC-WC53)*H(MX(3)-53)
WC54(n+1)=WC54+EPS*(ZC-WC54)*H(MX(3)-54)
WC55(n+1)=WC55+EPS*(ZC-WC55)*H(MX(3)-55)
WC56(n+1)=WC56+EPS*(ZC-WC56)*H(MX(3)-56)
WC57(n+1)=WC57+EPS*(ZC-WC57)*H(MX(3)-57)
WC58(n+1)=WC58+EPS*(ZC-WC58)*H(MX(3)-58)
WC59(n+1)=WC59+EPS*(ZC-WC59)*H(MX(3)-59)
WC60(n+1)=WC60+EPS*(ZC-WC60)*H(MX(3)-60)
WC61(n+1)=WC61+EPS*(ZC-WC61)*H(MX(3)-61)
WC62(n+1)=WC62+EPS*(ZC-WC62)*H(MX(3)-62)
WC63(n+1)=WC63+EPS*(ZC-WC63)*H(MX(3)-63)
WC64(n+1)=WC64+EPS*(ZC-WC64)*H(MX(3)-64)
WC65(n+1)=WC65+EPS*(ZC-WC65)*H(MX(3)-65)
WC66(n+1)=WC66+EPS*(ZC-WC66)*H(MX(3)-66)
WC67(n+1)=WC67+EPS*(ZC-WC67)*H(MX(3)-67)
WC68(n+1)=WC68+EPS*(ZC-WC68)*H(MX(3)-68)
WC69(n+1)=WC69+EPS*(ZC-WC69)*H(MX(3)-69)
WC70(n+1)=WC70+EPS*(ZC-WC70)*H(MX(3)-70)
WC71(n+1)=WC71+EPS*(ZC-WC71)*H(MX(3)-71)
WC72(n+1)=WC72+EPS*(ZC-WC72)*H(MX(3)-72)
WC73(n+1)=WC73+EPS*(ZC-WC73)*H(MX(3)-73)
WC74(n+1)=WC74+EPS*(ZC-WC74)*H(MX(3)-74)
WC75(n+1)=WC75+EPS*(ZC-WC75)*H(MX(3)-75)
WC76(n+1)=WC76+EPS*(ZC-WC76)*H(MX(3)-76)
WC77(n+1)=WC77+EPS*(ZC-WC77)*H(MX(3)-77)
WC78(n+1)=WC78+EPS*(ZC-WC78)*H(MX(3)-78)
WC79(n+1)=WC79+EPS*(ZC-WC79)*H(MX(3)-79)
WC80(n+1)=WC80+EPS*(ZC-WC80)*H(MX(3)-80)
WC81(n+1)=WC81+EPS*(ZC-WC81)*H(MX(3)-81)
WC82(n+1)=WC82+EPS*(ZC-WC82)*H(MX(3)-82)
WC83(n+1)=WC83+EPS*(ZC-WC83)*H(MX(3)-83)
WC84(n+1)=WC84+EPS*(ZC-WC84)*H(MX(3)-84)
WC85(n+1)=WC85+EPS*(ZC-WC85)*H(MX(3)-85)
WC86(n+1)=WC86+EPS*(ZC-WC86)*H(MX(3)-86)
WC87(n+1)=WC87+EPS*(ZC-WC87)*H(MX(3)-87)
WC88(n+1)=WC88+EPS*(ZC-WC88)*H(MX(3)-88)
WC89(n+1)=WC89+EPS*(ZC-WC89)*H(MX(3)-89)
WC90(n+1)=WC90+EPS*(ZC-WC90)*H(MX(3)-90)
WC91(n+1)=WC91+EPS*(ZC-WC91)*H(MX(3)-91)
WC92(n+1)=WC92+EPS*(ZC-WC92)*H(MX(3)-92)
WC93(n+1)=WC93+EPS*(ZC-WC93)*H(MX(3)-93)
WC94(n+1)=WC94+EPS*(ZC-WC94)*H(MX(3)-94)
WC95(n+1)=WC95+EPS*(ZC-WC95)*H(MX(3)-95)
WC96(n+1)=WC96+EPS*(ZC-WC96)*H(MX(3)-96)
WC97(n+1)=WC97+EPS*(ZC-WC97)*H(MX(3)-97)
WC98(n+1)=WC98+EPS*(ZC-WC98)*H(MX(3)-98)
WC99(n+1)=WC99+EPS*(ZC-WC99)*H(MX(3)-99)
IP(n+1)=MX(3)
ZZA=ZA
ZZB=ZB
ZZC=ZC

where ...
ZA = 2*RAN(1)-1   
ZB = BR*RAN(1)   
ZC = BC*RAN(1)   
D0 = SQRT((WA0-ZA)^2+((WB0-ZB))^2+(WC0-ZC)^2)   
D1 = SQRT((WA1-ZA)^2+((WB1-ZB))^2+(WC1-ZC)^2)   
D2 = SQRT((WA2-ZA)^2+((WB2-ZB))^2+(WC2-ZC)^2)   
D3 = SQRT((WA3-ZA)^2+((WB3-ZB))^2+(WC3-ZC)^2)   
D4 = SQRT((WA4-ZA)^2+((WB4-ZB))^2+(WC4-ZC)^2)   
D5 = SQRT((WA5-ZA)^2+((WB5-ZB))^2+(WC5-ZC)^2)   
D6 = SQRT((WA6-ZA)^2+((WB6-ZB))^2+(WC6-ZC)^2)   
D7 = SQRT((WA7-ZA)^2+((WB7-ZB))^2+(WC7-ZC)^2)   
D8 = SQRT((WA8-ZA)^2+((WB8-ZB))^2+(WC8-ZC)^2)   
D9 = SQRT((WA9-ZA)^2+((WB9-ZB))^2+(WC9-ZC)^2)   
D10 = SQRT((WA10-ZA)^2+((WB10-ZB))^2+(WC10-ZC)^2)   
D11 = SQRT((WA11-ZA)^2+((WB11-ZB))^2+(WC11-ZC)^2)   
D12 = SQRT((WA12-ZA)^2+((WB12-ZB))^2+(WC12-ZC)^2)   
D13 = SQRT((WA13-ZA)^2+((WB13-ZB))^2+(WC13-ZC)^2)   
D14 = SQRT((WA14-ZA)^2+((WB14-ZB))^2+(WC14-ZC)^2)   
D15 = SQRT((WA15-ZA)^2+((WB15-ZB))^2+(WC15-ZC)^2)   
D16 = SQRT((WA16-ZA)^2+((WB16-ZB))^2+(WC16-ZC)^2)   
D17 = SQRT((WA17-ZA)^2+((WB17-ZB))^2+(WC17-ZC)^2)   
D18 = SQRT((WA18-ZA)^2+((WB18-ZB))^2+(WC18-ZC)^2)   
D19 = SQRT((WA19-ZA)^2+((WB19-ZB))^2+(WC19-ZC)^2)   
D20 = SQRT((WA20-ZA)^2+((WB20-ZB))^2+(WC20-ZC)^2)   
D21 = SQRT((WA21-ZA)^2+((WB21-ZB))^2+(WC21-ZC)^2)   
D22 = SQRT((WA22-ZA)^2+((WB22-ZB))^2+(WC22-ZC)^2)   
D23 = SQRT((WA23-ZA)^2+((WB23-ZB))^2+(WC23-ZC)^2)   
D24 = SQRT((WA24-ZA)^2+((WB24-ZB))^2+(WC24-ZC)^2)   
D25 = SQRT((WA25-ZA)^2+((WB25-ZB))^2+(WC25-ZC)^2)   
D26 = SQRT((WA26-ZA)^2+((WB26-ZB))^2+(WC26-ZC)^2)   
D27 = SQRT((WA27-ZA)^2+((WB27-ZB))^2+(WC27-ZC)^2)   
D28 = SQRT((WA28-ZA)^2+((WB28-ZB))^2+(WC28-ZC)^2)   
D29 = SQRT((WA29-ZA)^2+((WB29-ZB))^2+(WC29-ZC)^2)   
D30 = SQRT((WA30-ZA)^2+((WB30-ZB))^2+(WC30-ZC)^2)   
D31 = SQRT((WA31-ZA)^2+((WB31-ZB))^2+(WC31-ZC)^2)   
D32 = SQRT((WA32-ZA)^2+((WB32-ZB))^2+(WC32-ZC)^2)   
D33 = SQRT((WA33-ZA)^2+((WB33-ZB))^2+(WC33-ZC)^2)   
D34 = SQRT((WA34-ZA)^2+((WB34-ZB))^2+(WC34-ZC)^2)   
D35 = SQRT((WA35-ZA)^2+((WB35-ZB))^2+(WC35-ZC)^2)   
D36 = SQRT((WA36-ZA)^2+((WB36-ZB))^2+(WC36-ZC)^2)   
D37 = SQRT((WA37-ZA)^2+((WB37-ZB))^2+(WC37-ZC)^2)   
D38 = SQRT((WA38-ZA)^2+((WB38-ZB))^2+(WC38-ZC)^2)   
D39 = SQRT((WA39-ZA)^2+((WB39-ZB))^2+(WC39-ZC)^2)   
D40 = SQRT((WA40-ZA)^2+((WB40-ZB))^2+(WC40-ZC)^2)   
D41 = SQRT((WA41-ZA)^2+((WB41-ZB))^2+(WC41-ZC)^2)   
D42 = SQRT((WA42-ZA)^2+((WB42-ZB))^2+(WC42-ZC)^2)   
D43 = SQRT((WA43-ZA)^2+((WB43-ZB))^2+(WC43-ZC)^2)   
D44 = SQRT((WA44-ZA)^2+((WB44-ZB))^2+(WC44-ZC)^2)   
D45 = SQRT((WA45-ZA)^2+((WB45-ZB))^2+(WC45-ZC)^2)   
D46 = SQRT((WA46-ZA)^2+((WB46-ZB))^2+(WC46-ZC)^2)   
D47 = SQRT((WA47-ZA)^2+((WB47-ZB))^2+(WC47-ZC)^2)   
D48 = SQRT((WA48-ZA)^2+((WB48-ZB))^2+(WC48-ZC)^2)   
D49 = SQRT((WA49-ZA)^2+((WB49-ZB))^2+(WC49-ZC)^2)   
D50 = SQRT((WA50-ZA)^2+((WB50-ZB))^2+(WC50-ZC)^2)   
D51 = SQRT((WA51-ZA)^2+((WB51-ZB))^2+(WC51-ZC)^2)   
D52 = SQRT((WA52-ZA)^2+((WB52-ZB))^2+(WC52-ZC)^2)   
D53 = SQRT((WA53-ZA)^2+((WB53-ZB))^2+(WC53-ZC)^2)   
D54 = SQRT((WA54-ZA)^2+((WB54-ZB))^2+(WC54-ZC)^2)   
D55 = SQRT((WA55-ZA)^2+((WB55-ZB))^2+(WC55-ZC)^2)   
D56 = SQRT((WA56-ZA)^2+((WB56-ZB))^2+(WC56-ZC)^2)   
D57 = SQRT((WA57-ZA)^2+((WB57-ZB))^2+(WC57-ZC)^2)   
D58 = SQRT((WA58-ZA)^2+((WB58-ZB))^2+(WC58-ZC)^2)   
D59 = SQRT((WA59-ZA)^2+((WB59-ZB))^2+(WC59-ZC)^2)   
D60 = SQRT((WA60-ZA)^2+((WB60-ZB))^2+(WC60-ZC)^2)   
D61 = SQRT((WA61-ZA)^2+((WB61-ZB))^2+(WC61-ZC)^2)   
D62 = SQRT((WA62-ZA)^2+((WB62-ZB))^2+(WC62-ZC)^2)   
D63 = SQRT((WA63-ZA)^2+((WB63-ZB))^2+(WC63-ZC)^2)   
D64 = SQRT((WA64-ZA)^2+((WB64-ZB))^2+(WC64-ZC)^2)   
D65 = SQRT((WA65-ZA)^2+((WB65-ZB))^2+(WC65-ZC)^2)   
D66 = SQRT((WA66-ZA)^2+((WB66-ZB))^2+(WC66-ZC)^2)   
D67 = SQRT((WA67-ZA)^2+((WB67-ZB))^2+(WC67-ZC)^2)   
D68 = SQRT((WA68-ZA)^2+((WB68-ZB))^2+(WC68-ZC)^2)   
D69 = SQRT((WA69-ZA)^2+((WB69-ZB))^2+(WC69-ZC)^2)   
D70 = SQRT((WA70-ZA)^2+((WB70-ZB))^2+(WC70-ZC)^2)   
D71 = SQRT((WA71-ZA)^2+((WB71-ZB))^2+(WC71-ZC)^2)   
D72 = SQRT((WA72-ZA)^2+((WB72-ZB))^2+(WC72-ZC)^2)   
D73 = SQRT((WA73-ZA)^2+((WB73-ZB))^2+(WC73-ZC)^2)   
D74 = SQRT((WA74-ZA)^2+((WB74-ZB))^2+(WC74-ZC)^2)   
D75 = SQRT((WA75-ZA)^2+((WB75-ZB))^2+(WC75-ZC)^2)   
D76 = SQRT((WA76-ZA)^2+((WB76-ZB))^2+(WC76-ZC)^2)   
D77 = SQRT((WA77-ZA)^2+((WB77-ZB))^2+(WC77-ZC)^2)   
D78 = SQRT((WA78-ZA)^2+((WB78-ZB))^2+(WC78-ZC)^2)   
D79 = SQRT((WA79-ZA)^2+((WB79-ZB))^2+(WC79-ZC)^2)   
D80 = SQRT((WA80-ZA)^2+((WB80-ZB))^2+(WC80-ZC)^2)   
D81 = SQRT((WA81-ZA)^2+((WB81-ZB))^2+(WC81-ZC)^2)   
D82 = SQRT((WA82-ZA)^2+((WB82-ZB))^2+(WC82-ZC)^2)   
D83 = SQRT((WA83-ZA)^2+((WB83-ZB))^2+(WC83-ZC)^2)   
D84 = SQRT((WA84-ZA)^2+((WB84-ZB))^2+(WC84-ZC)^2)   
D85 = SQRT((WA85-ZA)^2+((WB85-ZB))^2+(WC85-ZC)^2)   
D86 = SQRT((WA86-ZA)^2+((WB86-ZB))^2+(WC86-ZC)^2)   
D87 = SQRT((WA87-ZA)^2+((WB87-ZB))^2+(WC87-ZC)^2)   
D88 = SQRT((WA88-ZA)^2+((WB88-ZB))^2+(WC88-ZC)^2)   
D89 = SQRT((WA89-ZA)^2+((WB89-ZB))^2+(WC89-ZC)^2)   
D90 = SQRT((WA90-ZA)^2+((WB90-ZB))^2+(WC90-ZC)^2)   
D91 = SQRT((WA91-ZA)^2+((WB91-ZB))^2+(WC91-ZC)^2)   
D92 = SQRT((WA92-ZA)^2+((WB92-ZB))^2+(WC92-ZC)^2)   
D93 = SQRT((WA93-ZA)^2+((WB93-ZB))^2+(WC93-ZC)^2)   
D94 = SQRT((WA94-ZA)^2+((WB94-ZB))^2+(WC94-ZC)^2)   
D95 = SQRT((WA95-ZA)^2+((WB95-ZB))^2+(WC95-ZC)^2)   
D96 = SQRT((WA96-ZA)^2+((WB96-ZB))^2+(WC96-ZC)^2)   
D97 = SQRT((WA97-ZA)^2+((WB97-ZB))^2+(WC97-ZC)^2)   
D98 = SQRT((WA98-ZA)^2+((WB98-ZB))^2+(WC98-ZC)^2)   
D99 = SQRT((WA99-ZA)^2+((WB99-ZB))^2+(WC99-ZC)^2)   

User-defined functions:
H(U) = EXP(-SIG*(1-COS(2*PI*U*K/NTOT)))