Здравствуйте, у меня есть следующие данные (зависимость молярной теплоёмкость воды от температуры воды):
Код:
list = {{0., 4.217}, {1., 4.213}, {2., 4.21}, {3., 4.207}, {4., 4.205}, {5., 4.202}, {6., 4.2}, {7., 4.198}, {8., 4.196}, {9., 4.194}, {10., 4.192}, {11., 4.191}, {12., 4.189}, {13., 4.188}, {14., 4.187}, {15., 4.186}, {16., 4.185}, {17., 4.184}, {18., 4.183}, {19., 4.182}, {20., 4.182}, {21., 4.181}, {22., 4.181}, {23., 4.18}, {24., 4.18}, {25., 4.18}, {26., 4.179}, {27., 4.179}, {28., 4.179}, {29., 4.179}, {30., 4.178}, {31., 4.178}, {32., 4.178}, {33., 4.178}, {34., 4.178}, {35., 4.178}, {36., 4.178}, {37., 4.178}, {38., 4.178}, {39., 4.179}, {40., 4.179}, {41., 4.179}, {42., 4.179}, {43., 4.179}, {44., 4.179}, {45., 4.18}, {46., 4.18}, {47., 4.18}, {48., 4.18}, {49., 4.181}, {50., 4.181}, {51., 4.181}, {52., 4.182}, {53., 4.182}, {54., 4.182}, {55., 4.183}, {56., 4.183}, {57., 4.183}, {58., 4.184}, {59., 4.184}, {60., 4.185}, {61., 4.185}, {62., 4.186}, {63., 4.186}, {64., 4.187}, {65., 4.187}, {66., 4.188}, {67., 4.188}, {68., 4.189}, {69., 4.189}, {70., 4.19}, {71., 4.19}, {72., 4.191}, {73., 4.192}, {74., 4.192}, {75., 4.193}, {76., 4.194}, {77., 4.194}, {78., 4.195}, {79., 4.196}, {80., 4.196}, {81., 4.197}, {82., 4.198}, {83., 4.199}, {84., 4.2}, {85., 4.2}, {86., 4.201}, {87., 4.202}, {88., 4.203}, {89., 4.204}, {90., 4.205}, {91., 4.206}, {92., 4.207}, {93., 4.208}, {94., 4.209}, {95., 4.21}, {96., 4.211}, {97., 4.212}, {98., 4.213}, {99., 4.214}, {100., 4.216}};
Я добавляю к этим данным единицы измерения и интерполирую:
Код:
f = Interpolation@MapAt[Quantity[#1, ("Kilojoules")/("Kelvins" "Moles")] &, MapAt[Quantity[#1, "DegreesCelsius"] &, list, {;; , 1}], {;; , 2}];
Теперь интегрирую двумя способами:
Надеюсь на Вашу помощь. Есть идея, что подобное поведение в первом случае связано с разной точностью у входных данных, т. к.: