The code runs at: 2025-07-09 19:22:43
Baseflow: solved with bvp4c @ eta=[0.0, 10.0000], y=[0.0, 10.0000] 
Baseflow: extracted at i=   1, Re= 378.30, x= 378.30 
Baseflow: extracted at i=  11, Re= 428.07, x= 484.39 
Baseflow: extracted at i=  21, Re= 472.63, x= 590.47 
Baseflow: extracted at i=  31, Re= 513.33, x= 696.56 
Baseflow: extracted at i=  41, Re= 551.04, x= 802.64 
Baseflow: extracted at i=  51, Re= 586.32, x= 908.73 
Baseflow: extracted at i=  61, Re= 619.60, x=1014.81 
Baseflow: extracted at i=  71, Re= 651.18, x=1120.90 
Baseflow: extracted at i=  81, Re= 681.30, x=1226.98 
Baseflow: extracted at i=  91, Re= 710.14, x=1333.07 
Baseflow: extracted at i= 101, Re= 737.86, x=1439.15 
Baseflow: extracted at i= 111, Re= 764.57, x=1545.24 
Baseflow: extracted at i= 121, Re= 790.38, x=1651.32 
Baseflow: extracted at i= 131, Re= 815.37, x=1757.41 
Baseflow: extracted at i= 141, Re= 839.62, x=1863.49 
Baseflow: extracted at i= 151, Re= 863.19, x=1969.58 
Baseflow: extracted at i= 161, Re= 886.13, x=2075.66 
Baseflow: extracted at i= 171, Re= 908.49, x=2181.75 
Baseflow: extracted at i= 181, Re= 930.32, x=2287.83 
Baseflow: extracted at i= 191, Re= 951.64, x=2393.92 
Baseflow: extracted at i= 201, Re= 972.50, x=2500.00 

NPSE: initialization 
NPSE strategy: Standard 
Mode (1,-1), amp=2.5e-06, a= 4.414e-02 +6.299e-03 error=1.4e-02 
Mode (1,+1), amp=2.5e-06, a= 4.414e-02 +6.299e-03 error=1.4e-02 
Mode (2,-1), amp=2.5e-06, a= 9.065e-02 +3.046e-03 error=2.6e-02 
Mode (2,+0), amp=7.6e-03, a= 1.049e-01 +3.187e-03 error=2.5e-02 
Mode (2,+1), amp=2.5e-06, a= 9.065e-02 +3.046e-03 error=2.6e-02 

NPSE: Marching start. 
...Strategy: Fmn is NOT iterated. 
NPSE: station 2/201, Re=383.6/972.5 (7.58 % active)
  Fmn(4,+0)=2.1e-05,  Fmn(0,+0)=5.9e-06,  Fmn(4,-1)=1.9e-08,
  Fmn(4,+1)=1.9e-08,  Fmn(0,-1)=9.6e-09,  Fmn(0,+1)=9.6e-09,
  Fmn(3,-1)=8.1e-09,  Fmn(5,+1)=8.1e-09,  Fmn(3,+1)=8.1e-09,

  res_a(2,+0)=8.33e-04, res_n=0.00e+00, a=1.0446192e-01 2.4943477e-03 i 
  res_a(2,+0)=1.57e-04, res_n=0.00e+00, a=1.0452649e-01 2.6375640e-03 i 
  res_a(2,+0)=2.97e-05, res_n=0.00e+00, a=1.0451903e-01 2.6087763e-03 i 
  res_a(2,+0)=5.63e-06, res_n=0.00e+00, a=1.0451951e-01 2.6143889e-03 i 
  res_a(2,+0)=1.07e-06, res_n=0.00e+00, a=1.0451960e-01 2.6133248e-03 i 
  res_a(2,+0)=2.03e-07, res_n=0.00e+00, a=1.0451955e-01 2.6135210e-03 i 
  res_a(2,+0)=3.84e-08, res_n=0.00e+00, a=1.0451956e-01 2.6134859e-03 i 
  res_a(2,+0)=7.29e-09, res_n=0.00e+00, a=1.0451956e-01 2.6134920e-03 i 
  res_a(2,+0)=1.38e-09, res_n=0.00e+00, a=1.0451956e-01 2.6134910e-03 i 
  res_a(2,+0)=2.62e-10, res_n=0.00e+00, a=1.0451956e-01 2.6134911e-03 i 
  res_a(2,+0)=4.95e-11, res_n=0.00e+00, a=1.0451956e-01 2.6134911e-03 i 
  res_a(2,+0)=9.36e-12, res_n=0.00e+00, a=1.0451956e-01 2.6134911e-03 i 
  res_a(2,+0)=1.77e-12, res_n=0.00e+00, a=1.0451956e-01 2.6134911e-03 i 
  res_a(2,+0)=3.34e-13, res_n=0.00e+00, a=1.0451956e-01 2.6134911e-03 i 
  Mode(0 ,-2) Started 
  Mode(0 ,-1) Started 
  Mode(0 ,0) Started 
  Mode(0 ,1) Started 
  Mode(0 ,2) Started 
  Mode(1 ,-2) Started 
  Mode(1 ,0) Started 
  Mode(1 ,2) Started 
  Mode(2 ,-2) Started 
  Mode(2 ,2) Started 
  Mode(3 ,-2) Started 
  Mode(3 ,-1) Started 
  Mode(3 ,0) Started 
  Mode(3 ,1) Started 
  Mode(3 ,2) Started 
  Mode(4 ,-2) Started 
  Mode(4 ,-1) Started 
  Mode(4 ,0) Started 
  Mode(4 ,1) Started 
  Mode(4 ,2) Started 
  Mode(5 ,-2) Started 
  Mode(5 ,-1) Started 
  Mode(5 ,0) Started 
  Mode(5 ,1) Started 
  Mode(5 ,2) Started 

