
 # subdomains  in x-direction nsx =     1
 # subdomains  in y-direction nsy =     1

 # control subdomain ncontrol =     0

 Global # elements  in x-direction ngx =   400
 Global # elements  in y-direction ngy =   400

 Local # elements  in x-direction nx =   400
 Local # elements  in y-direction ny =   400

  processor grid size in the x-dir  hpx =        200.00000    
  processor grid size in the y-dir  hpy  =        200.00000    

  domain size xsize =       0.20000000E+03 m

  domain size ysize =       0.20000000E+03 m

   maxlow (max.iter. in low freq range ) =   4000 perclow ( perc. of(low) freq. solved with maxlow =       0.50000000E-01

  rounded xsize =       0.20000000E+03 m ; 
 rounded ysize =       0.20000000E+03 m 


  mesh size h=       0.50000000E+00 m

  number of snapshots in freq. domain=     1

  freq. to compute  snapshots in freq.domain
    1         0.30000000E-01 kHz

  number of snapshots in time domain=     4

  times to compute  snapshots in time domain
    1         0.10000000E+02 msec

  times to compute  snapshots in time domain
    2         0.20000000E+02 msec

  times to compute  snapshots in time domain
    3         0.25000000E+02 msec

  times to compute  snapshots in time domain
    4         0.30000000E+02 msec

  quadrature choice=1,standard,=2,mid point
  keyqua =     2

  number of points to compute solution= 
 in the frequency domain    npf  =   200

   key==1 for type 1 source (tsang&rader paper) ;
       key==2 for type 2 source(gaussian derivative)             ;    2



   frequency  freq=       0.70000000E+00 khz



  w0 =       0.43982297E+01  alfa =       0.11060000E+01



  total simulation time  stime =        0.12000000E+03 msec


   parameters of  source band-pass filter

 set equal to 0. if not filter is desired

 f1  =       0.00000000E+00   khz 

 f2  =       0.00000000E+00   khz 

 f3  =       0.15000000E+01   khz 

 f4  =       0.16000000E+01   khz 



   maxiter =  max. # inner iterations  4000

   reduc =  relative error required        0.10000000E-03

  number of materials (nmat) =     5

  mat. number   density (ros)    bulk modulus (xkr)   Shear modulus (xmur)


     1          0.260000E+07     0.133000E+07          0.850000E+06

     2          0.260000E+07     0.133000E+07          0.850000E+06

     3          0.260000E+07     0.133000E+07          0.850000E+06

     4          0.260000E+07     0.133000E+07          0.850000E+06

     5          0.260000E+07     0.133000E+07          0.850000E+06


  mat. number  deltp (p-log.decr.)     delts (s-log.decr.)   


     1          0.000000E+00          0.000000E+00

     2          0.000000E+00          0.000000E+00

     3          0.000000E+00          0.000000E+00

     4          0.000000E+00          0.000000E+00

     5          0.000000E+00          0.000000E+00


  mat. number   density (ros)    bulk modulus (xkr)   Shear modulus (xmur)


     1          0.260000E+07     0.133000E+07          0.850000E+06

     2          0.260000E+07     0.133000E+07          0.850000E+06

     3          0.260000E+07     0.133000E+07          0.850000E+06

     4          0.260000E+07     0.133000E+07          0.850000E+06

     5          0.260000E+07     0.133000E+07          0.850000E+06


  bulmin-umin-bulmax-umax   

     0.13300E+07    0.85000E+06    0.13300E+07    0.85000E+06


 minimum  phase pp velocity   =  0.97336E+00 m/msec

 maximum  phase pp velocity  =  0.97336E+00 m/msec



 minimum  phase ss velocity   =  0.57177E+00 m/msec

 maximum  phase ss velocity  =  0.57177E+00 m/msec



 phase  velocity  pp1 =  0.97336E+00 m/msec

 phase velocity pp2 =  0.97336E+00 m/msec

 phase velocity ss1 =  0.57177E+00 m/msec

 phase velocity ss2 =  0.57177E+00 m/msec



   mat. number      fq1                     fq2    (Liu model)  


     1        0.10000000E-06 kHz       0.10000000E+08 kHz

     2        0.10000000E-06 kHz       0.10000000E+08 kHz

     3        0.10000000E-06 kHz       0.10000000E+08 kHz

     4        0.10000000E-06 kHz       0.10000000E+08 kHz

     5        0.10000000E-06 kHz       0.10000000E+08 kHz

  material number     qp(i)  (Liu model)


     1              0.00000000E+00

     2              0.00000000E+00

     3              0.00000000E+00

     4              0.00000000E+00

     5              0.00000000E+00

  material number     const(i)       tau1(i)(msec)          tau2 (i)(msec) 


     1             0.000000E+00      0.159155E+07           0.159155E-07

     2             0.000000E+00      0.159155E+07           0.159155E-07

     3             0.000000E+00      0.159155E+07           0.159155E-07

     4             0.000000E+00      0.159155E+07           0.159155E-07

     5             0.000000E+00      0.159155E+07           0.159155E-07

 Material  Number               Ref. Speed ve    relaxed  plane wave modulus xmr


            1                   0.97336E+00         0.24633E+07

            2                   0.97336E+00         0.24633E+07

            3                   0.97336E+00         0.24633E+07

            4                   0.97336E+00         0.24633E+07

            5                   0.97336E+00         0.24633E+07

  time step  dt =       0.51368287E-01 msec



  number of time steps  nt =   2336



  kdel=0 : delta, kdel = 1:  d/dx(delta). kdel =      1

   mute in msec  tmute =       0.00000000E+00
 mute in samples   ntmute =     0


  jsou =     5  ksou =     5

  xsou =       0.20000000E+01 m ysou =       0.20000000E+01 m

  number of receivers, nrec =    80
  xinic    2.0000000000000000     

  yinic    2.0000000000000000     

  xfin    190.00000000000000     

  yfin    2.0000000000000000     

  deltarecx    2.3797468354430378     

  deltarecy    0.0000000000000000     


 receiver i =    1   irx & iry =     5    5
               rx(i) & ry(i) =   2.0000000      2.0000000      m

 receiver i =    2   irx & iry =     9    5
               rx(i) & ry(i) =   4.3797468      2.0000000      m

 receiver i =    3   irx & iry =    14    5
               rx(i) & ry(i) =   6.7594937      2.0000000      m

 receiver i =    4   irx & iry =    19    5
               rx(i) & ry(i) =   9.1392405      2.0000000      m

 receiver i =    5   irx & iry =    24    5
               rx(i) & ry(i) =   11.518987      2.0000000      m

 receiver i =    6   irx & iry =    28    5
               rx(i) & ry(i) =   13.898734      2.0000000      m

 receiver i =    7   irx & iry =    33    5
               rx(i) & ry(i) =   16.278481      2.0000000      m

 receiver i =    8   irx & iry =    38    5
               rx(i) & ry(i) =   18.658228      2.0000000      m

 receiver i =    9   irx & iry =    43    5
               rx(i) & ry(i) =   21.037975      2.0000000      m

 receiver i =   10   irx & iry =    47    5
               rx(i) & ry(i) =   23.417722      2.0000000      m

 receiver i =   11   irx & iry =    52    5
               rx(i) & ry(i) =   25.797468      2.0000000      m

 receiver i =   12   irx & iry =    57    5
               rx(i) & ry(i) =   28.177215      2.0000000      m

 receiver i =   13   irx & iry =    62    5
               rx(i) & ry(i) =   30.556962      2.0000000      m

 receiver i =   14   irx & iry =    66    5
               rx(i) & ry(i) =   32.936709      2.0000000      m

 receiver i =   15   irx & iry =    71    5
               rx(i) & ry(i) =   35.316456      2.0000000      m

 receiver i =   16   irx & iry =    76    5
               rx(i) & ry(i) =   37.696203      2.0000000      m

 receiver i =   17   irx & iry =    81    5
               rx(i) & ry(i) =   40.075949      2.0000000      m

 receiver i =   18   irx & iry =    85    5
               rx(i) & ry(i) =   42.455696      2.0000000      m

 receiver i =   19   irx & iry =    90    5
               rx(i) & ry(i) =   44.835443      2.0000000      m

 receiver i =   20   irx & iry =    95    5
               rx(i) & ry(i) =   47.215190      2.0000000      m

 receiver i =   21   irx & iry =   100    5
               rx(i) & ry(i) =   49.594937      2.0000000      m

 receiver i =   22   irx & iry =   104    5
               rx(i) & ry(i) =   51.974684      2.0000000      m

 receiver i =   23   irx & iry =   109    5
               rx(i) & ry(i) =   54.354430      2.0000000      m

 receiver i =   24   irx & iry =   114    5
               rx(i) & ry(i) =   56.734177      2.0000000      m

 receiver i =   25   irx & iry =   119    5
               rx(i) & ry(i) =   59.113924      2.0000000      m

 receiver i =   26   irx & iry =   123    5
               rx(i) & ry(i) =   61.493671      2.0000000      m

 receiver i =   27   irx & iry =   128    5
               rx(i) & ry(i) =   63.873418      2.0000000      m

 receiver i =   28   irx & iry =   133    5
               rx(i) & ry(i) =   66.253165      2.0000000      m

 receiver i =   29   irx & iry =   138    5
               rx(i) & ry(i) =   68.632911      2.0000000      m

 receiver i =   30   irx & iry =   143    5
               rx(i) & ry(i) =   71.012658      2.0000000      m

 receiver i =   31   irx & iry =   147    5
               rx(i) & ry(i) =   73.392405      2.0000000      m

 receiver i =   32   irx & iry =   152    5
               rx(i) & ry(i) =   75.772152      2.0000000      m

 receiver i =   33   irx & iry =   157    5
               rx(i) & ry(i) =   78.151899      2.0000000      m

 receiver i =   34   irx & iry =   162    5
               rx(i) & ry(i) =   80.531646      2.0000000      m

 receiver i =   35   irx & iry =   166    5
               rx(i) & ry(i) =   82.911392      2.0000000      m

 receiver i =   36   irx & iry =   171    5
               rx(i) & ry(i) =   85.291139      2.0000000      m

 receiver i =   37   irx & iry =   176    5
               rx(i) & ry(i) =   87.670886      2.0000000      m

 receiver i =   38   irx & iry =   181    5
               rx(i) & ry(i) =   90.050633      2.0000000      m

 receiver i =   39   irx & iry =   185    5
               rx(i) & ry(i) =   92.430380      2.0000000      m

 receiver i =   40   irx & iry =   190    5
               rx(i) & ry(i) =   94.810127      2.0000000      m

 receiver i =   41   irx & iry =   195    5
               rx(i) & ry(i) =   97.189873      2.0000000      m

 receiver i =   42   irx & iry =   200    5
               rx(i) & ry(i) =   99.569620      2.0000000      m

 receiver i =   43   irx & iry =   204    5
               rx(i) & ry(i) =   101.94937      2.0000000      m

 receiver i =   44   irx & iry =   209    5
               rx(i) & ry(i) =   104.32911      2.0000000      m

 receiver i =   45   irx & iry =   214    5
               rx(i) & ry(i) =   106.70886      2.0000000      m

 receiver i =   46   irx & iry =   219    5
               rx(i) & ry(i) =   109.08861      2.0000000      m

 receiver i =   47   irx & iry =   223    5
               rx(i) & ry(i) =   111.46835      2.0000000      m

 receiver i =   48   irx & iry =   228    5
               rx(i) & ry(i) =   113.84810      2.0000000      m

 receiver i =   49   irx & iry =   233    5
               rx(i) & ry(i) =   116.22785      2.0000000      m

 receiver i =   50   irx & iry =   238    5
               rx(i) & ry(i) =   118.60759      2.0000000      m

 receiver i =   51   irx & iry =   242    5
               rx(i) & ry(i) =   120.98734      2.0000000      m

 receiver i =   52   irx & iry =   247    5
               rx(i) & ry(i) =   123.36709      2.0000000      m

 receiver i =   53   irx & iry =   252    5
               rx(i) & ry(i) =   125.74684      2.0000000      m

 receiver i =   54   irx & iry =   257    5
               rx(i) & ry(i) =   128.12658      2.0000000      m

 receiver i =   55   irx & iry =   262    5
               rx(i) & ry(i) =   130.50633      2.0000000      m

 receiver i =   56   irx & iry =   266    5
               rx(i) & ry(i) =   132.88608      2.0000000      m

 receiver i =   57   irx & iry =   271    5
               rx(i) & ry(i) =   135.26582      2.0000000      m

 receiver i =   58   irx & iry =   276    5
               rx(i) & ry(i) =   137.64557      2.0000000      m

 receiver i =   59   irx & iry =   281    5
               rx(i) & ry(i) =   140.02532      2.0000000      m

 receiver i =   60   irx & iry =   285    5
               rx(i) & ry(i) =   142.40506      2.0000000      m

 receiver i =   61   irx & iry =   290    5
               rx(i) & ry(i) =   144.78481      2.0000000      m

 receiver i =   62   irx & iry =   295    5
               rx(i) & ry(i) =   147.16456      2.0000000      m

 receiver i =   63   irx & iry =   300    5
               rx(i) & ry(i) =   149.54430      2.0000000      m

 receiver i =   64   irx & iry =   304    5
               rx(i) & ry(i) =   151.92405      2.0000000      m

 receiver i =   65   irx & iry =   309    5
               rx(i) & ry(i) =   154.30380      2.0000000      m

 receiver i =   66   irx & iry =   314    5
               rx(i) & ry(i) =   156.68354      2.0000000      m

 receiver i =   67   irx & iry =   319    5
               rx(i) & ry(i) =   159.06329      2.0000000      m

 receiver i =   68   irx & iry =   323    5
               rx(i) & ry(i) =   161.44304      2.0000000      m

 receiver i =   69   irx & iry =   328    5
               rx(i) & ry(i) =   163.82278      2.0000000      m

 receiver i =   70   irx & iry =   333    5
               rx(i) & ry(i) =   166.20253      2.0000000      m

 receiver i =   71   irx & iry =   338    5
               rx(i) & ry(i) =   168.58228      2.0000000      m

 receiver i =   72   irx & iry =   342    5
               rx(i) & ry(i) =   170.96203      2.0000000      m

 receiver i =   73   irx & iry =   347    5
               rx(i) & ry(i) =   173.34177      2.0000000      m

 receiver i =   74   irx & iry =   352    5
               rx(i) & ry(i) =   175.72152      2.0000000      m

 receiver i =   75   irx & iry =   357    5
               rx(i) & ry(i) =   178.10127      2.0000000      m

 receiver i =   76   irx & iry =   361    5
               rx(i) & ry(i) =   180.48101      2.0000000      m

 receiver i =   77   irx & iry =   366    5
               rx(i) & ry(i) =   182.86076      2.0000000      m

 receiver i =   78   irx & iry =   371    5
               rx(i) & ry(i) =   185.24051      2.0000000      m

 receiver i =   79   irx & iry =   376    5
               rx(i) & ry(i) =   187.62025      2.0000000      m

 receiver i =   80   irx & iry =   380    5
               rx(i) & ry(i) =   190.00000      2.0000000      m

  nyquist frequency =       0.97336319E+01 khz

  dt  =       0.51368287E-01  msec

 delf =       0.80000000E-02  khz

 delw =       0.50265482E-01  2 pi*khz



   parameters of  source band-pass filter

 set equal to 0. if not filter is desired

 f1  =       0.00000000E+00   khz 

 f2  =       0.00000000E+00   khz 

 f3  =       0.15000000E+01   khz 

 f4  =       0.16000000E+01   khz 



  freq. indices for snapshots 
    1      4

  time steps for snapshots 
    1    194

  time steps for snapshots 
    2    389

  time steps for snapshots 
    3    486

  time steps for snapshots 
    4    584

  if1 if2 if3 if4 =     1     1   188   201

 jw1 - jw2 - jwlow       1    200     10

  pulse length in samples     105

 psi=      0.39200000E+01width=      0.71428571E+00msec

 ts=      0.17857143E+01msec. totpul=      0.53571429E+01msec

  pulse length in samples    2336
  writes in 20 source in time
  writes in 21 amplitude source
  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.00000000E+00 m  
 up to yf                     =       0.35000000E+02 m 
 material type number  numat  =     1

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =     1
 up to yf  (y-global grid #)    =    70
 material type number  numat    =     1

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.35000000E+02 m  
 up to yf                     =       0.65000000E+02 m 
 material type number  numat  =     3

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =    71
 up to yf  (y-global grid #)    =   130
 material type number  numat    =     3

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.65000000E+02 m  
 up to yf                     =       0.71000000E+02 m 
 material type number  numat  =     4

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   131
 up to yf  (y-global grid #)    =   142
 material type number  numat    =     4

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.71000000E+02 m  
 up to yf                     =       0.12500000E+03 m 
 material type number  numat  =     3

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   143
 up to yf  (y-global grid #)    =   250
 material type number  numat    =     3

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.12500000E+03 m  
 up to yf                     =       0.13100000E+03 m 
 material type number  numat  =     5

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   251
 up to yf  (y-global grid #)    =   262
 material type number  numat    =     5

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.13100000E+03 m  
 up to yf                     =       0.16500000E+03 m 
 material type number  numat  =     3

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   263
 up to yf  (y-global grid #)    =   330
 material type number  numat    =     3

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.16500000E+03 m  
 up to yf                     =       0.20000000E+03 m 
 material type number  numat  =     2

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   331
 up to yf  (y-global grid #)    =   400
 material type number  numat    =     2

  from  xi                     =      -0.10000000E+01 m  
 up to xf                     =      -0.10000000E+01 m 
 from  yi                     =      -0.10000000E+01 m  
 up to yf                     =      -0.10000000E+01 m 
 material type number  numat  =    -1


  jw1 , jw2  in main =     1      200

  frequency number   200 frequency :       0.15960000E+01
 jw  real sour, imag sour =    200     -0.53985E-03      0.39211E-03

  global error =        0.8448732E-02 iter. number =   100 denomin. =        0.3412115E-04

  global error =        0.6842986E-02 iter. number =   200 denomin. =        0.4716775E-04

  global error =        0.5427133E-02 iter. number =   300 denomin. =        0.6045024E-04

  global error =        0.4823721E-02 iter. number =   400 denomin. =        0.7273025E-04

  global error =        0.4169478E-02 iter. number =   500 denomin. =        0.8639193E-04

  global error =        0.2659167E-02 iter. number =   600 denomin. =        0.9588958E-04

  global error =        0.1745674E-02 iter. number =   700 denomin. =        0.9841321E-04

  global error =        0.6678606E-03 iter. number =   800 denomin. =        0.9981875E-04

  global error =        0.5885918E-03 iter. number =   900 denomin. =        0.1004619E-03

  global error =        0.5941031E-03 iter. number =  1000 denomin. =        0.1006395E-03

  global error =        0.5289645E-03 iter. number =  1100 denomin. =        0.1008089E-03

  global error =        0.3335453E-03 iter. number =  1200 denomin. =        0.1011929E-03

  global error =        0.1874430E-03 iter. number =  1300 denomin. =        0.1016489E-03

  global error =        0.1163631E-03 iter. number =  1400 denomin. =        0.1017399E-03

  global error =        0.1007112E-03 iter. number =  1500 denomin. =        0.1016467E-03

  convergence  for freq. #   200 achieved in  1505 inner iterations 
 relative error value =       0.99866329E-04


  frequency number   199 frequency :       0.15880000E+01
 jw  real sour, imag sour =    199     -0.12154E-02      0.72596E-03

  global error =        0.4656672E-02 iter. number =   100 denomin. =        0.2241465E-03

  global error =        0.3686921E-02 iter. number =   200 denomin. =        0.2685923E-03

  global error =        0.2996516E-02 iter. number =   300 denomin. =        0.3118308E-03

  global error =        0.2628593E-02 iter. number =   400 denomin. =        0.3545788E-03

  global error =        0.2269972E-02 iter. number =   500 denomin. =        0.4016694E-03

  global error =        0.1425812E-02 iter. number =   600 denomin. =        0.4437397E-03

  global error =        0.9211428E-03 iter. number =   700 denomin. =        0.4742278E-03

  global error =        0.3979722E-03 iter. number =   800 denomin. =        0.4969506E-03

  global error =        0.3433793E-03 iter. number =   900 denomin. =        0.5166181E-03

  global error =        0.3187141E-03 iter. number =  1000 denomin. =        0.5327511E-03

  global error =        0.2683547E-03 iter. number =  1100 denomin. =        0.5466909E-03

  global error =        0.1664571E-03 iter. number =  1200 denomin. =        0.5553350E-03

  global error =        0.1004046E-03 iter. number =  1300 denomin. =        0.5607215E-03

  convergence  for freq. #   199 achieved in  1301 inner iterations 
 relative error value =       0.99456960E-04


  frequency number   198 frequency :       0.15800000E+01
 jw  real sour, imag sour =    198     -0.20292E-02      0.97691E-03

  global error =        0.4488558E-02 iter. number =   100 denomin. =        0.8080931E-03

  global error =        0.3462452E-02 iter. number =   200 denomin. =        0.9021540E-03

  global error =        0.2917606E-02 iter. number =   300 denomin. =        0.9758201E-03

  global error =        0.2447628E-02 iter. number =   400 denomin. =        0.1036645E-02

  global error =        0.1975572E-02 iter. number =   500 denomin. =        0.1094620E-02

  global error =        0.1193529E-02 iter. number =   600 denomin. =        0.1122358E-02

  global error =        0.7542314E-03 iter. number =   700 denomin. =        0.1120771E-02

  global error =        0.4919524E-03 iter. number =   800 denomin. =        0.1082204E-02

  global error =        0.4298299E-03 iter. number =   900 denomin. =        0.1054192E-02

  global error =        0.3533353E-03 iter. number =  1000 denomin. =        0.1041070E-02

  global error =        0.2677501E-03 iter. number =  1100 denomin. =        0.1034360E-02

  global error =        0.1697622E-03 iter. number =  1200 denomin. =        0.1032235E-02

  global error =        0.1230427E-03 iter. number =  1300 denomin. =        0.1030826E-02

  global error =        0.1044217E-03 iter. number =  1400 denomin. =        0.1032055E-02

  convergence  for freq. #   198 achieved in  1431 inner iterations 
 relative error value =       0.99887297E-04


  frequency number   197 frequency :       0.15720000E+01
 jw  real sour, imag sour =    197     -0.29806E-02      0.11183E-02

  global error =        0.4121444E-02 iter. number =   100 denomin. =        0.1436400E-02

  global error =        0.3296099E-02 iter. number =   200 denomin. =        0.1628234E-02

  global error =        0.2687333E-02 iter. number =   300 denomin. =        0.1877268E-02

  global error =        0.2256409E-02 iter. number =   400 denomin. =        0.2150472E-02

  global error =        0.1848072E-02 iter. number =   500 denomin. =        0.2447420E-02

  global error =        0.1198802E-02 iter. number =   600 denomin. =        0.2679096E-02

  global error =        0.7721286E-03 iter. number =   700 denomin. =        0.2828130E-02

  global error =        0.3585814E-03 iter. number =   800 denomin. =        0.2877607E-02

  global error =        0.3174552E-03 iter. number =   900 denomin. =        0.2893055E-02

  global error =        0.2914166E-03 iter. number =  1000 denomin. =        0.2897977E-02
