sib-yi5b.rr の tryb6() の出力.  (ODE*dt の bug fix 版)
tmpb.py は Data/2021_07_09_tryb6_tmpb.py として保存.
graph は Data/2021_07_09_tryb6_fig1.png

maxLEt= 8.656684050147129e+29
res_1.x= [0.060145405402867516, -1.20074804549872, 9.438660716835336, -29.022737131194667, -46.606911486598264, 587.9594544508735, -808.9438085694859, -6769.394034019031, 35050.636639119104, -66736.73622440665, 37633.340028566956, 29467.77387996437, -6048.122536295094, -16923.59111153512, -12748.901726971548, -6689.469492772272, -2878.1021418826467, -1086.2619274418348, -372.94734545976127, -119.14781206671961, -36.278775153736234, -8.677863612903192, -4.07044156492677, 0.2519178884835531, -0.3881719338093061, -1.1539945701702408, 0.12398100749043622, -0.3155873643236216, -0.0052641099570091615, 0.0067706012834766165]
res_1.cost(loss functionの値)= 3.853369149670236e-07
res_1.optimality= 2.892216991024088e-09
// For asir input
MaxLEt= 8.656684050147129e+29 ;
Fcoef= [0.060145405402867516, -1.20074804549872, 9.438660716835336, -29.022737131194667, -46.606911486598264, 587.9594544508735, -808.9438085694859, -6769.394034019031, 35050.636639119104, -66736.73622440665, 37633.340028566956, 29467.77387996437, -6048.122536295094, -16923.59111153512, -12748.901726971548, -6689.469492772272, -2878.1021418826467, -1086.2619274418348, -372.94734545976127, -119.14781206671961, -36.278775153736234, -8.677863612903192, -4.07044156492677, 0.2519178884835531, -0.3881719338093061, -1.1539945701702408, 0.12398100749043622, -0.3155873643236216, -0.0052641099570091615, 0.0067706012834766165] ;
if (length(Fcoef)>=M) print(show_errors(L,B,Pt,Q,M,EE,T0,T1,H,Fcoef,MaxLEt));
//
QStandardPaths: XDG_RUNTIME_DIR not set, defaulting to '/tmp/runtime-nobuki'
Show relative errors with simulation values
QFileSystemWatcher::removePaths: list is empty
QFileSystemWatcher::removePaths: list is empty

real	0m54.003s
user	0m17.147s
sys	0m0.227s

--- on orange3m  python3
17.090u 0.274s 1:39.22 17.4%	0+0k 98192+0io 235pf+0w
relative error が win wsl ubuntsu より大きい. 原因不明.
------------------------
[[1.00646081599953700578746732469409857865e+79,13430548124698817961.7552819096280274774,[0.000120177830991900032913640950289920188842,0.000257528980893160008324116317200358429638,0.000267115798450930200974397716090608079738,0.000108281251831504054848202803002285956114,0.000153862169044879427884797580968975500709,0.000271142715706900325328042808999067067767,6.36125854393596046254495901248247977848e-05,0.000108422138002485829939047550999334972377,0.000317928325032889031440697569313405968644,8.61598008765699615891886655217614919597e-05]],[int value of (Lf-B)^2,int value of ((Lf-B)/MaxLEt)^2,list of value diffs at given points]]
Main terms of L/MaxLEt is (-4.71561e-52*t^29+2.33584e-51*t^27-2.94135e-51*t^25+1.37367e-51*t^23-5.81702e-52*t^21)*dt^11+(-7.77887e-22*t^35+1.98088e-21*t^33-2.62004e-21*t^31+2.02442e-21*t^29-1.04233e-21*t^27+4.19652e-22*t^25)*dt