NPSE: station 3/201, Re=388.8/972.5 (45.45 % active)
  Fmn(4,+0)=8.4e-06,  Fmn(0,+0)=5.6e-06,  Fmn(4,-1)=1.6e-08,
  Fmn(4,+1)=1.6e-08,  Fmn(0,-1)=8.8e-09,  Fmn(0,+1)=8.8e-09,
  Fmn(3,-1)=7.9e-09,  Fmn(5,+1)=7.9e-09,  Fmn(3,+1)=7.9e-09,

  res_a(2,+0)=2.21e-04, res_n=0.00e+00, a=1.0422052e-01 2.2340189e-03 i 
  res_a(2,+0)=4.17e-05, res_n=0.00e+00, a=1.0420694e-01 2.1946376e-03 i 
  res_a(2,+0)=7.85e-06, res_n=0.00e+00, a=1.0420823e-01 2.2023859e-03 i 
  res_a(2,+0)=1.48e-06, res_n=0.00e+00, a=1.0420823e-01 2.2009035e-03 i 
  res_a(2,+0)=2.80e-07, res_n=0.00e+00, a=1.0420818e-01 2.2011795e-03 i 
  res_a(2,+0)=5.29e-08, res_n=0.00e+00, a=1.0420820e-01 2.2011295e-03 i 
  res_a(2,+0)=1.00e-08, res_n=0.00e+00, a=1.0420820e-01 2.2011383e-03 i 
  res_a(2,+0)=1.89e-09, res_n=0.00e+00, a=1.0420820e-01 2.2011368e-03 i 
  res_a(2,+0)=3.57e-10, res_n=0.00e+00, a=1.0420820e-01 2.2011370e-03 i 
  res_a(2,+0)=6.74e-11, res_n=0.00e+00, a=1.0420820e-01 2.2011370e-03 i 
  res_a(2,+0)=1.27e-11, res_n=0.00e+00, a=1.0420820e-01 2.2011370e-03 i 
  res_a(2,+0)=2.39e-12, res_n=0.00e+00, a=1.0420820e-01 2.2011370e-03 i 
  res_a(2,+0)=4.49e-13, res_n=0.00e+00, a=1.0420820e-01 2.2011370e-03 i 

