state_space_matrices.cpp

#include "state_space_matrices.h"
#include <algorithm>
#include <iostream>

// Function to linearly interpolate between two values
double interpolate(double x1, double x2, double y1, double y2, double x) {
    return y1 + (x - x1) * (y2 - y1) / (x2 - x1);
}

void interpolateMatrices(const std::vector<double>& A1, const std::vector<double>& A2, const std::vector<double>& B1, const std::vector<double>& B2, double ratio, std::vector<double>& A, std::vector<double>& B) {
    size_t sizeA = A1.size();
    size_t sizeB = B1.size();

    A.resize(sizeA);
    B.resize(sizeB);

    for (size_t i = 0; i < sizeA; ++i) {
        A[i] = A1[i] * (1 - ratio) + A2[i] * ratio;
    }

    for (size_t i = 0; i < sizeB; ++i) {
        B[i] = B1[i] * (1 - ratio) + B2[i] * ratio;
    }
}

void getStateSpaceMatrices(double speed, std::vector<double>& A, std::vector<double>& B) {
    // Define matrices for specific speeds
    std::vector<double> A5 = {
        -9.61335633071093e-13, 4.99999999991687, 1.00000000005217, 1.12249207912694e-10,
        1.45212601715332e-12, -5.13981176158834e-14, -1.17800813245335e-10, 1.00000000101190,
        0.000840272043916469, -3.04085570568224e-07, -0.101668288331192, -2.70007795164826,
        4.37929332288747e-06, -0.0149473707639360, -0.00572337288021148, -0.244530705224609
    };
    std::vector<double> B5 = {
        -1.69674752574204e-12, -1.10484134201392e-17, 2.75706811981441e-12,
        -7.22760359972355e-13, -6.12418309548146e-18, -4.20492416927654e-12,
        -0.00261112713870843, -2.21991902099372e-06, 0.00164151357821695,
        0.000330132364452263, 2.40305021375118e-16, 7.27066721427474e-05
    };

    std::vector<double> A4_95 = {
        -9.01283788107966e-13, 4.94999999991807, 1.00000000005150, 1.02225340570838e-10,
        1.37317803068419e-12, -2.86969784271035e-14, -1.15281382177819e-10, 1.00000000098300,
        0.000823550970431623, -3.04087182572284e-07, -0.100503174078730, -2.67724319766960,
        4.29214637054594e-06, -0.0149473707644157, -0.00567344878374934, -0.241977131658948
    };
    std::vector<double> B4_95 = {
        -1.61006977645766e-12, -1.06544146598581e-17, 2.66967681685285e-12,
        -6.91978221078683e-13, -5.99632733161866e-18, -4.07505080975970e-12,
        -0.00255916571644820, -2.21991902104794e-06, 0.00160884742992723,
        0.000323562575167797, 2.35330901458283e-16, 7.12597969421443e-05
    };

    std::vector<double> A4_90 = {
        -8.51536884045135e-13, 4.89999999991898, 1.00000000005069, 9.23009772626311e-11,
        1.29580609987106e-12, 3.25171493261161e-13, -1.12789854274247e-10, 1.00000000095469,
        0.000806997947184077, -3.04180718538697e-07, -0.0993363431449647, -2.65444934899444,
        4.20587515693409e-06, -0.0149473708068564, -0.00562387064064128, -0.239408859052273
    };
    std::vector<double> B4_90 = {
        -1.52697588496482e-12, -1.02706000003297e-17, 2.58420943064636e-12,
        -6.62192528869874e-13, -5.87004517738357e-18, -3.94788388991593e-12,
        -0.00250772651913888, -2.21991902111140e-06, 0.00157650958580106,
        0.000317058813935240, 2.30411006763557e-16, 6.98274633723678e-05
    };

    std::vector<double> A4_85 = {
        -8.09411890578019e-13, 4.84999999991995, 1.00000000004998, 8.21180526120949e-11,
        1.22008348438622e-12, -1.37884025164158e-14, -1.10325748017571e-10, 1.00000000092697,
        0.000790612974199054, -3.04090189686791e-07, -0.0981776113377637, -2.63147696377808,
        4.12047968321949e-06, -0.0149473707650344, -0.00557377565075599, -0.236857591329661
    };
    std::vector<double> B4_85 = {
        -1.44744228451458e-12, -9.89522983762018e-18, 2.50063587458466e-12,
        -6.33383318384642e-13, -5.74480330111493e-18, -3.82339603364036e-12,
        -0.00245680954679795, -2.21991902118488e-06, 0.00154450004580024,
        0.000310621080764036, 2.25531887560148e-16, 6.84096714324699e-05
    };

    std::vector<double> A4_80 = {
        -7.34385696386519e-13, 4.79999999992032, 1.00000000004926, 7.13630831093136e-11,
        1.14604906589345e-12, -5.64980035888962e-15, -1.07890003608118e-10, 1.00000000089993,
        0.000774396051483462, -3.04091482463981e-07, -0.0970240931186455, -2.60832406365490,
        4.03595995063572e-06, -0.0149473707655361, -0.00552314966151221, -0.234319612513039
    };
    std::vector<double> B4_80 = {
        -1.37128170056585e-12, -9.53234680342427e-18, 2.41892416470366e-12,
        -6.05527557051140e-13, -5.62086441546745e-18, -3.70155881450731e-12,
        -0.00240641479940192, -2.21991902123147e-06, 0.00151281880986238,
        0.000304249375650748, 2.20706376118664e-16, 6.70064211183090e-05
    };

    std::vector<double> A4_75 = {
        -7.10500661197627e-13, 4.74999999992122, 1.00000000004852, 6.04081737911871e-11,
        1.07368047033182e-12, -1.16644335498305e-14, -1.05482470828804e-10, 1.00000000087342,
        0.000758347179076054, -3.04092730803388e-07, -0.0958711576940379, -2.58514989886509,
        3.95231596029939e-06, -0.0149473707657502, -0.00547262415438826, -0.231775268864480
    };
    std::vector<double> B4_75 = {
        -1.29840427824320e-12, -9.17921355699394e-18, 2.33907680822088e-12,
        -5.78604020085229e-13, -5.49844249091608e-18, -3.58234465312399e-12,
        -0.00235654227702166, -2.21991902130067e-06, 0.00148146587798612,
        0.000297943698586811, 2.15933612235164e-16, 6.56177124244894e-05
    };

    std::vector<double> A4_70 = {
        -6.36946054900172e-13, 4.69999999992179, 1.00000000004776, 5.09670001074028e-11,
        1.00270246715853e-12, 3.19865131867704e-13, -1.03101737580554e-10, 1.00000000084756,
        0.000742466356983197, -3.04176814318684e-07, -0.0947242840277329, -2.56184437538917,
        3.86954771359133e-06, -0.0149473708064573, -0.00542176745050590, -0.229241212198348
    };
    std::vector<double> B4_70 = {
        -1.22866925778143e-12, -8.83417269357864e-18, 2.26104071311120e-12,
        -5.52590706050294e-13, -5.37715381153143e-18, -3.46572532635260e-12,
        -0.00230719197964646, -2.21991902136798e-06, 0.00145044125011801,
        0.000291704049581678, 2.11204883035566e-16, 6.42435453498958e-05
    };

    std::vector<double> A4_65 = {
        -6.08075351441128e-13, 4.64999999992280, 1.00000000004696, 4.48498325364050e-11,
        9.33912512656165e-13, 3.36074554266266e-13, -1.00751139580150e-10, 1.00000000082220,
        0.000726753585226173, -3.04175781502065e-07, -0.0935755333308265, -2.53857473388569,
        3.78765521173668e-06, -0.0149473708059661, -0.00537123612402460, -0.226692664118000
    };
    std::vector<double> B4_65 = {
        -1.16198449140292e-12, -8.49917273872745e-18, 2.18481022119452e-12,
        -5.27467946313950e-13, -5.25724224944738e-18, -3.35167264493667e-12,
        -0.00225836390727870, -2.21991902141998e-06, 0.00141974492621474,
        0.000285530428621664, 2.06532091399675e-16, 6.28839198879649e-05
    };

    std::vector<double> A4_60 = {
        -5.44038431403657e-13, 4.59999999992374, 1.00000000004625, 3.28672582471383e-11,
        8.66407071770115e-13, -1.10541214032152e-14, -9.84250250815014e-11, 1.00000000079745,
        0.000711208863827049, -3.04095794639885e-07, -0.0924332958216358, -2.51510038398956,
        3.70663845546712e-06, -0.0149473707666786, -0.00532007914061475, -0.224161192842129
    };
    std::vector<double> B4_60 = {
        -1.09828884498768e-12, -8.16746445207661e-18, 2.11036127673514e-12,
        -5.03214340399976e-13, -5.13899319882119e-18, -3.24015886007684e-12,
        -0.00221005805993403, -2.21991902145577e-06, 0.00138937690624019,
        0.000279422835723957, 2.01910787377551e-16, 6.15388360393821e-05
    };

    std::vector<double> A4_55 = {
        -5.09863129663126e-13, 4.54999999992488, 1.00000000004549, 2.76491854268955e-11,
        8.00922220930493e-13, 1.61063026509714e-14, -9.61295947027062e-11, 1.00000000077322,
        0.000695832192815888, -3.04096738984186e-07, -0.0912927609142237, -2.49161503832512,
        3.62649744678456e-06, -0.0149473707668653, -0.00526906639904487, -0.221623997982454
    };
    std::vector<double> B4_55 = {
        -1.03741838699244e-12, -7.85243500456000e-18, 2.03766521613812e-12,
        -4.79810762013297e-13, -5.02163555100139e-18, -3.13115598580454e-12,
        -0.00216227443766832, -2.21991902151617e-06, 0.00135933719018618,
        0.000273381270872339, 1.97336892838584e-16, 6.02082937970477e-05
    };

    std::vector<double> A4_50 = {
        -4.80907710152366e-13, 4.49999999992566, 1.00000000004466, 1.98170032110000e-11,
        7.36698094957368e-13, 1.55569229749384e-14, -9.38663634155683e-11, 1.00000000074948,
        0.000680623572193629, -3.04097641262347e-07, -0.0901539016336029, -2.46802488856239,
        3.54723218644193e-06, -0.0149473707672943, -0.00521782001340744, -0.219088950021781
    };
    std::vector<double> B4_50 = {
        -9.79314866929540e-13, -7.54115076169229e-18, 1.96669842412490e-12,
        -4.57234223814677e-13, -4.90582309308642e-18, -3.02463520605888e-12,
        -0.00211501304045585, -2.21991902156017e-06, 0.00132962577799247,
        0.000267405734076291, 1.92821780687276e-16, 5.88922931597468e-05
    };

    std::vector<double> A4_45 = {
        -4.29849801081019e-13, 4.44999999992607, 1.00000000004388, 1.36214722022456e-11,
        6.73804154281341e-13, 3.26261824043706e-13, -9.16258932726659e-11, 1.00000000072632,
        0.000665583002001082, -3.04170528559918e-07, -0.0890180073966229, -2.44436172685189,
        3.46884267583229e-06, -0.0149473708055136, -0.00516647079518563, -0.216555104316174
    };
    std::vector<double> B4_45 = {
        -9.23856580226852e-13, -7.23634041402454e-18, 1.89743545582575e-12,
        -4.35468795752644e-13, -4.79143117422719e-18, -2.92056897532845e-12,
        -0.00206827386833815, -2.21991902161591e-06, 0.00130024266964242,
        0.000261496225329907, 1.88356691533841e-16, 5.75908341228137e-05
    };

    std::vector<double> A4_40 = {
        -3.96300209902776e-13, 4.39999999992696, 1.00000000004308, 6.75432380297426e-12,
        6.12932002182860e-13, 1.05403093769439e-14, -8.94142980076183e-11, 1.00000000070365,
        0.000650710482247267, -3.04098867051907e-07, -0.0878862342780060, -2.42060656193699,
        3.39132891588047e-06, -0.0149473707677530, -0.00511494473156342, -0.214025617157660
    };
    std::vector<double> B4_40 = {
        -8.70983933992964e-13, -6.94494367347455e-18, 1.82986086135969e-12,
        -4.14492205354249e-13, -4.67831133911401e-18, -2.81892929070398e-12,
        -0.00202205692131589, -2.21991902166944e-06, 0.00127118786509411,
        0.000255652744627290, 1.83939548767301e-16, 5.63039166816053e-05
    };

    std::vector<double> A4_35 = {
        -3.59862967864474e-13, 4.34999999992815, 1.00000000004230, 1.01423051463133e-12,
        5.53673115712191e-13, 1.77503306331594e-14, -8.72300751969881e-11, 1.00000000068151,
        0.000636006012946361, -3.04099452877514e-07, -0.0867581066630601, -2.39677478526517,
        3.31469090823753e-06, -0.0149473707680722, -0.00506330279492391, -0.211498551285175
    };
    std::vector<double> B4_35 = {
        -8.20535754564521e-13, -6.65563641028675e-18, 1.76395166919321e-12,
        -3.94285483090783e-13, -4.56650892455699e-18, -2.71968738694661e-12,
        -0.00197636219938931, -2.21991902170547e-06, 0.00124246136430606,
        0.000249875291977952, 1.79570550052882e-16, 5.50315408348316e-05
    };

    std::vector<double> A4_30 = {
        -3.14630068374561e-13, 4.29999999992847, 1.00000000004153, -5.90399432581650e-12,
        4.95557188797085e-13, 3.26480275543406e-13, -8.50750592359144e-11, 1.00000000065993,
        0.000621469594131326, -3.04166147958872e-07, -0.0856301146747578, -2.37289325863506,
        3.23892865390931e-06, -0.0149473708047130, -0.00501164565327045, -0.208966919787890
    };
    std::vector<double> B4_30 = {
        -7.72548340226338e-13, -6.37358336485659e-18, 1.69967365919524e-12,
        -3.74828052129556e-13, -4.45621423824834e-18, -2.62281473598666e-12,
        -0.00193118970259903, -2.21991902175395e-06, 0.00121406316726283,
        0.000244163867376011, 1.75259927394926e-16, 5.37737065777762e-05
    };

    std::vector<double> A4_25 = {
        -2.95597045179901e-13, 4.24999999992948, 1.00000000004068, -1.16153789267665e-11,
        4.39442764307309e-13, 3.15983403211076e-13, -8.29502104296538e-11, 1.00000000063875,
        0.000607101225814062, -3.04164752228469e-07, -0.0845067921089167, -2.34889236252383,
        3.16404215427862e-06, -0.0149473708046275, -0.00495970261955136, -0.206442691225694
    };
    std::vector<double> B4_25 = {
        -7.26801795813780e-13, -6.10907703224775e-18, 1.63701300778807e-12,
        -3.56106112938662e-13, -4.34724647914621e-18, -2.52828379231269e-12,
        -0.00188653943095817, -2.21991902181593e-06, 0.00118599327393210,
        0.000238518470823313, 1.70998656434083e-16, 5.25304139075445e-05
    };

    std::vector<double> A4_20 = {
        -2.69739387081996e-13, 4.19999999993043, 1.00000000003996, -1.49106900902859e-11,
        3.84707942824626e-13, 1.52964498469088e-13, -8.08512956654831e-11, 1.00000000061810,
        0.000592900908011078, -3.04132021001925e-07, -0.0833891833350070, -2.32472452239075,
        3.09003141043808e-06, -0.0149473707864835, -0.00490728431762176, -0.203931280768879
    };
    std::vector<double> B4_20 = {
        -6.83307372451525e-13, -5.84249341409771e-18, 1.57594315201848e-12,
        -3.38094211333534e-13, -4.23944380793372e-18, -2.43606575680449e-12,
        -0.00184241138443889, -2.21991902184347e-06, 0.00115825168425643,
        0.000232939102311442, 1.66788347053030e-16, 5.13016628189262e-05
    };

    std::vector<double> A4_15 = {
        -2.26365385238438e-13, 4.14999999993133, 1.00000000004173, -1.06243914476169e-10,
        3.32616778123344e-13, 9.94947331116554e-14, -8.59273705540049e-11, 1.00000000082165,
        0.000578772704645207, -3.04094164614538e-07, -0.0822698685070506, -2.30062406106911,
        3.03948596057035e-06, -0.0149473687422088, -0.00485532734510313, -0.201403175081947
    };
    std::vector<double> B4_15 = {
        -4.43518062630060e-13, -6.10326770258080e-18, 1.51145561995067e-12,
        -8.24593727351371e-13, -4.48839667346175e-18, -2.35327112165317e-12,
        -0.00179880501898382, -2.21991902181777e-06, 0.00113083839966913,
        0.000227425600125897, 1.76642697500876e-16, 5.00874534311768e-05
    };

    std::vector<double> A4_10 = {
        -1.71591729937411e-13, 4.09999999993196, 1.00000000003823, -2.30902769467705e-11,
        2.79724046994734e-13, 3.32780503772012e-14, -7.67347513014687e-11, 1.00000000057834,
        0.000565004424037502, -3.04101391139061e-07, -0.0811543491409987, -2.27642039118906,
        2.94463719576559e-06, -0.0149473707693059, -0.00480313024544702, -0.198879492417194
    };
    std::vector<double> B4_10 = {
        -6.02613806557536e-13, -5.33640976278357e-18, 1.45849664675212e-12,
        -3.04137679561917e-13, -4.02800129771091e-18, -2.25845519586168e-12,
        -0.00175572196693653, -2.21991902193962e-06, 0.00110375341582392,
        0.000221978449431704, 1.58523435700598e-16, 4.88877853773449e-05
    };

    std::vector<double> A4_05 = {
        -1.25351206124835e-13, 4.04999999993229, 1.00000000003759, -3.18079604862804e-11,
        2.29222475413386e-13, 1.67582076309644e-13, -7.47238416904712e-11, 1.00000000055911,
        0.000551308257899748, -3.04130048551607e-07, -0.0800452469222076, -2.25200998101794,
        2.87325372729240e-06, -0.0149473707863342, -0.00475030782566760, -0.196373149448497
    };
    std::vector<double> B4_05 = {
        -5.65240437285463e-13, -5.10129748475905e-18, 1.40207923802994e-12,
        -2.88156709810724e-13, -3.92442588855310e-18, -2.17300573287236e-12,
        -0.00171316059593679, -2.21991902195780e-06, 0.00107699673697938,
        0.000216597165056060, 1.54469623271095e-16, 4.77026590160665e-05
    };

    std::vector<double> A4_00 = {
        -1.20130999577785e-13, 3.99999999993383, 1.00000000003665, -3.13606971749740e-11,
        1.81285027547738e-13, 3.17606899670818e-13, -7.27340267136783e-11, 1.00000000054033,
        0.000537780142357921, -3.04156771644445e-07, -0.0789370670217614, -2.22759187666301,
        2.80274601999740e-06, -0.0149473708031882, -0.00469764127356482, -0.193859754669469
    };
    std::vector<double> B4_00 = {
        -5.29785003373010e-13, -4.85979511513432e-18, 1.34715406419159e-12,
        -2.72817168125866e-13, -3.82230972321250e-18, -2.08975554010860e-12,
        -0.00167112145016075, -2.21991902200819e-06, 0.00105056836170180,
        0.000211281908728735, 1.50467048145063e-16, 4.65320742243753e-05
    };
    std::vector<double> A3_95 = {
        -8.38879451973883e-14, 3.94999999993448, 1.00000000003592, -3.59728437762498e-11,
        1.33992037944711e-13, 1.77967074930255e-13, -7.07762567264371e-11, 1.00000000052210,
        0.000524420077433023, -3.04128558605071e-07, -0.0778343474913371, -2.20301702096228,
        2.73311407491194e-06, -0.0149473707861847, -0.00464453996298439, -0.191357977879329
    };
    std::vector<double> B3_95 = {
        -4.96153817315486e-13, -4.64563782961434e-18, 1.29372776597539e-12,
        -2.58099696016061e-13, -3.72131332334172e-18, -2.00867616573371e-12,
        -0.00162960452959253, -2.21991902205642e-06, 0.00102446828994419,
        0.000206032680438813, 1.46515832626729e-16, 4.53760309965550e-05
    };

    std::vector<double> A3_90 = {
        -7.08096462873331e-14, 3.89999999993529, 1.00000000003501, -3.74440221893866e-11,
        8.83984657112187e-14, 1.62367028164840e-13, -6.88431018552365e-11, 1.00000000050422,
        0.000511228063124853, -3.04127720291887e-07, -0.0767309605286780, -2.17845549556317,
        2.66435789355151e-06, -0.0149473707864187, -0.00459167579237172, -0.188845238074267
    };
    std::vector<double> B3_90 = {
        -4.64285105954456e-13, -4.42912710640091e-18, 1.24174704806829e-12,
        -2.43987946654815e-13, -3.62183860334375e-18, -1.92973842099234e-12,
        -0.00158860983421589, -2.21991902206455e-06, 0.000998696521660729,
        0.000200849480198465, 1.42619355275743e-16, 4.42345293317750e-05
    };

    std::vector<double> A3_85 = {
        -3.10062965632281e-14, 3.84999999993609, 1.00000000003434, -4.42195295835289e-11,
        4.43721424406109e-14, 1.68512692914096e-13, -6.69393473683235e-11, 1.00000000048670,
        0.000498204099481443, -3.04126897475938e-07, -0.0756331524043924, -2.15375598382889,
        2.59647747698272e-06, -0.0149473707861329, -0.00453845269645144, -0.186342694018917
    };
    std::vector<double> B3_85 = {
        -4.34052866720778e-13, -4.21970299967415e-18, 1.19121157165711e-12,
        -2.30467010990168e-13, -3.52308126217456e-18, -1.85291509071587e-12,
        -0.00154813736410918, -2.21991902212669e-06, 0.000973253056864856,
        0.000195732307992946, 1.38772071843206e-16, 4.31075692234765e-05
    };

    std::vector<double> A3_80 = {
        2.04841254552793e-14, 3.79999999993675, 1.00000000003364, -4.60264085506388e-11,
        1.49327879357590e-15, 1.69925741381328e-13, -6.50675283065473e-11, 1.00000000046985,
        0.000485348186495237, -3.04126093770865e-07, -0.0745358808190186, -2.12902228302665,
        2.52947282684493e-06, -0.0149473707861846, -0.00448527792891425, -0.183834812219759
    };
    std::vector<double> B3_80 = {
        -4.05455883059858e-13, -4.00960377116235e-18, 1.14209179223265e-12,
        -2.17516921257025e-13, -3.42625773848026e-18, -1.77817630852260e-12,
        -0.00150818711921493, -2.21991902214750e-06, 0.000948137895485899,
        0.000190681163833172, 1.34980470849132e-16, 4.19951506705809e-05
    };

    std::vector<double> A3_75 = {
        1.01819641086105e-14, 3.74999999993772, 1.00000000003275, -4.54253795849307e-11,
        -3.95029792066643e-14, 1.67056271899159e-13, -6.32186205588797e-11, 1.00000000045314,
        0.000472660324193826, -3.04125215419307e-07, -0.0734433088016021, -2.10416528304654,
        2.46334394432031e-06, -0.0149473707861980, -0.00443180150662047, -0.181334705119356
    };
    std::vector<double> B3_75 = {
        -3.78425319657249e-13, -3.81550824891711e-18, 1.09436385734712e-12,
        -2.05123149092441e-13, -3.33109430846002e-18, -1.70549376657142e-12,
        -0.00146875909958376, -2.21991902218520e-06, 0.000923351037521288,
        0.000185696047713399, 1.31239703111035e-16, 4.08972736684625e-05
    };

    std::vector<double> A3_70 = {
        6.84668360779129e-14, 3.69999999993788, 1.00000000003180, -5.15362975568181e-11,
        -7.96503720351989e-14, 1.56817075809444e-13, -6.13977760375147e-11, 1.00000000043707,
        0.000460140512599520, -3.04124287562347e-07, -0.0723531792046976, -2.07922150411463,
        2.39809083064692e-06, -0.0149473707861614, -0.00437816568323339, -0.178836257387762
    };
    std::vector<double> B3_70 = {
        -3.52880990322738e-13, -3.61900858149558e-18, 1.04800608295708e-12,
        -1.93268762195699e-13, -3.23721542275853e-18, -1.63483878432675e-12,
        -0.00142985330522538, -2.21991902224124e-06, 0.000898892482943580,
        0.000180776959633020, 1.27553622330231e-16, 3.98139382135401e-05
    };

    std::vector<double> A3_65 = {
        8.45044925870338e-14, 3.64999999993856, 1.00000000003107, -5.20352622963858e-11,
        -1.18181582265110e-13, 1.58102614668155e-13, -5.96061953617280e-11, 1.00000000042138,
        0.000447788751718305, -3.04123429338922e-07, -0.0712657821149120, -2.05421947025212,
        2.33371348734260e-06, -0.0149473707860753, -0.00432448596534239, -0.176337455177220
    };
    std::vector<double> B3_65 = {
        -3.28699377453967e-13, -3.43148122218586e-18, 1.00300518758589e-12,
        -1.81938642114483e-13, -3.14426585584781e-18, -1.56618297771451e-12,
        -0.00139146973610854, -2.21991902225244e-06, 0.000874762231700232,
        0.000175923899591443, 1.23918605721970e-16, 3.87451443022532e-05
    };

    std::vector<double> A3_60 = {
        9.71245016353355e-14, 3.59999999994033, 1.00000000003033, -5.08980477588319e-11,
        -1.54206781051474e-13, 1.07801803378058e-13, -5.78478342803789e-11, 1.00000000040602,
        0.000435605041581371, -3.04112156005630e-07, -0.0701824167084689, -2.02909070901393,
        2.27021191577967e-06, -0.0149473707785411, -0.00427048935382220, -0.173845817918863
    };
    std::vector<double> B3_60 = {
        -3.05924353327592e-13, -3.25160536857479e-18, 9.59332439876004e-13,
        -1.71115340816327e-13, -3.05299041316327e-18, -1.49949714950903e-12,
        -0.00135360839229657, -2.21991902230393e-06, 0.000850960283798518,
        0.000171136867588089, 1.20337229309853e-16, 3.76908919310719e-05
    };

    std::vector<double> A3_55 = {
        1.33478689131089e-13, 3.54999999994097, 1.00000000002960, -5.80982708759753e-11,
        -1.89953981410919e-13, 1.07352296180381e-13, -5.61103837728662e-11, 1.00000000039118,
        0.000423589382191330, -3.04111582418524e-07, -0.0691022516295529, -2.00388536436178,
        2.20758611717361e-06, -0.0149473707785817, -0.00421637602957460, -0.171356000511093
    };
    std::vector<double> B3_55 = {
        -2.84446612980054e-13, -3.09050499418502e-18, 9.16966697885047e-13,
        -1.60785231781351e-13, -2.96255877588796e-18, -1.43475311045380e-12,
        -0.00131626927374410, -2.21991902230791e-06, 0.000827486639178234,
        0.000166415863623669, 1.16807404473474e-16, 3.66511810967027e-05
    };

    std::vector<double> A3_50 = {
        1.24205924808632e-13, 3.49999999994185, 1.00000000002883, -5.59444943305279e-11,
        -2.23873009836693e-13, 1.59202245580764e-13, -5.44040060702302e-11, 1.00000000037654,
        0.000411741773581632, -3.04120703202617e-07, -0.0680232545347252, -1.97863379459139,
        2.14583609310436e-06, -0.0149473707857273, -0.00416226391013546, -0.168862709469597
    };
    std::vector<double> B3_50 = {
        -2.64185892766037e-13, -2.91929984939772e-18, 8.75890614378261e-13,
        -1.50932530376500e-13, -2.87414521732783e-18, -1.37192158709018e-12,
        -0.00127945238051391, -2.21991902235432e-06, 0.000804341297847681,
        0.000161760887696341, 1.13332612606478e-16, 3.56260117952877e-05
    };

    std::vector<double> A3_45 = {
        1.65793196432262e-13, 3.44999999994253, 1.00000000002799, -5.83828409054535e-11,
        -2.56621152842878e-13, 8.71144559021800e-14, -5.27216357066812e-11, 1.00000000036241,
        0.000400062215755562, -3.04110460419934e-07, -0.0669479298579427, -1.95329136418706,
        2.08496184482896e-06, -0.0149473707784967, -0.00410797858600565, -0.166373076884703
    };
    std::vector<double> B3_45 = {
        -2.45134368170002e-13, -2.75757157892653e-18, 8.36080336907852e-13,
        -1.41540497379583e-13, -2.78682302689946e-18, -1.31097378876829e-12,
        -0.00124315771257367, -2.21991902237450e-06, 0.000781524259755848,
        0.000157171939805558, 1.09909499925642e-16, 3.46153840232924e-05
    };

    std::vector<double> A3_40 = {
        1.75160897267003e-13, 3.39999999994363, 1.00000000002721, -5.73576737441742e-11,
        -2.87615125795551e-13, 8.39421244355712e-14, -5.10673284073749e-11, 1.00000000034867,
        0.000388550708742301, -3.04109838963711e-07, -0.0658771512481716, -1.92780973946420,
        2.02496337391495e-06, -0.0149473707786119, -0.00405332725817195, -0.163892352185450
    };
    std::vector<double> B3_40 = {
        -2.27189807018032e-13, -2.59902687961548e-18, 7.97517976137974e-13,
        -1.32597145403132e-13, -2.70145517790651e-18, -1.25188106883385e-12,
        -0.00120738526996528, -2.21991902242839e-06, 0.000759035524898984,
        0.000152649019953355, 1.06539881842952e-16, 3.36192977777363e-05
    };

    std::vector<double> A3_35 = {
        2.24243171380183e-13, 3.34999999994405, 1.00000000002647, -6.30080060156786e-11,
        -3.17910114071461e-13, 1.05470235200054e-13, -4.94447311424573e-11, 1.00000000033528,
        0.000377207252547174, -3.04109265016710e-07, -0.0648083380620230, -1.90228614333664,
        1.96584068161310e-06, -0.0149473735841698, -0.00399869597653415, -0.161408868344375
    };
    std::vector<double> B3_35 = {
        -2.10347341874861e-13, -2.46047915865617e-18, 7.60177511216487e-13,
        -1.24085774378661e-13, -2.61661699344587e-18, -1.19461465357764e-12,
        -0.00117213505266265, -2.21991902245506e-06, 0.000736875093231406,
        0.000148192128135178, 1.03222776610986e-16, 3.26377530543961e-05
    };

    std::vector<double> A3_30 = {
        2.36132797773404e-13, 3.29999999994514, 1.00000000002573, -5.85998515721656e-11,
        -3.46169453586957e-13, 9.83382980064343e-14, -4.78452842988695e-11, 1.00000000032232,
        0.000366031847192015, -3.04108600408211e-07, -0.0637407432276166, -1.87671715465140,
        1.90759376959556e-06, -0.0149473735841731, -0.00394406936085231, -0.158921908470985
    };
    std::vector<double> B3_30 = {
        -1.94513129279092e-13, -2.30721310253654e-18, 7.24045321853599e-13,
        -1.15992819053900e-13, -2.53388290159955e-18, -1.13914512131025e-12,
        -0.00113740706068684, -2.21991902247849e-06, 0.000715042964737512,
        0.000143801264353434, 9.99606220897767e-17, 3.16707498498363e-05
    };

    std::vector<double> A3_25 = {
        2.34901317058384e-13, 3.24999999994611, 1.00000000002497, -5.54698417945056e-11,
        -3.73113096907824e-13, 9.67644316756506e-14, -4.62767249739067e-11, 1.00000000030958,
        0.000355024492689806, -3.04107999796581e-07, -0.0626780198060388, -1.85101078233263,
        1.85022263925205e-06, -0.0149473721815078, -0.00388908473143915, -0.156444141644023
    };
    std::vector<double> B3_25 = {
        -1.79654443364753e-13, -2.16820961153514e-18, 6.89088001069871e-13,
        -1.08304536793835e-13, -2.45207805555843e-18, -1.08544387937574e-12,
        -0.00110320129403838, -2.21991902248487e-06, 0.000693539139389219,
        0.000139476428608177, 9.67510756864602e-17, 3.07182881607794e-05
    };

    std::vector<double> A3_20 = {
        2.63322143198711e-13, 3.19999999994632, 1.00000000002414, -5.74885358649953e-11,
        -3.99548698593500e-13, 1.10337284016709e-13, -4.47370825041659e-11, 1.00000000029704,
        0.000344185189063968, -3.04107416648321e-07, -0.0616173202053234, -1.82525364020527,
        1.79372729189593e-06, -0.0149473707786435, -0.00383408483655709, -0.153964435830106
    };
    std::vector<double> B3_20 = {
        -1.65784885601757e-13, -2.04983901072478e-18, 6.55294842368849e-13,
        -1.01006211362404e-13, -2.37141499345651e-18, -1.03348225359127e-12,
        -0.00106951775274458, -2.21991902251482e-06, 0.000672363617176128,
        0.000135217620895621, 9.35949419817981e-17, 2.97803679830372e-05
    };

    std::vector<double> A3_15 = {
        2.73372939151543e-13, 3.14999999994789, 1.00000000002354, -5.96538486620966e-11,
        -4.23337838521050e-13, 8.90393658658432e-14, -4.32236231493813e-11, 1.00000000028526,
        0.000333513936328038, -3.04106779868259e-07, -0.0605618553480810, -1.79934909227681,
        1.73810772919627e-06, -0.0149473721814877, -0.00377868734083453, -0.151495256684123
    };
    std::vector<double> B3_15 = {
        -1.52753711093722e-13, -1.92597675196815e-18, 6.22643484080726e-13,
        -9.40852322424995e-14, -2.29222746010671e-18, -9.83230854031150e-13,
        -0.00103635643679870, -2.21991902254173e-06, 0.000651516398066776,
        0.000131024841217679, 9.04912570339310e-17, 2.88569893138597e-05
    };

    std::vector<double> A3_10 = {
        2.73630564265106e-13, 3.09999999994898, 1.00000000002272, -5.49726743817928e-11,
        -4.46291996492909e-13, 9.65356917237226e-14, -4.17377002564079e-11, 1.00000000027352,
        0.000323010734496475, -3.04106154309873e-07, -0.0595077135865780, -1.77340397971468,
        1.68336395265383e-06, -0.0149473735841239, -0.00372331418085484, -0.149022335508694
    };
    std::vector<double> B3_10 = {
        -1.40564340381603e-13, -1.78424619036518e-18, 5.91116240157563e-13,
        -8.75271744363098e-14, -2.21501133395099e-18, -9.34660551232325e-13,
        -0.00100371734621403, -2.21991902258036e-06, 0.000630997482043114,
        0.000126898089571492, 8.74421796400372e-17, 2.79481521488649e-05
    };

    std::vector<double> A3_05 = {
        2.99293429996257e-13, 3.04999999994943, 1.00000000002230, -9.22053390107744e-11,
        -4.68485824024097e-13, 1.17241574708443e-13, -4.07287678890015e-11, 1.00000000032097,
        0.000312675583604042, -3.04109531553727e-07, -0.0584157465285275, -1.74550395250290,
        1.62949596386932e-06, -0.0149473707794795, -0.00367408310474138, -0.146549081186278
    };
    std::vector<double> B3_05 = {
        -1.28523286361650e-13, -1.70816432501727e-18, 5.47323417057072e-13,
        -8.28419003722725e-14, -2.16373689381424e-18, -8.66644733068882e-13,
        -0.000971600480758073, -2.21991902256551e-06, 0.000610806866361585,
        0.000122837366129456, 8.54606765605554e-17, 2.70538539415640e-05
    };

    std::vector<double> A3_00 = {
        3.05149366265809e-13, 2.99999999995042, 1.00000000002202, -8.98866598616259e-11,
        -4.88911138505945e-13, 1.11939952492970e-13, -3.98812210801150e-11, 1.00000000030044,
        0.000302508483661247, -3.04108159029544e-07, -0.0573145567519485, -1.71689842575197,
        1.57650376403890e-06, -0.0149473707793941, -0.00362671168783090, -0.144079843752658
    };
    std::vector<double> B3_00 = {
        -1.16907095927519e-13, -1.61094923018803e-18, 5.01453262594098e-13,
        -7.86090098343763e-14, -2.12127340733339e-18, -7.95064248744720e-13,
        -0.000940005840695205, -2.21991902259168e-06, 0.000590944553226427,
        0.000118842670749934, 8.38150316060650e-17, 2.61740966119779e-05
    };
    std::vector<double> A2_95 = {
        3.27903505011285e-13, 2.94999999995087, 1.00000000002152, -8.22004022534317e-11,
        -5.08598272286772e-13, 1.11593213437539e-13, -3.90171048528620e-11, 1.00000000028069,
        0.000292509434668889, -3.04106908425986e-07, -0.0562191179925915, -1.68825465727950,
        1.52438735510081e-06, -0.0149473707793319, -0.00357876091730501, -0.141617329621690
    };
    std::vector<double> B2_95 = {
        -1.06119122890711e-13, -1.52197969532159e-18, 4.58444588714857e-13,
        -7.45255940978773e-14, -2.07796940249752e-18, -7.27882770739817e-13,
        -0.000908933426032108, -2.21991902259905e-06, 0.000571410543626081,
        0.000114914003381515, 8.21283024938920e-17, 2.53088811143889e-05
    };

    std::vector<double> A2_90 = {
        3.09311894575587e-13, 2.89999999995217, 1.00000000002132, -7.97235134314035e-11,
        -4.98855108030089e-13, 9.49811129304789e-14, -3.81417023652614e-11, 1.00000000026210,
        0.000282678433379412, -3.04105393013320e-07, -0.0551266428800224, -1.65965918351381,
        1.47314345210058e-06, -0.0149473721816993, -0.00353057335858804, -0.139150494314873
    };
    std::vector<double> B2_90 = {
        -9.61097030177016e-14, -1.44506481889794e-18, 4.18189246252643e-13,
        -7.05889389647638e-14, -2.03346516198521e-18, -6.64917860862828e-13,
        -0.000878383236784850, -2.21991902261979e-06, 0.000552204837501107,
        0.000111051364022283, 8.04037771585227e-17, 2.44582074262866e-05
    };

    std::vector<double> A2_85 = {
        2.99295617765142e-13, 2.84999999995236, 1.00000000002073, -6.96655026231128e-11,
        -4.72526632300071e-13, 7.34080418346467e-14, -3.72460604130607e-11, 1.00000000024416,
        0.000273015481353273, -3.04103643311372e-07, -0.0540395522985808, -1.63105129852794,
        1.42277337174304e-06, -0.0149473721816455, -0.00348190969720583, -0.136687712010372
    };
    std::vector<double> B2_85 = {
        -8.68610035705899e-14, -1.36077833746836e-18, 3.80563365109249e-13,
        -6.67979158255946e-14, -1.98828268908207e-18, -6.05992370643165e-13,
        -0.000848355272955952, -2.21991902265533e-06, 0.000533327434786657,
        0.000107254752678907, 7.86429055392872e-17, 2.36220755272316e-05
    };

    std::vector<double> A2_80 = {
        2.73928365701321e-13, 2.79999999995363, 1.00000000002029, -6.66951138888106e-11,
        -4.46213328029151e-13, 8.85056064913600e-14, -3.63528200036429e-11, 1.00000000022683,
        0.000263520580880157, -3.04102210628128e-07, -0.0529573515481892, -1.60244627824872,
        1.37327935745499e-06, -0.0149473735839246, -0.00343283039878630, -0.134227029489206
    };
    std::vector<double> B2_80 = {
        -7.82730521329852e-14, -1.28334248433810e-18, 3.45451867388462e-13,
        -6.31478873054000e-14, -1.94230792014395e-18, -5.50933335952998e-13,
        -0.000818849534527629, -2.21991902265229e-06, 0.000514778335408130,
        0.000103524169344853, 7.68487851655483e-17, 2.28004853948522e-05
    };

    std::vector<double> A2_75 = {
        2.60994591626722e-13, 2.74999999995419, 1.00000000001997, -6.80595548867290e-11,
        -4.21418067087981e-13, 5.84755685703314e-14, -3.54417337181303e-11, 1.00000000021080,
        0.000254193731955906, -3.04100805750921e-07, -0.0518804275305642, -1.57382153655478,
        1.32466140700349e-06, -0.0149473721810988, -0.00338324532353888, -0.131770871061938
    };
    std::vector<double> B2_75 = {
        -7.03790618106306e-14, -1.19633903380558e-18, 3.12737743958947e-13,
        -5.96361435054710e-14, -1.89583414481668e-18, -4.99570086006203e-13,
        -0.000789866021522563, -2.21991902266336e-06, 0.000496557539319090,
        9.98596140287675e-05, 7.50247230359885e-17, 2.19934370099219e-05
    };

    std::vector<double> A2_70 = {
        2.18987049617074e-13, 2.69999999995529, 1.00000000001954, -4.91318201178844e-11,
        -3.97162577779128e-13, 4.85116362841178e-14, -3.45168672142503e-11, 1.00000000019529,
        0.000245034934562529, -3.04099412869708e-07, -0.0508074306397321, -1.54522346942990,
        1.27691951833524e-06, -0.0149473721808688, -0.00333333844602999, -0.129313510170114
    };
    std::vector<double> B2_70 = {
        -6.30793269764717e-14, -1.13330262491934e-18, 2.82314861592414e-13,
        -5.62614745362399e-14, -1.84823580059915e-18, -4.51736867042914e-13,
        -0.000761404733943123, -2.21991902269018e-06, 0.000478665046463169,
        9.62610867252765e-05, 7.31734565143581e-17, 2.12009303513632e-05
    };

    std::vector<double> A2_65 = {
        2.36347242216405e-13, 2.64999999995579, 1.00000000001903, -5.51623758163527e-11,
        -3.74723360413025e-13, 6.21773202782399e-14, -3.35935398628541e-11, 1.00000000018051,
        0.000236044188688285, -3.04098310742857e-07, -0.0497399487125223, -1.51659489469514,
        1.23005368840810e-06, -0.0149473707778186, -0.00328288285877862, -0.126861903051396
    };
    std::vector<double> B2_65 = {
        -5.64208180316709e-14, -1.05767871465958e-18, 2.54072208150796e-13,
        -5.30183460840711e-14, -1.80023638393875e-18, -4.07270867683966e-13,
        -0.000733465671784876, -2.21991902271409e-06, 0.000461100856782390,
        9.27285874371294e-05, 7.12980169054505e-17, 2.04229653997269e-05
    };

    std::vector<double> A2_60 = {
        2.08972607399084e-13, 2.59999999995707, 1.00000000001837, -4.70333973149745e-11,
        -3.52271900233835e-13, 1.26554833118697e-14, -3.26530597599557e-11, 1.00000000016665,
        0.000227221494318109, -3.04097117068239e-07, -0.0486768556256378, -1.48798890405930,
        1.18406391534412e-06, -0.0149473721803908, -0.00323208998356500, -0.124410059507814
    };
    std::vector<double> B2_60 = {
        -5.02688785270684e-14, -9.73419719595673e-19, 2.27908241324376e-13,
        -4.99061285153458e-14, -1.75165533071907e-18, -3.66013638769872e-13,
        -0.000706048835050094, -2.21991902273474e-06, 0.000443864970225729,
        8.92621161640491e-05, 6.94005402113028e-17, 1.96595421360004e-05
    };

    std::vector<double> A2_55 = {
        2.07238978796908e-13, 2.54999999995724, 1.00000000001797, -4.32577493168313e-11,
        -3.31535697373319e-13, 6.09855802377692e-14, -3.17187862084861e-11, 1.00000000015332,
        0.000218566851442594, -3.04098121613368e-07, -0.0476193156294486, -1.45936128752666,
        1.13895019656679e-06, -0.0149473721827156, -0.00318078425253205, -0.121963272430573
    };
    std::vector<double> B2_55 = {
        -4.46577158373456e-14, -9.09142595081378e-19, 2.03720775058397e-13,
        -4.69221963057836e-14, -1.70278665344203e-18, -3.27810265386117e-13,
        -0.000679154223727487, -2.21991902270763e-06, 0.000426957386736209,
        8.58616729081957e-05, 6.74840200254806e-17, 1.89106605414147e-05
    };

    std::vector<double> A2_50 = {
        1.85948903624946e-13, 2.49999999995830, 1.00000000001762, -4.05860426948557e-11,
        -3.11206288268641e-13, 5.60805316838911e-14, -3.07759690948322e-11, 1.00000000014077,
        0.000210080260060588, -3.04097138943220e-07, -0.0465659728293839, -1.43074663245396,
        1.09471252971974e-06, -0.0149473721825508, -0.00312910206602502, -0.119516802103210
    };
    std::vector<double> B2_50 = {
        -3.95311499359261e-14, -8.48001604627045e-19, 1.81412773258421e-13,
        -4.40612725791688e-14, -1.65347419793285e-18, -2.92509470624694e-13,
        -0.000652781837852994, -2.21991902274173e-06, 0.000410378106289095,
        8.25272576694322e-05, 6.55516870162711e-17, 1.81763205976742e-05
    };

    std::vector<double> A2_45 = {
        1.75756195149775e-13, 2.44999999995930, 1.00000000001696, -3.49433061170217e-11,
        -2.91494775388937e-13, 4.31622596803590e-14, -2.98261612584921e-11, 1.00000000012903,
        0.000201761720143844, -3.04096153842942e-07, -0.0455167965171715, -1.40214639147396,
        1.05135091246084e-06, -0.0149473721823008, -0.00307704920863906, -0.117070484770529
    };
    std::vector<double> B2_45 = {
        -3.48624282993734e-14, -7.82345909840705e-19, 1.60876899909373e-13,
        -4.13226303124468e-14, -1.60402134264426e-18, -2.59964262505521e-13,
        -0.000626931677401734, -2.21991902274915e-06, 0.000394127128823889,
        7.92588704476940e-05, 6.36063924284446e-17, 1.74565222868442e-05
    };

    std::vector<double> A2_40 = {
        1.73305206180524e-13, 2.39999999996024, 1.00000000001645, -3.38791285436148e-11,
        -2.72916647527139e-13, 3.97243619996734e-14, -2.88792776116198e-11, 1.00000000011782,
        0.000193611231692934, -3.04095299227230e-07, -0.0444732036795088, -1.37352713084624,
        1.00886534232291e-06, -0.0149473721820811, -0.00302449405740272, -0.114629045104024
    };
    std::vector<double> B2_40 = {
        -3.06148033881485e-14, -7.08303781626188e-19, 1.42027096463812e-13,
        -3.87036263651907e-14, -1.55456666502240e-18, -2.30031295525351e-13,
        -0.000601603742385905, -2.21991902276382e-06, 0.000378204454305799,
        7.60565112441804e-05, 6.16502111587375e-17, 1.67512655917416e-05
    };

    std::vector<double> A2_35 = {
        1.54363847813419e-13, 2.34999999996107, 1.00000000001598, -2.79513505245852e-11,
        -2.55004731579788e-13, 3.92838048562457e-14, -2.79262239334557e-11, 1.00000000010734,
        0.000185628794693582, -3.04094464278993e-07, -0.0434348110770961, -1.34489736445509,
        9.67255817020432e-07, -0.0149473721818902, -0.00297147008183728, -0.112191249820603
    };
    std::vector<double> B2_35 = {
        -2.67513005581351e-14, -6.57270129936102e-19, 1.24772357701884e-13,
        -3.62003954903788e-14, -1.50454437218164e-18, -2.02571080298430e-13,
        -0.000576798032797407, -2.21991902274782e-06, 0.000362610082689644,
        7.29201800594318e-05, 5.96862887311921e-17, 1.60605504955045e-05
    };

    std::vector<double> A2_30 = {
        1.64434911513247e-13, 2.29999999996166, 1.00000000001541, -2.62492948166055e-11,
        -2.38323094531263e-13, 3.37056147201143e-14, -2.69784814668815e-11, 1.00000000009726,
        0.000177814409134776, -3.04093821616920e-07, -0.0424005896958447, -1.31628314581384,
        9.26522334065719e-07, -0.0149473721820062, -0.00291808001137052, -0.109753518338592
    };
    std::vector<double> B2_30 = {
        -2.32817789783106e-14, -6.01063301334541e-19, 1.09017748170862e-13,
        -3.38106806315676e-14, -1.45449057914188e-18, -1.77447717502124e-13,
        -0.000552514548655085, -2.21991902276485e-06, 0.000347344013949441,
        6.98498768952533e-05, 5.77174685733231e-17, 1.53843769820193e-05
    };

    std::vector<double> A2_25 = {
        1.42428305405216e-13, 2.24999999996228, 1.00000000001490, -2.35625819279581e-11,
        -2.22482951456457e-13, 2.72772917707596e-14, -2.60313526798241e-11, 1.00000000008802,
        0.000170168075015595, -3.04093132249071e-07, -0.0413712068622961, -1.28766807908978,
        8.86664891159054e-07, -0.0149473721817521, -0.00286425923907816, -0.107318130537131
    };
    std::vector<double> B2_25 = {
        -2.01520453605847e-14, -5.45252144850987e-19, 9.46800553540059e-14,
        -3.15319068393359e-14, -1.40465016474775e-18, -1.54529374763086e-13,
        -0.000528753289957488, -2.21991902280529e-06, 0.000332406248048675,
        6.68456017489424e-05, 5.57464041174891e-17, 1.47227450346519e-05
    };

    std::vector<double> A2_20 = {
        1.10661569213953e-13, 2.19999999996317, 1.00000000001441, -1.87765694268018e-11,
        -2.06874203779360e-13, 1.70833252749846e-14, -2.50853838527138e-11, 1.00000000007924,
        0.000162689792313891, -3.04092536638412e-07, -0.0403466727132332, -1.25905316420605,
        8.47683486089927e-07, -0.0149473721815235, -0.00281001181387304, -0.104885015780602
    };
    std::vector<double> B2_20 = {
        -1.73352375577272e-14, -5.03971665645716e-19, 8.16780318510166e-14,
        -2.93610221953055e-14, -1.35458282490435e-18, -1.33688035392874e-13,
        -0.000505514256689637, -2.21991902279844e-06, 0.000317796784944467,
        6.39073546229101e-05, 5.37757050589095e-17, 1.40756546383629e-05
    };

    std::vector<double> A2_15 = {
        1.06082898084662e-13, 2.14999999996430, 1.00000000001382, -1.70017256299665e-11,
        -1.91994696267945e-13, 2.65621429499670e-14, -2.41449722520313e-11, 1.00000000007110,
        0.000155379561024426, -3.04092002389481e-07, -0.0393268440143640, -1.23043931850272,
        8.09578116125987e-07, -0.0149473721813218, -0.00275534112350788, -0.102453904010533
    };
    std::vector<double> B2_15 = {
        -1.48238403704978e-14, -4.37745028118276e-19, 6.99227053211141e-14,
        -2.72950058334344e-14, -1.30471247657784e-18, -1.14799343738262e-13,
        -0.000482797448866847, -2.21991902279954e-06, 0.000303515624615062,
        6.10351355176253e-05, 5.18082762473946e-17, 1.34431057781544e-05
    };

    std::vector<double> A2_10 = {
        1.27177213203931e-13, 2.09999999996511, 1.00000000001319, -1.49371614251011e-11,
        -1.77880349109302e-13, 1.85996537360634e-14, -2.32097839602341e-11, 1.00000000006348,
        0.000148237381139325, -3.04091546322648e-07, -0.0383116380508542, -1.20182481966393,
        7.72348779103959e-07, -0.0149473721813131, -0.00270024005375505, -0.100024829093278
    };
    std::vector<double> B2_10 = {
        -1.25702827602701e-14, -3.94720832407028e-19, 5.93429249967289e-14,
        -2.53315082450900e-14, -1.25520570289774e-18, -9.77429074330361e-14,
        -0.000460602866490894, -2.21991902280983e-06, 0.000289562767032215,
        5.82289444341799e-05, 4.98467354497396e-17, 1.28250984396351e-05
    };

    std::vector<double> A2_05 = {
        8.37339359278806e-14, 2.04999999996589, 1.00000000001276, -1.11690545839321e-11,
        -1.64776265009528e-13, 1.81195417749229e-14, -2.22788892286466e-11, 1.00000000005627,
        0.000141263252653640, -3.04091111718952e-07, -0.0373011457208380, -1.17321826065141,
        7.35995472600053e-07, -0.0149473721810962, -0.00264474340873531, -0.0975971890641141
    };
    std::vector<double> B2_05 = {
        -1.05841561364818e-14, -3.53030657784654e-19, 4.98576467654679e-14,
        -2.34676828895511e-14, -1.20607352522389e-18, -8.24022409181032e-14,
        -0.000438930509566870, -2.21991902284776e-06, 0.000275938212171716,
        5.54887813714048e-05, 4.78936410608404e-17, 1.22216326083833e-05
    };

    std::vector<double> A2_00 = {
        9.05649584650802e-14, 1.99999999996664, 1.00000000001223, -9.46672150067366e-12,
        -1.52771014498061e-13, 1.59098002320443e-14, -2.13546326904026e-11, 1.00000000004975,
        0.000134457175550217, -3.04090725576200e-07, -0.0362957126237417, -1.14460153339465,
        7.00518194354120e-07, -0.0149473721811751, -0.00258877899584769, -0.0951729752223792
    };
    std::vector<double> B2_00 = {
        -8.82092113070538e-15, -3.07813370562036e-19, 4.13947034202250e-14,
        -2.17012272259813e-14, -1.15658929146909e-18, -6.86645644696887e-14,
        -0.000417780378079634, -2.21991902282821e-06, 0.000262641959998515,
        5.28146463316389e-05, 4.59516732654799e-17, 1.16327082712193e-05
    };
    std::vector<double> A1_95 = {
        9.15840728514946e-14, 1.94999999996758, 1.00000000001160, -7.60888149707488e-12,
        -1.40355184761343e-13, 1.12745767315957e-14, -2.04431905166481e-11, 1.00000000004369,
        0.000127819149825327, -3.04090391495434e-07, -0.0352951356276948, -1.11598302327911,
        6.65916941849231e-07, -0.0149473721809319, -0.00253238015317309, -0.0927512074527795
    };
    std::vector<double> B1_95 = {
        -7.27317054969638e-15, -2.73069215278863e-19, 3.38812501275362e-14,
        -2.00279573269219e-14, -1.10757627992792e-18, -5.64207473554853e-14,
        -0.000397152472050992, -2.21991902285496e-06, 0.000249674010502622,
        5.02065393130351e-05, 4.40235377727568e-17, 1.10583254145083e-05
    };

    std::vector<double> A1_90 = {
        8.71160440225692e-14, 1.89999999996834, 1.00000000001108, -8.38238796015633e-12,
        -1.29302200823931e-13, 1.65451281073330e-14, -1.95447918639052e-11, 1.00000000003804,
        0.000121349175472679, -3.04090140910133e-07, -0.0342988053367930, -1.08737835025607,
        6.32191712651915e-07, -0.0149473721809124, -0.00247560848643963, -0.0903297492186989
    };
    std::vector<double> B1_90 = {
        -5.92183729637818e-15, -2.39170663608892e-19, 2.72500663527075e-14,
        -1.84466698999658e-14, -1.05951821472826e-18, -4.55657392920599e-14,
        -0.000377046791472444, -2.21991902287621e-06, 0.000237034363656723,
        4.76644603178976e-05, 4.21117418865086e-17, 1.04984840259971e-05
    };

    std::vector<double> A1_85 = {
        6.74130516917269e-14, 1.84999999996918, 1.00000000001060, -5.21144545272349e-12,
        -1.18666035593196e-13, 1.74746151479040e-14, -1.86438243442255e-11, 1.00000000003292,
        0.000115047252481678, -3.04089840543751e-07, -0.0333080898684701, -1.05875329987668,
        5.99342504711769e-07, -0.0149473721807309, -0.00241832912331823, -0.0879133252389154
    };
    std::vector<double> B1_85 = {
        -4.74854401055139e-15, -2.08954323065007e-19, 2.14347447126742e-14,
        -1.69541831166366e-14, -1.01212976312916e-18, -3.59981718099731e-14,
        -0.000357463336342166, -2.21991902287018e-06, 0.000224723019439405,
        4.51884093453062e-05, 4.02189497168375e-17, 9.95318409318907e-06
    };

    std::vector<double> A1_80 = {
        6.14011469904079e-14, 1.79999999996998, 1.00000000001008, -4.64229355086277e-12,
        -1.09100679469877e-13, 9.49840751193476e-15, -1.77701934846157e-11, 1.00000000002812,
        0.000108913380845800, -3.04089352281514e-07, -0.0323217510997670, -1.03013832428725,
        5.67369315419341e-07, -0.0149473721798776, -0.00236066206125906, -0.0854977028164335
    };
    std::vector<double> B1_80 = {
        -3.72900879825884e-15, -1.71870559212470e-19, 1.63714189232132e-14,
        -1.55474111319355e-14, -9.64316896529636e-19, -2.76208465550386e-14,
        -0.000338402106668359, -2.21991902287636e-06, 0.000212739977837172,
        4.27783863965677e-05, 3.83475385367806e-17, 9.42242560454389e-06
    };

    std::vector<double> A1_75 = {
        5.96432032495061e-14, 1.74999999997098, 1.00000000000954, -2.72902776346390e-12,
        -9.91407390299487e-14, -5.69390185874216e-16, -1.68960519731232e-11, 1.00000000002376,
        0.000102947560557056, -3.04089166950811e-07, -0.0313401765188081, -1.00152366523347,
        5.36272142528513e-07, -0.0149473721800334, -0.00230256882764496, -0.0830842258322976
    };
    std::vector<double> B1_75 = {
        -2.86071486860739e-15, -1.39547893358732e-19, 1.19980861725856e-14,
        -1.42239119192161e-14, -9.18301374693662e-19, -2.03396581668165e-14,
        -0.000319863102452450, -2.21991902289697e-06, 0.000201085238833825,
        4.04343914720083e-05, 3.65002378620315e-17, 8.90620854877053e-06
    };

    std::vector<double> A1_70 = {
        5.37937424990081e-14, 1.64999999997253, 1.00000000000851, -1.09069807931791e-12,
        -8.24971210836142e-14, -8.64807369402798e-15, -1.51974294397220e-11, 1.00000000001616,
        9.15200739950942e-05, -3.04088902002311e-07, -0.0293912460471435, -0.944298399633515,
        4.76705836449900e-07, -0.0149473714785268, -0.00218511618417612, -0.0782633472351850
    };
    std::vector<double> B1_70 = {
        -1.49962567359814e-15, -9.40862605277232e-20, 5.09436107664837e-15,
        -1.18157591956638e-14, -8.27134232915892e-19, -8.71066660099694e-15,
        -0.000284351770388096, -2.21991902290101e-06, 0.000178760668558853,
        3.59444856944480e-05, 3.28875384020817e-17, 7.91739869242931e-06
    };

    std::vector<double> A1_65 = {
        5.37937424990081e-14, 1.64999999997253, 1.00000000000851, -1.09069807931791e-12,
        -8.24971210836142e-14, -8.64807369402798e-15, -1.51974294397220e-11, 1.00000000001616,
        9.15200739950942e-05, -3.04088902002311e-07, -0.0293912460471435, -0.944298399633515,
        4.76705836449900e-07, -0.0149473714785268, -0.00218511618417612, -0.0782633472351850
    };
    std::vector<double> B1_65 = {
        -1.49962567359814e-15, -9.40862605277232e-20, 5.09436107664837e-15,
        -1.18157591956638e-14, -8.27134232915892e-19, -8.71066660099694e-15,
        -0.000284351770388096, -2.21991902290101e-06, 0.000178760668558853,
        3.59444856944480e-05, 3.28875384020817e-17, 7.91739869242931e-06
    };

    std::vector<double> A1_60 = {
        3.88426477024941e-14, 1.59999999997338, 1.00000000000801, 1.08810448872257e-13,
        -7.48444538434014e-14, 1.23377426946556e-15, -1.43770009195517e-11, 1.00000000001297,
        8.60584077089651e-05, -3.04088816756785e-07, -0.0284241323530611, -0.915681898460560,
        4.48236698564514e-07, -0.0149473721798300, -0.00212573355167943, -0.0758567683475418
    };
    std::vector<double> B1_60 = {
        -9.84293292107709e-16, -6.52962437706438e-20, 2.45232978889065e-15,
        -1.07257722821632e-14, -7.82680486349141e-19, -4.19451562513277e-15,
        -0.000267379442545625, -2.21991902290757e-06, 0.000168090837263961,
        3.37985748419248e-05, 3.11269762637431e-17, 7.44480587168822e-06
    };

    std::vector<double> A1_55 = {
        3.13706925608613e-14, 1.54999999997408, 1.00000000000761, 3.63846098946323e-13,
        -6.83461295128662e-14, -1.37768638230678e-14, -1.35688600737595e-11, 1.00000000001007,
        8.07647927460569e-05, -3.04088736915713e-07, -0.0274617729648659, -0.887063853963513,
        4.20643567714302e-07, -0.0149473714785851, -0.00206591723512642, -0.0734524783523806
    };
    std::vector<double> B1_55 = {
        -5.63614331234141e-16, -5.86936326050528e-20, 2.82637953839815e-16,
        -9.70782229915911e-15, -7.38909460699029e-19, -4.37721216152605e-16,
        -0.000250929340164623, -2.21991902291807e-06, 0.000157749308516027,
        3.17186920154431e-05, 2.94001347666399e-17, 6.98675444328391e-06
    };
    
    std::vector<double> A1_50 = {
        3.14632311223029e-14, 1.49999999997496, 1.00000000000702, 5.31511703253300e-13,
        -6.09799274533346e-14, 1.71425383341976e-14, -1.26650601812356e-11, 1.00000000000855,
        7.56392290959861e-05, -3.04088733165472e-07, -0.0265663361672598, -0.857446258266571,
        3.93926441449104e-07, -0.0149473721796996, -0.00200568475095520, -0.0710501350586491
    };
    std::vector<double> B1_50 = {
        -3.05885976026752e-16, -3.70446557991115e-20, 2.43098455861984e-17,
        -8.60758637905969e-15, -6.89283439291824e-19, 5.54806580451567e-18,
        -0.000235001463249952, -2.21991902293386e-06, 0.000147736082367683,
        2.97048371941878e-05, 2.74446817303784e-17, 6.54324471362276e-06
    };

    std::vector<double> A1_45 = {
        2.69196425809927e-14, 1.44999999997584, 1.00000000000644, 9.53835562533040e-13,
        -5.46527286912673e-14, 6.39048592940548e-15, -1.17712771836577e-11, 1.00000000000757,
        7.06817167559298e-05, -3.04088704409907e-07, -0.0256807615830026, -0.827673668865574,
        3.68085317583577e-07, -0.0149473718292888, -0.00194502511874101, -0.0686498812605787
    };
    std::vector<double> B1_45 = {
        -1.19000100120535e-16, -2.47315156033385e-20, 2.24717840906566e-17,
        -7.58385636844918e-15, -6.41808816571545e-19, 1.27901975035399e-18,
        -0.000219595811794553, -2.21991902292178e-06, 0.000138051158725847,
        2.77570104010535e-05, 2.55266120986784e-17, 6.11427633890188e-06
    };

    std::vector<double> A1_40 = {
        2.02110976564499e-14, 1.39999999997673, 1.00000000000592, 6.83947230799185e-13,
        -4.85790987297279e-14, -1.52037790634013e-14, -1.09165471495194e-11, 1.00000000000668,
        6.58922557193833e-05, -3.04088677596717e-07, -0.0247952388031277, -0.797901002242620,
        3.43120193748512e-07, -0.0149473714786064, -0.00188393788072617, -0.0662519031294182
    };
    std::vector<double> B1_40 = {
        3.12783356433910e-17, -5.65382682073561e-21, 2.04981973235541e-17,
        -6.65612631267408e-15, -5.94993910454603e-19, -2.54321642413104e-18,
        -0.000204712385802577, -2.21991902291193e-06, 0.000128694537580558,
        2.58752116367492e-05, 2.36891232433764e-17, 5.69984926488434e-06
    };

    std::vector<double> A1_35 = {
        3.87468228001021e-14, 1.34999999997747, 1.00000000000540, 1.55992719845448e-12,
        -4.40474846860519e-14, -6.54669218400177e-15, -1.01044836667703e-11, 1.00000000000569,
        6.12708459834240e-05, -3.04088676156276e-07, -0.0239096404061216, -0.768130577554475,
        3.19031067651548e-07, -0.0149473718291645, -0.00182243210467539, -0.0638558337349560
    };
    std::vector<double> B1_35 = {
        1.40255617071960e-16, 1.27187494468340e-20, 1.86347248395946e-17,
        -5.81802319843342e-15, -5.51362519451931e-19, -5.79899931513975e-18,
        -0.000190351185277104, -2.21991902292444e-06, 0.000119666218932870,
        2.40594409009756e-05, 2.19317020331092e-17, 5.29996349162653e-06
    };

    std::vector<double> A1_30 = {
        3.67603924683566e-14, 1.29999999997830, 1.00000000000489, 1.21172789678897e-12,
        -3.86643488651810e-14, -4.53143048256564e-16, -9.33474070108497e-12, 1.00000000000485,
        5.68174875409009e-05, -3.04088700210543e-07, -0.0230240710886659, -0.738358989136828,
        2.95817936744363e-07, -0.0149473718292963, -0.00176049429818721, -0.0614621010073135
    };
    std::vector<double> B1_30 = {
        2.17545633890072e-16, 1.95542862462201e-20, 1.71270061934820e-17,
        -5.06333325900066e-15, -5.09205747917302e-19, -8.53438633917911e-18,
        -0.000176512210217772, -2.21991902294384e-06, 0.000110966202781739,
        2.23096981946750e-05, 2.02538049897597e-17, 4.91461901929573e-06
    };

    std::vector<double> A1_25 = {
        3.06803016880095e-14, 1.24999999997912, 1.00000000000454, 2.70995286446666e-13,
        -3.53590490452972e-14, 4.36963516882542e-15, -8.58836652890134e-12, 1.00000000000425,
        5.25321803972969e-05, -3.04088668932929e-07, -0.0221386052035571, -0.708584415306682,
        2.73480798881457e-07, -0.0149473718290189, -0.00169811728343699, -0.0590709585357931
    };
    std::vector<double> B1_25 = {
        2.69021514123591e-16, 2.24857203623222e-20, 1.59919820441612e-17,
        -4.38632798364771e-15, -4.68749912377429e-19, -1.08430805160391e-17,
        -0.000163195460622485, -2.21991902292681e-06, 0.000102594489125091,
        2.06259835167868e-05, 1.86547525634244e-17, 4.54381584774419e-06
    };

    std::vector<double> A1_20 = {
        1.94438974384020e-14, 1.19999999998004, 1.00000000000405, 1.13182622735720e-12,
        -2.95104129397437e-14, -3.01028059375345e-14, -7.88375142462160e-12, 1.00000000000354,
        4.84149245303581e-05, -3.04088662612035e-07, -0.0212530461767396, -0.678812621881867,
        2.52019652230749e-07, -0.0149473718289740, -0.00163532386202727, -0.0566816558306366
    };
    std::vector<double> B1_20 = {
        2.99413709918186e-16, 3.21089612051495e-20, 1.49797776721923e-17,
        -3.78144479492794e-15, -4.30601148167830e-19, -1.27309051972321e-17,
        -0.000150400936498987, -2.21991902296179e-06, 9.45510779670986e-05,
        1.90082968672435e-05, 1.71339200573535e-17, 4.18755397700924e-06
    };

    std::vector<double> A1_15 = {
        8.17655268025388e-15, 1.14999999998088, 1.00000000000367, 1.02302873935295e-12,
        -2.60054955212335e-14, -2.17090536844577e-14, -7.21812518966502e-12, 1.00000000000299,
        4.44657199511634e-05, -3.04088673914618e-07, -0.0203674838662985, -0.649040293433090,
        2.31434493781577e-07, -0.0149473720043479, -0.00157210087309447, -0.0542945932512460
    };
    std::vector<double> B1_15 = {
        3.12835505668306e-16, 3.36741628489864e-20, 1.40236892329790e-17,
        -3.24282345642850e-15, -3.94251577996147e-19, -1.41962721959821e-17,
        -0.000138128637841684, -2.21991902295666e-06, 8.68359693036053e-05,
        1.74566382458948e-05, 1.56905628258762e-17, 3.84583340712374e-06
    };

    std::vector<double> A1_10 = {
        1.75623483643728e-14, 1.09999999998162, 1.00000000000329, 6.15630127150444e-13,
        -2.27394175443636e-14, -1.25977144951247e-14, -6.58764339944501e-12, 1.00000000000240,
        4.06845666448409e-05, -3.04088689634945e-07, -0.0194819750852631, -0.619267176932321,
        2.11725321599843e-07, -0.0149473720044743, -0.00150844741935143, -0.0519098685685635
    };
    std::vector<double> B1_10 = {
        3.11068550495063e-16, 4.01048999156157e-20, 1.28878709902702e-17,
        -2.76559067132383e-15, -3.59530253823372e-19, -1.53948238163469e-17,
        -0.000126378564654005, -2.21991902294364e-06, 7.94491631361823e-05,
        1.59710076527135e-05, 1.43238034730165e-17, 3.51865413809732e-06
    };

    std::vector<double> A1_05 = {
        1.32134189701959e-14, 1.04999999998244, 1.00000000000294, 5.26575242860989e-13,
        -2.05101900854871e-14, 2.61826096508779e-15, -5.99601215396308e-12, 1.00000000000203,
        3.70714646182181e-05, -3.04088709251698e-07, -0.0185964381076265, -0.589494664496493,
        1.92892133376922e-07, -0.0149473720043920, -0.00144436893266569, -0.0495272544748810
    };
    std::vector<double> B1_05 = {
        2.99745677625605e-16, 3.67586607954055e-20, 1.21888041109894e-17,
        -2.34444391378034e-15, -3.26934869822356e-19, -1.62270336870243e-17,
        -0.000115150716937624, -2.21991902293651e-06, 7.23906594653495e-05,
        1.45514050870865e-05, 1.30328994270395e-17, 3.20601616985688e-06
    };
    std::vector<double> A1_00 = {
        2.23054448126036e-14, 0.999999999983276, 1.00000000000266, 2.86717885668046e-13,
        -1.86712384030688e-14, 9.46669152879898e-15, -5.44042692945485e-12, 1.00000000000172,
        3.36264138607990e-05, -3.04088708669801e-07, -0.0177108910144620, -0.559722309422723,
        1.74934927132404e-07, -0.0149473720044049, -0.00137986365271151, -0.0471468129506809
    };
    std::vector<double> B1_00 = {
        2.80553677930112e-16, 3.80576165732420e-20, 1.14900019295046e-17,
        -1.97494550194419e-15, -2.96659755119512e-19, -1.67966633440528e-17,
        -0.000104445094693318, -2.21991902293735e-06, 6.56604582911099e-05,
        1.31978305491429e-05, 1.18166620493253e-17, 2.90791950247454e-06
    };
    std::vector<double> A0_95 = {
        2.81765219760708e-14, 0.949999999984115, 1.00000000000224, 7.52388680069159e-13,
        -1.53428775400487e-14, -5.11748962582201e-15, -4.91220038116722e-12, 1.00000000000126,
        3.03494143714445e-05, -3.04088719562547e-07, -0.0168253328887104, -0.529950219190337,
        1.57853700431562e-07, -0.0149473720044402, -0.00131493201047327, -0.0447685337011653
    };
    std::vector<double> B0_95 = {
        2.56325892883034e-16, 3.83418999024824e-20, 1.04354459101730e-17,
        -1.65215169248193e-15, -2.68898299999025e-19, -1.71444762480191e-17,
        -9.42616979225867e-05, -2.21991902296806e-06, 5.92585596139766e-05,
        1.19102840386191e-05, 1.06742001052426e-17, 2.62436413593423e-06
    };

    std::vector<double> A0_90 = {
        1.15450414164574e-14, 0.899999999984989, 1.00000000000204, -7.79682694039008e-13,
        -1.24446690287617e-14, 2.22241445327327e-14, -4.40598675192542e-12, 1.00000000000094,
        2.72404661492861e-05, -3.04088708432438e-07, -0.0159398019313084, -0.500177519719206,
        1.41648450752593e-07, -0.0149473720044938, -0.00124957056733269, -0.0423925418025225
    };
    std::vector<double> B0_90 = {
        2.28746677229889e-16, 4.00327164639288e-20, 9.54817889792410e-18,
        -1.37205074606063e-15, -2.41147025112716e-19, -1.73234119581265e-17,
        -8.46005266228099e-05, -2.21991902295363e-06, 5.31849634319056e-05,
        1.06887655552615e-05, 9.60428195582752e-18, 2.35535007021971e-06
    };

    std::vector<double> A0_85 = {
        5.18865570204514e-15, 0.849999999985863, 1.00000000000177, 9.23643912344203e-14,
        -9.67133850999817e-15, 9.43217722010991e-15, -3.95198905551515e-12, 1.00000000000062,
        2.42995691906679e-05, -3.04088714912692e-07, -0.0150542705022785, -0.470404853993152,
        1.26319176748058e-07, -0.0149473720045853, -0.00118378185390699, -0.0400187458664940
    };
    std::vector<double> B0_85 = {
        2.80553677930112e-16, 3.80576165732420e-20, 1.14900019295046e-17,
        -1.97494550194419e-15, -2.96659755119512e-19, -1.67966633440528e-17,
        -0.000104445094693318, -2.21991902293735e-06, 6.56604582911099e-05,
        1.31978305491429e-05, 1.18166620493253e-17, 2.90791950247454e-06
    };

    std::vector<double> A0_80 = {
        2.30841159199512e-14, 0.799999999986717, 1.00000000000145, 1.06516117351754e-12,
        -8.17181032729845e-15, -5.66773662447807e-14, -3.49312205876576e-12, 1.00000000000059,
        2.15267234920454e-05, -3.04088686906005e-07, -0.0141686977241166, -0.440633076541411,
        1.11865875731301e-07, -0.0149473718293867, -0.00111756922441416, -0.0376470188627491
    };
    std::vector<double> B0_80 = {
        1.70615939518478e-16, 3.37897340071290e-20, 8.20669046771010e-18,
        -9.23208303556196e-16, -1.93522482563141e-19, -1.70931248755918e-17,
        -6.68448604456031e-05, -2.21991902295831e-06, 4.20226785575650e-05,
        8.44381266988425e-06, 7.67664317896266e-18, 1.86094584136560e-06
    };

    std::vector<double> A0_75 = {
        -8.58763405628733e-15, 0.749999999987507, 1.00000000000126, 5.24637570439804e-13,
        -7.16407172715310e-15, 3.03371185631579e-15, -3.13088333744857e-12, 1.00000000000028,
        1.89219290532415e-05, -3.04088724622658e-07, -0.0132831682878979, -0.410860335184786,
        9.82885454201937e-08, -0.0149473718294534, -0.00105092540048681, -0.0352776308623462
    };
    std::vector<double> B0_75 = {
        1.42660695974616e-16, 3.05287869837089e-20, 7.43694927689276e-18,
        -7.47056833700591e-16, -1.71976130921692e-19, -1.67966633440528e-17,
        -5.87503655686894e-05, -2.21991902297209e-06, 3.69339898650402e-05,
        7.42037826731075e-06, 6.81561599511456e-18, 1.63555567818956e-06
    };

    std::vector<double> A0_70 = {
        -1.28728753653248e-15, 0.699999999988350, 1.00000000000100, 6.70145320801319e-13,
        -5.56918872550381e-15, 4.60400246809677e-15, -2.77707283109399e-12, 1.00000000000020,
        1.64851858730402e-05, -3.04088725010739e-07, -0.0123976130369639, -0.381088218709117,
        8.55871840117989e-08, -0.0149473720047365, -0.000983856633089700, -0.0329103548960700
    };
    std::vector<double> B0_70 = {
        1.16570136665316e-16, 3.50135742652209e-20, 6.64073830647372e-18,
        -5.98196675954962e-16, -1.48378527700891e-19, -1.64260864296290e-17,
        -5.11780961662260e-05, -2.21991902297510e-06, 3.21736036686134e-05,
        6.46297189120608e-06, 6.02109335742837e-18, 1.42470681584066e-06
    };

    std::vector<double> A0_65 = {
        -7.70979403983022e-15, 0.649999999989156, 1.00000000000087, -5.18349350355977e-13,
        -4.85164475577985e-15, -1.06871829222554e-14, -2.40873519107290e-12, 1.00000000000013,
        1.42164939487926e-05, -3.04088716216424e-07, -0.0115120730729249, -0.351315700310901,
        7.37617891903504e-08, -0.0149473719170363, -0.000916358874814622, -0.0305453328145497
    };
    std::vector<double> B0_65 = {
        9.23097315574223e-17, 1.62932760822367e-20, 6.11134271443978e-18,
        -4.73678759193084e-16, -1.30378407234157e-19, -1.59411600673259e-17,
        -4.41280522385503e-05, -2.21991902296228e-06, 2.77415199682927e-05,
        5.57159354174907e-06, 5.29112874806112e-18, 1.22839925434790e-06
    };

    std::vector<double> A0_60 = {
        7.31765471935343e-15, 0.599999999990009, 1.00000000000057, 9.18559467541714e-13,
        -3.80135876809668e-15, -8.17170684579071e-15, -2.11833068278245e-12, 0.999999999999985,
        1.21158532793202e-05, -3.04088719890130e-07, -0.0106265346302429, -0.321543150667797,
        6.28123596628794e-08, -0.0149473719171854, -0.000848433589276270, -0.0281825148509628
    };
    std::vector<double> B0_60 = {
        7.18340932596561e-17, 2.87131212532013e-20, 5.23678119639971e-18,
        -3.70365826257214e-16, -1.14175765869462e-19, -1.54287051342371e-17,
        -3.76002337867158e-05, -2.21991902297133e-06, 2.36377387645675e-05,
        4.74624321867716e-06, 4.62296205428252e-18, 1.04663299371330e-06
    };

    std::vector<double> A0_55 = {
        -1.15577146119396e-14, 0.549999999990845, 1.00000000000054, 5.43141152142618e-13,
        -2.41136179425222e-15, 3.40942261895449e-15, -1.82998100789672e-12, 0.999999999999801,
        1.01832638645940e-05, -3.04088715936021e-07, -0.00974098606057166, -0.291770872577873,
        5.27388927958163e-08, -0.0149473719609998, -0.000780081971177142, -0.0258218598091893
    };
    std::vector<double> B0_55 = {
        5.37642489370515e-17, 1.76458620797045e-20, 4.67641596223179e-18,
        -2.85544003875344e-16, -9.77647433254584e-20, -1.48781337185218e-17,
        -3.15946408101810e-05, -2.21991902296127e-06, 1.98622600569533e-05,
        3.98692092178220e-06, 4.01411214354175e-18, 8.79408033912782e-07
    };

    std::vector<double> A0_50 = {
        1.66021359753550e-14, 0.499999999991675, 1.00000000000043, 5.98325410415346e-15,
        -1.82482858067303e-15, -3.05656571592921e-15, -1.57611568909862e-12, 0.999999999999921,
        8.41872570293982e-06, -3.04088711561443e-07, -0.00885544552149369, -0.261998360994712,
        4.35413861528086e-08, -0.0149473719609644, -0.000711302025464217, -0.0234634346893551
    };
    std::vector<double> B0_50 = {
        3.87319603271525e-17, 1.53034272997656e-20, 4.00117188459250e-18,
        -2.16626733818965e-16, -8.67603596665949e-20, -1.43148568085977e-17,
        -2.61112733097069e-05, -2.21991902296800e-06, 1.64150838458078e-05,
        3.29362665106758e-06, 3.46211912818249e-18, 7.26724374963059e-07
    };



    if (speed >= 5 && speed < 6) {
        A = A5;
        B = B5;
    } else if (speed >= 4.95 && speed < 5) {
        double ratio = (speed - 4.95) / (5.0 - 4.95);
        interpolateMatrices(A4_95, A5, B4_95, B5, ratio, A, B);
    } else if (speed >= 4.90 && speed < 4.95) {
        double ratio = (speed - 4.90) / (4.95 - 4.90);
        interpolateMatrices(A4_90, A4_95, B4_90, B4_95, ratio, A, B);
    } else if (speed >= 4.85 && speed < 4.90) {
        double ratio = (speed - 4.85) / (4.90 - 4.85);
        interpolateMatrices(A4_85, A4_90, B4_85, B4_90, ratio, A, B);
    } else if (speed >= 4.80 && speed < 4.85) {
        double ratio = (speed - 4.80) / (4.85 - 4.80);
        interpolateMatrices(A4_80, A4_85, B4_80, B4_85, ratio, A, B);
    } else if (speed >= 4.75 && speed < 4.80) {
        double ratio = (speed - 4.75) / (4.80 - 4.75);
        interpolateMatrices(A4_75, A4_80, B4_75, B4_80, ratio, A, B);
    } else if (speed >= 4.70 && speed < 4.75) {
        double ratio = (speed - 4.70) / (4.75 - 4.70);
        interpolateMatrices(A4_70, A4_75, B4_70, B4_75, ratio, A, B);
    } else if (speed >= 4.65 && speed < 4.70) {
        double ratio = (speed - 4.65) / (4.70 - 4.65);
        interpolateMatrices(A4_65, A4_70, B4_65, B4_70, ratio, A, B);
    } else if (speed >= 4.60 && speed < 4.65) {
        double ratio = (speed - 4.60) / (4.65 - 4.60);
        interpolateMatrices(A4_60, A4_65, B4_60, B4_65, ratio, A, B);
    } else if (speed >= 4.55 && speed < 4.60) {
        double ratio = (speed - 4.55) / (4.60 - 4.55);
        interpolateMatrices(A4_55, A4_60, B4_55, B4_60, ratio, A, B);
    } else if (speed >= 4.50 && speed < 4.55) {
        double ratio = (speed - 4.50) / (4.55 - 4.50);
        interpolateMatrices(A4_50, A4_55, B4_50, B4_55, ratio, A, B);
    } else if (speed >= 4.45 && speed < 4.50) {
        double ratio = (speed - 4.45) / (4.50 - 4.45);
        interpolateMatrices(A4_45, A4_50, B4_45, B4_50, ratio, A, B);
    } else if (speed >= 4.40 && speed < 4.45) {
        double ratio = (speed - 4.40) / (4.45 - 4.40);
        interpolateMatrices(A4_40, A4_45, B4_40, B4_45, ratio, A, B);
    } else if (speed >= 4.35 && speed < 4.40) {
        double ratio = (speed - 4.35) / (4.40 - 4.35);
        interpolateMatrices(A4_35, A4_40, B4_35, B4_40, ratio, A, B);
    } else if (speed >= 4.30 && speed < 4.35) {
        double ratio = (speed - 4.30) / (4.35 - 4.30);
        interpolateMatrices(A4_30, A4_35, B4_30, B4_35, ratio, A, B);
    } else if (speed >= 4.25 && speed < 4.30) {
        double ratio = (speed - 4.25) / (4.30 - 4.25);
        interpolateMatrices(A4_25, A4_30, B4_25, B4_30, ratio, A, B);
    } else if (speed >= 4.20 && speed < 4.25) {
        double ratio = (speed - 4.20) / (4.25 - 4.20);
        interpolateMatrices(A4_20, A4_25, B4_20, B4_25, ratio, A, B);
    } else if (speed >= 4.15 && speed < 4.20) {
        double ratio = (speed - 4.15) / (4.20 - 4.15);
        interpolateMatrices(A4_15, A4_20, B4_15, B4_20, ratio, A, B);
    } else if (speed >= 4.10 && speed < 4.15) {
        double ratio = (speed - 4.10) / (4.15 - 4.10);
        interpolateMatrices(A4_10, A4_15, B4_10, B4_15, ratio, A, B);
    } else if (speed >= 4.05 && speed < 4.10) {
        double ratio = (speed - 4.05) / (4.10 - 4.05);
        interpolateMatrices(A4_05, A4_10, B4_05, B4_10, ratio, A, B);
    } else if (speed >= 4.00 && speed < 4.05) {
        double ratio = (speed - 4.00) / (4.05 - 4.00);
        interpolateMatrices(A4_00, A4_05, B4_00, B4_05, ratio, A, B);
    } else if (speed >= 3.95 && speed < 4.00) {
        double ratio = (speed - 3.95) / (4.00 - 3.95);
        interpolateMatrices(A3_95, A4_00, B3_95, B4_00, ratio, A, B);
    } else if (speed >= 3.90 && speed < 3.95) {
        double ratio = (speed - 3.90) / (3.95 - 3.90);
        interpolateMatrices(A3_90, A3_95, B3_90, B3_95, ratio, A, B);
    } else if (speed >= 3.85 && speed < 3.90) {
        double ratio = (speed - 3.85) / (3.90 - 3.85);
        interpolateMatrices(A3_85, A3_90, B3_85, B3_90, ratio, A, B);
    } else if (speed >= 3.80 && speed < 3.85) {
        double ratio = (speed - 3.80) / (3.85 - 3.80);
        interpolateMatrices(A3_80, A3_85, B3_80, B3_85, ratio, A, B);
    } else if (speed >= 3.75 && speed < 3.80) {
        double ratio = (speed - 3.75) / (3.80 - 3.75);
        interpolateMatrices(A3_75, A3_80, B3_75, B3_80, ratio, A, B);
    } else if (speed >= 3.70 && speed < 3.75) {
        double ratio = (speed - 3.70) / (3.75 - 3.70);
        interpolateMatrices(A3_70, A3_75, B3_70, B3_75, ratio, A, B);
    } else if (speed >= 3.65 && speed < 3.70) {
        double ratio = (speed - 3.65) / (3.70 - 3.65);
        interpolateMatrices(A3_65, A3_70, B3_65, B3_70, ratio, A, B);
    } else if (speed >= 3.60 && speed < 3.65) {
        double ratio = (speed - 3.60) / (3.65 - 3.60);
        interpolateMatrices(A3_60, A3_65, B3_60, B3_65, ratio, A, B);
    } else if (speed >= 3.55 && speed < 3.60) {
        double ratio = (speed - 3.55) / (3.60 - 3.55);
        interpolateMatrices(A3_55, A3_60, B3_55, B3_60, ratio, A, B);
    } else if (speed >= 3.50 && speed < 3.55) {
        double ratio = (speed - 3.50) / (3.55 - 3.50);
        interpolateMatrices(A3_50, A3_55, B3_50, B3_55, ratio, A, B);
    } else if (speed >= 3.45 && speed < 3.50) {
        double ratio = (speed - 3.45) / (3.50 - 3.45);
        interpolateMatrices(A3_45, A3_50, B3_45, B3_50, ratio, A, B);
    } else if (speed >= 3.40 && speed < 3.45) {
        double ratio = (speed - 3.40) / (3.45 - 3.40);
        interpolateMatrices(A3_40, A3_45, B3_40, B3_45, ratio, A, B);
    } else if (speed >= 3.35 && speed < 3.40) {
        double ratio = (speed - 3.35) / (3.40 - 3.35);
        interpolateMatrices(A3_35, A3_40, B3_35, B3_40, ratio, A, B);
    } else if (speed >= 3.30 && speed < 3.35) {
        double ratio = (speed - 3.30) / (3.35 - 3.30);
        interpolateMatrices(A3_30, A3_35, B3_30, B3_35, ratio, A, B);
    } else if (speed >= 3.25 && speed < 3.30) {
        double ratio = (speed - 3.25) / (3.30 - 3.25);
        interpolateMatrices(A3_25, A3_30, B3_25, B3_30, ratio, A, B);
    } else if (speed >= 3.20 && speed < 3.25) {
        double ratio = (speed - 3.20) / (3.25 - 3.20);
        interpolateMatrices(A3_20, A3_25, B3_20, B3_25, ratio, A, B);
    } else if (speed >= 3.15 && speed < 3.20) {
        double ratio = (speed - 3.15) / (3.20 - 3.15);
        interpolateMatrices(A3_15, A3_20, B3_15, B3_20, ratio, A, B);
    } else if (speed >= 3.10 && speed < 3.15) {
        double ratio = (speed - 3.10) / (3.15 - 3.10);
        interpolateMatrices(A3_10, A3_15, B3_10, B3_15, ratio, A, B);
    } else if (speed >= 3.05 && speed < 3.10) {
        double ratio = (speed - 3.05) / (3.10 - 3.05);
        interpolateMatrices(A3_05, A3_10, B3_05, B3_10, ratio, A, B);
    } else if (speed >= 3.00 && speed < 3.05) {
        double ratio = (speed - 3.00) / (3.05 - 3.00);
        interpolateMatrices(A3_00, A3_05, B3_00, B3_05, ratio, A, B);
    } else if (speed >= 2.95 && speed < 3.00) {
        double ratio = (speed - 2.95) / (3.00 - 2.95);
        interpolateMatrices(A2_95, A3_00, B2_95, B3_00, ratio, A, B);
    } else if (speed >= 2.90 && speed < 2.95) {
        double ratio = (speed - 2.90) / (2.95 - 2.90);
        interpolateMatrices(A2_90, A2_95, B2_90, B2_95, ratio, A, B);
    } else if (speed >= 2.85 && speed < 2.90) {
        double ratio = (speed - 2.85) / (2.90 - 2.85);
        interpolateMatrices(A2_85, A2_90, B2_85, B2_90, ratio, A, B);
    } else if (speed >= 2.80 && speed < 2.85) {
        double ratio = (speed - 2.80) / (2.85 - 2.80);
        interpolateMatrices(A2_80, A2_85, B2_80, B2_85, ratio, A, B);
    } else if (speed >= 2.75 && speed < 2.80) {
        double ratio = (speed - 2.75) / (2.80 - 2.75);
        interpolateMatrices(A2_75, A2_80, B2_75, B2_80, ratio, A, B);
    } else if (speed >= 2.70 && speed < 2.75) {
        double ratio = (speed - 2.70) / (2.75 - 2.70);
        interpolateMatrices(A2_70, A2_75, B2_70, B2_75, ratio, A, B);
    } else if (speed >= 2.65 && speed < 2.70) {
        double ratio = (speed - 2.65) / (2.70 - 2.65);
        interpolateMatrices(A2_65, A2_70, B2_65, B2_70, ratio, A, B);
    } else if (speed >= 2.60 && speed < 2.65) {
        double ratio = (speed - 2.60) / (2.65 - 2.60);
        interpolateMatrices(A2_60, A2_65, B2_60, B2_65, ratio, A, B);
    } else if (speed >= 2.55 && speed < 2.60) {
        double ratio = (speed - 2.55) / (2.60 - 2.55);
        interpolateMatrices(A2_55, A2_60, B2_55, B2_60, ratio, A, B);
    } else if (speed >= 2.50 && speed < 2.55) {
        double ratio = (speed - 2.50) / (2.55 - 2.50);
        interpolateMatrices(A2_50, A2_55, B2_50, B2_55, ratio, A, B);
    } else if (speed >= 2.45 && speed < 2.50) {
        double ratio = (speed - 2.45) / (2.50 - 2.45);
        interpolateMatrices(A2_45, A2_50, B2_45, B2_50, ratio, A, B);
    } else if (speed >= 2.40 && speed < 2.45) {
        double ratio = (speed - 2.40) / (2.45 - 2.40);
        interpolateMatrices(A2_40, A2_45, B2_40, B2_45, ratio, A, B);
    } else if (speed >= 2.35 && speed < 2.40) {
        double ratio = (speed - 2.35) / (2.40 - 2.35);
        interpolateMatrices(A2_35, A2_40, B2_35, B2_40, ratio, A, B);
    } else if (speed >= 2.30 && speed < 2.35) {
        double ratio = (speed - 2.30) / (2.35 - 2.30);
        interpolateMatrices(A2_30, A2_35, B2_30, B2_35, ratio, A, B);
    } else if (speed >= 2.25 && speed < 2.30) {
        double ratio = (speed - 2.25) / (2.30 - 2.25);
        interpolateMatrices(A2_25, A2_30, B2_25, B2_30, ratio, A, B);
    } else if (speed >= 2.20 && speed < 2.25) {
        double ratio = (speed - 2.20) / (2.25 - 2.20);
        interpolateMatrices(A2_20, A2_25, B2_20, B2_25, ratio, A, B);
    } else if (speed >= 2.15 && speed < 2.20) {
        double ratio = (speed - 2.15) / (2.20 - 2.15);
        interpolateMatrices(A2_15, A2_20, B2_15, B2_20, ratio, A, B);
    } else if (speed >= 2.10 && speed < 2.15) {
        double ratio = (speed - 2.10) / (2.15 - 2.10);
        interpolateMatrices(A2_10, A2_15, B2_10, B2_15, ratio, A, B);
    } else if (speed >= 2.05 && speed < 2.10) {
        double ratio = (speed - 2.05) / (2.10 - 2.05);
        interpolateMatrices(A2_05, A2_10, B2_05, B2_10, ratio, A, B);
    } else if (speed >= 2.00 && speed < 2.05) {
        double ratio = (speed - 2.00) / (2.05 - 2.00);
        interpolateMatrices(A2_00, A2_05, B2_00, B2_05, ratio, A, B);
    } else if (speed >= 1.95 && speed < 2.00) {
        double ratio = (speed - 1.95) / (2.00 - 1.95);
        interpolateMatrices(A1_95, A2_00, B1_95, B2_00, ratio, A, B);
    } else if (speed >= 1.90 && speed < 1.95) {
        double ratio = (speed - 1.90) / (1.95 - 1.90);
        interpolateMatrices(A1_90, A1_95, B1_90, B1_95, ratio, A, B);
    } else if (speed >= 1.85 && speed < 1.90) {
        double ratio = (speed - 1.85) / (1.90 - 1.85);
        interpolateMatrices(A1_85, A1_90, B1_85, B1_90, ratio, A, B);
    } else if (speed >= 1.80 && speed < 1.85) {
        double ratio = (speed - 1.80) / (1.85 - 1.80);
        interpolateMatrices(A1_80, A1_85, B1_80, B1_85, ratio, A, B);
    } else if (speed >= 1.75 && speed < 1.80) {
        double ratio = (speed - 1.75) / (1.80 - 1.75);
        interpolateMatrices(A1_75, A1_80, B1_75, B1_80, ratio, A, B);
    } else if (speed >= 1.70 && speed < 1.75) {
        double ratio = (speed - 1.70) / (1.75 - 1.70);
        interpolateMatrices(A1_70, A1_75, B1_70, B1_75, ratio, A, B);
    } else if (speed >= 1.65 && speed < 1.70) {
        double ratio = (speed - 1.65) / (1.70 - 1.65);
        interpolateMatrices(A1_65, A1_70, B1_65, B1_70, ratio, A, B);
    } else if (speed >= 1.60 && speed < 1.65) {
        double ratio = (speed - 1.60) / (1.65 - 1.60);
        interpolateMatrices(A1_60, A1_65, B1_60, B1_65, ratio, A, B);
    } else if (speed >= 1.55 && speed < 1.60) {
        double ratio = (speed - 1.55) / (1.60 - 1.55);
        interpolateMatrices(A1_55, A1_60, B1_55, B1_60, ratio, A, B);
    } else if (speed >= 1.50 && speed < 1.55) {
        double ratio = (speed - 1.50) / (1.55 - 1.50);
        interpolateMatrices(A1_50, A1_55, B1_50, B1_55, ratio, A, B);
    } else if (speed >= 1.45 && speed < 1.50) {
        double ratio = (speed - 1.45) / (1.50 - 1.45);
        interpolateMatrices(A1_45, A1_50, B1_45, B1_50, ratio, A, B);
    } else if (speed >= 1.40 && speed < 1.45) {
        double ratio = (speed - 1.40) / (1.45 - 1.40);
        interpolateMatrices(A1_40, A1_45, B1_40, B1_45, ratio, A, B);
    } else if (speed >= 1.35 && speed < 1.40) {
        double ratio = (speed - 1.35) / (1.40 - 1.35);
        interpolateMatrices(A1_35, A1_40, B1_35, B1_40, ratio, A, B);
    } else if (speed >= 1.30 && speed < 1.35) {
        double ratio = (speed - 1.30) / (1.35 - 1.30);
        interpolateMatrices(A1_30, A1_35, B1_30, B1_35, ratio, A, B);
    } else if (speed >= 1.25 && speed < 1.30) {
        double ratio = (speed - 1.25) / (1.30 - 1.25);
        interpolateMatrices(A1_25, A1_30, B1_25, B1_30, ratio, A, B);
    } else if (speed >= 1.20 && speed < 1.25) {
        double ratio = (speed - 1.20) / (1.25 - 1.20);
        interpolateMatrices(A1_20, A1_25, B1_20, B1_25, ratio, A, B);
    } else if (speed >= 1.15 && speed < 1.20) {
        double ratio = (speed - 1.15) / (1.20 - 1.15);
        interpolateMatrices(A1_15, A1_20, B1_15, B1_20, ratio, A, B);
    } else if (speed >= 1.10 && speed < 1.15) {
        double ratio = (speed - 1.10) / (1.15 - 1.10);
        interpolateMatrices(A1_10, A1_15, B1_10, B1_15, ratio, A, B);
    } else if (speed >= 1.05 && speed < 1.10) {
        double ratio = (speed - 1.05) / (1.10 - 1.05);
        interpolateMatrices(A1_05, A1_10, B1_05, B1_10, ratio, A, B);
    } else if (speed >= 1.00 && speed < 1.05) {
        double ratio = (speed - 1.00) / (1.05 - 1.00);
        interpolateMatrices(A1_00, A1_05, B1_00, B1_05, ratio, A, B);
    } else if (speed >= 0.95 && speed < 1.00) {
        double ratio = (speed - 0.95) / (1.00 - 0.95);
        interpolateMatrices(A0_95, A1_00, B0_95, B1_00, ratio, A, B);
    } else if (speed >= 0.90 && speed < 0.95) {
        double ratio = (speed - 0.90) / (0.95 - 0.90);
        interpolateMatrices(A0_90, A0_95, B0_90, B0_95, ratio, A, B);
    } else if (speed >= 0.85 && speed < 0.90) {
        double ratio = (speed - 0.85) / (0.90 - 0.85);
        interpolateMatrices(A0_85, A0_90, B0_85, B0_90, ratio, A, B);
    } else if (speed >= 0.80 && speed < 0.85) {
        double ratio = (speed - 0.80) / (0.85 - 0.80);
        interpolateMatrices(A0_80, A0_85, B0_80, B0_85, ratio, A, B);
    } else if (speed >= 0.75 && speed < 0.80) {
        double ratio = (speed - 0.75) / (0.80 - 0.75);
        interpolateMatrices(A0_75, A0_80, B0_75, B0_80, ratio, A, B);
    } else if (speed >= 0.70 && speed < 0.75) {
        double ratio = (speed - 0.70) / (0.75 - 0.70);
        interpolateMatrices(A0_70, A0_75, B0_70, B0_75, ratio, A, B);
    } else if (speed >= 0.65 && speed < 0.70) {
        double ratio = (speed - 0.65) / (0.70 - 0.65);
        interpolateMatrices(A0_65, A0_70, B0_65, B0_70, ratio, A, B);
    } else if (speed >= 0.60 && speed < 0.65) {
        double ratio = (speed - 0.60) / (0.65 - 0.60);
        interpolateMatrices(A0_60, A0_65, B0_60, B0_65, ratio, A, B);
    } else if (speed >= 0.55 && speed < 0.60) {
        double ratio = (speed - 0.55) / (0.60 - 0.55);
        interpolateMatrices(A0_55, A0_60, B0_55, B0_60, ratio, A, B);
    } else if (speed >= 0.50 && speed < 0.55) {
        double ratio = (speed - 0.50) / (0.55 - 0.50);
        interpolateMatrices(A0_50, A0_55, B0_50, B0_55, ratio, A, B);
    } else {
        A = A0_50;
        B = B0_50;
    }

}