NPSE: station 4/201, Re=393.9/972.5 (45.45 % active)
  Fmn(4,+0)=2.5e-05,  Fmn(0,+0)=5.3e-06,  Fmn(2,+0)=7.3e-07,
  Fmn(4,-1)=1.1e-08,  Fmn(4,+1)=1.1e-08,  Fmn(0,-1)=9.1e-09,
  Fmn(0,+1)=9.1e-09,  Fmn(3,+1)=6.9e-09,  Fmn(5,-1)=6.9e-09,

  res_a(2,+0)=5.89e-05, res_n=0.00e+00, a=1.0390550e-01 1.8470836e-03 i 
  res_a(2,+0)=1.11e-05, res_n=0.00e+00, a=1.0390570e-01 1.8359743e-03 i 
  res_a(2,+0)=2.10e-06, res_n=0.00e+00, a=1.0390532e-01 1.8380360e-03 i 
  res_a(2,+0)=3.96e-07, res_n=0.00e+00, a=1.0390546e-01 1.8376636e-03 i 
  res_a(2,+0)=7.47e-08, res_n=0.00e+00, a=1.0390542e-01 1.8377289e-03 i 
  res_a(2,+0)=1.41e-08, res_n=0.00e+00, a=1.0390543e-01 1.8377179e-03 i 
  res_a(2,+0)=2.66e-09, res_n=0.00e+00, a=1.0390543e-01 1.8377197e-03 i 
  res_a(2,+0)=5.01e-10, res_n=0.00e+00, a=1.0390543e-01 1.8377194e-03 i 
  res_a(2,+0)=9.42e-11, res_n=0.00e+00, a=1.0390543e-01 1.8377194e-03 i 
  res_a(2,+0)=1.77e-11, res_n=0.00e+00, a=1.0390543e-01 1.8377194e-03 i 
  res_a(2,+0)=3.32e-12, res_n=0.00e+00, a=1.0390543e-01 1.8377194e-03 i 
  res_a(2,+0)=6.22e-13, res_n=0.00e+00, a=1.0390543e-01 1.8377194e-03 i 
  Mode(1 ,-5) Started 
  Mode(1 ,-3) Started 
  Mode(1 ,3) Started 
  Mode(1 ,5) Started 
  Mode(2 ,-5) Started 
  Mode(2 ,-3) Started 
  Mode(2 ,3) Started 
  Mode(2 ,5) Started 
  Mode(3 ,-5) Started 
  Mode(3 ,-3) Started 
  Mode(3 ,3) Started 
  Mode(3 ,5) Started 
  Mode(4 ,-5) Started 
  Mode(4 ,-3) Started 
  Mode(4 ,3) Started 
  Mode(4 ,5) Started 
  Mode(5 ,-5) Started 
  Mode(5 ,-3) Started 
  Mode(5 ,3) Started 
  Mode(5 ,5) Started 

NPSE: station 5/201, Re=399.0/972.5 (75.76 % active)
  Fmn(4,+0)=2.2e-05,  Fmn(0,+0)=5.0e-06,  Fmn(2,+0)=1.2e-06,
  Fmn(4,-1)=1.5e-08,  Fmn(4,+1)=1.5e-08,  Fmn(0,-1)=8.4e-09,
  Fmn(0,+1)=8.4e-09,  Fmn(3,+1)=7.1e-09,  Fmn(5,-1)=7.1e-09,

  res_a(2,+0)=1.76e-04, res_n=0.00e+00, a=1.0371957e-01 1.6065009e-03 i 
  res_a(2,+0)=3.32e-05, res_n=0.00e+00, a=1.0370187e-01 1.5783888e-03 i 
  res_a(2,+0)=6.26e-06, res_n=0.00e+00, a=1.0370431e-01 1.5841534e-03 i 
  res_a(2,+0)=1.18e-06, res_n=0.00e+00, a=1.0370403e-01 1.5830055e-03 i 
  res_a(2,+0)=2.23e-07, res_n=0.00e+00, a=1.0370405e-01 1.5832281e-03 i 
  res_a(2,+0)=4.23e-08, res_n=0.00e+00, a=1.0370405e-01 1.5831860e-03 i 
  res_a(2,+0)=8.00e-09, res_n=0.00e+00, a=1.0370405e-01 1.5831938e-03 i 
  res_a(2,+0)=1.52e-09, res_n=0.00e+00, a=1.0370405e-01 1.5831924e-03 i 
  res_a(2,+0)=2.87e-10, res_n=0.00e+00, a=1.0370405e-01 1.5831926e-03 i 
  res_a(2,+0)=5.44e-11, res_n=0.00e+00, a=1.0370405e-01 1.5831926e-03 i 
  res_a(2,+0)=1.03e-11, res_n=0.00e+00, a=1.0370405e-01 1.5831926e-03 i 
  res_a(2,+0)=1.94e-12, res_n=0.00e+00, a=1.0370405e-01 1.5831926e-03 i 
  res_a(2,+0)=3.66e-13, res_n=0.00e+00, a=1.0370405e-01 1.5831926e-03 i 
  Mode(0 ,-5) Started 
  Mode(0 ,-3) Started 
  Mode(0 ,3) Started 
  Mode(0 ,5) Started 

NPSE: station 6/201, Re=404.0/972.5 (81.82 % active)
  Fmn(0,+0)=4.9e-06,  Fmn(4,+0)=4.7e-06,  Fmn(2,+0)=1.8e-06,
  Fmn(4,-1)=1.1e-08,  Fmn(4,+1)=1.1e-08,  Fmn(0,-1)=7.5e-09,
  Fmn(0,+1)=7.5e-09,  Fmn(3,-1)=6.5e-09,  Fmn(5,+1)=6.5e-09,

  res_a(2,+0)=2.03e-04, res_n=0.00e+00, a=1.0332418e-01 1.2313723e-03 i 
  res_a(2,+0)=3.78e-05, res_n=0.00e+00, a=1.0335412e-01 1.2544619e-03 i 
  res_a(2,+0)=7.05e-06, res_n=0.00e+00, a=1.0334929e-01 1.2493256e-03 i 
  res_a(2,+0)=1.32e-06, res_n=0.00e+00, a=1.0335003e-01 1.2504158e-03 i 
  res_a(2,+0)=2.46e-07, res_n=0.00e+00, a=1.0334992e-01 1.2501926e-03 i 
  res_a(2,+0)=4.61e-08, res_n=0.00e+00, a=1.0334994e-01 1.2502369e-03 i 
  res_a(2,+0)=8.64e-09, res_n=0.00e+00, a=1.0334994e-01 1.2502284e-03 i 
  res_a(2,+0)=1.62e-09, res_n=0.00e+00, a=1.0334994e-01 1.2502300e-03 i 
  res_a(2,+0)=3.05e-10, res_n=0.00e+00, a=1.0334994e-01 1.2502297e-03 i 
  res_a(2,+0)=5.73e-11, res_n=0.00e+00, a=1.0334994e-01 1.2502297e-03 i 
  res_a(2,+0)=1.08e-11, res_n=0.00e+00, a=1.0334994e-01 1.2502297e-03 i 
  res_a(2,+0)=2.02e-12, res_n=0.00e+00, a=1.0334994e-01 1.2502297e-03 i 
  res_a(2,+0)=3.79e-13, res_n=0.00e+00, a=1.0334994e-01 1.2502297e-03 i 

NPSE: station 7/201, Re=408.9/972.5 (81.82 % active)
{操作在以下过程中被用户终止 <a href="matlab:matlab.internal.language.introspective.errorDocCallback('NPSE_Dvector')" style="font-weight:bold">NPSE_Dvector</a>


位置 <a href="matlab:matlab.internal.language.introspective.errorDocCallback('NPSE_Fmn', 'C:\Users\admin\Desktop\code\NPSE_Fmn.m', 38)" style="font-weight:bold">NPSE_Fmn</a> (<a href="matlab: opentoline('C:\Users\admin\Desktop\code\NPSE_Fmn.m',38,0)">第 38 行</a>)
        NPSE_Dvector(xi,m,n,Mesh,Case,alf,Fai);

位置 <a href="matlab:matlab.internal.language.introspective.errorDocCallback('main_secondary_Ma0', 'C:\Users\admin\Desktop\code\main_secondary_Ma0.m', 237)" style="font-weight:bold">main_secondary_Ma0</a> (<a href="matlab: opentoline('C:\Users\admin\Desktop\code\main_secondary_Ma0.m',237,0)">第 237 行</a>)
                    [Fmn] = NPSE_Fmn(xi,Case,Mesh,Baseflow,alf,Fai,modes);
} 
