%typeset_mode True
var('xi, eta') x_nodes = [-1, 1] xi_LGL = [-1.0, \ -0.87174014851, \ -0.591700181433,\ -0.209299217902,\ 0.209299217902, \ 0.591700181433, \ 0.87174014851, \ 1.0] eta_LGL = [-1.0, \ -0.87174014851, \ -0.591700181433,\ -0.209299217902,\ 0.209299217902, \ 0.591700181433, \ 0.87174014851, \ 1.0]
(ξ, η)
#L_i(\xi) L_0_xi = ((xi - (xi_LGL[1])) / (xi_LGL[0] - xi_LGL[1]))\ * ((xi - (xi_LGL[2])) / (xi_LGL[0] - xi_LGL[2]))\ * ((xi - (xi_LGL[3])) / (xi_LGL[0] - xi_LGL[3]))\ * ((xi - (xi_LGL[4])) / (xi_LGL[0] - xi_LGL[4]))\ * ((xi - (xi_LGL[5])) / (xi_LGL[0] - xi_LGL[5]))\ * ((xi - (xi_LGL[6])) / (xi_LGL[0] - xi_LGL[6]))\ * ((xi - (xi_LGL[7])) / (xi_LGL[0] - xi_LGL[7])) L_1_xi = ((xi - (xi_LGL[0])) / (xi_LGL[1] - xi_LGL[0]))\ * ((xi - (xi_LGL[2])) / (xi_LGL[1] - xi_LGL[2]))\ * ((xi - (xi_LGL[3])) / (xi_LGL[1] - xi_LGL[3]))\ * ((xi - (xi_LGL[4])) / (xi_LGL[1] - xi_LGL[4]))\ * ((xi - (xi_LGL[5])) / (xi_LGL[1] - xi_LGL[5]))\ * ((xi - (xi_LGL[6])) / (xi_LGL[1] - xi_LGL[6]))\ * ((xi - (xi_LGL[7])) / (xi_LGL[1] - xi_LGL[7])) L_2_xi = ((xi - (xi_LGL[0])) / (xi_LGL[2] - xi_LGL[0])) \ * ((xi - (xi_LGL[1])) / (xi_LGL[2] - xi_LGL[1])) \ * ((xi - (xi_LGL[3])) / (xi_LGL[2] - xi_LGL[3])) \ * ((xi - (xi_LGL[4])) / (xi_LGL[2] - xi_LGL[4])) \ * ((xi - (xi_LGL[5])) / (xi_LGL[2] - xi_LGL[5])) \ * ((xi - (xi_LGL[6])) / (xi_LGL[2] - xi_LGL[6])) \ * ((xi - (xi_LGL[7])) / (xi_LGL[2] - xi_LGL[7])) L_3_xi = ((xi - (xi_LGL[0])) / (xi_LGL[3] - xi_LGL[0])) \ * ((xi - (xi_LGL[1])) / (xi_LGL[3] - xi_LGL[1])) \ * ((xi - (xi_LGL[2])) / (xi_LGL[3] - xi_LGL[2])) \ * ((xi - (xi_LGL[4])) / (xi_LGL[3] - xi_LGL[4])) \ * ((xi - (xi_LGL[5])) / (xi_LGL[3] - xi_LGL[5])) \ * ((xi - (xi_LGL[6])) / (xi_LGL[3] - xi_LGL[6])) \ * ((xi - (xi_LGL[7])) / (xi_LGL[3] - xi_LGL[7])) L_4_xi = ((xi - (xi_LGL[0])) / (xi_LGL[4] - xi_LGL[0])) \ * ((xi - (xi_LGL[1])) / (xi_LGL[4] - xi_LGL[1])) \ * ((xi - (xi_LGL[2])) / (xi_LGL[4] - xi_LGL[2])) \ * ((xi - (xi_LGL[3])) / (xi_LGL[4] - xi_LGL[3])) \ * ((xi - (xi_LGL[5])) / (xi_LGL[4] - xi_LGL[5])) \ * ((xi - (xi_LGL[6])) / (xi_LGL[4] - xi_LGL[6])) \ * ((xi - (xi_LGL[7])) / (xi_LGL[4] - xi_LGL[7])) L_5_xi = ((xi - (xi_LGL[0])) / (xi_LGL[5] - xi_LGL[0])) \ * ((xi - (xi_LGL[1])) / (xi_LGL[5] - xi_LGL[1])) \ * ((xi - (xi_LGL[2])) / (xi_LGL[5] - xi_LGL[2])) \ * ((xi - (xi_LGL[3])) / (xi_LGL[5] - xi_LGL[3])) \ * ((xi - (xi_LGL[4])) / (xi_LGL[5] - xi_LGL[4])) \ * ((xi - (xi_LGL[6])) / (xi_LGL[5] - xi_LGL[6])) \ * ((xi - (xi_LGL[7])) / (xi_LGL[5] - xi_LGL[7])) L_6_xi = ((xi - (xi_LGL[0])) / (xi_LGL[6] - xi_LGL[0])) \ * ((xi - (xi_LGL[1])) / (xi_LGL[6] - xi_LGL[1])) \ * ((xi - (xi_LGL[2])) / (xi_LGL[6] - xi_LGL[2])) \ * ((xi - (xi_LGL[3])) / (xi_LGL[6] - xi_LGL[3])) \ * ((xi - (xi_LGL[4])) / (xi_LGL[6] - xi_LGL[4])) \ * ((xi - (xi_LGL[5])) / (xi_LGL[6] - xi_LGL[5])) \ * ((xi - (xi_LGL[7])) / (xi_LGL[6] - xi_LGL[7])) L_7_xi = ((xi - (xi_LGL[0])) / (xi_LGL[7] - xi_LGL[0])) \ * ((xi - (xi_LGL[1])) / (xi_LGL[7] - xi_LGL[1])) \ * ((xi - (xi_LGL[2])) / (xi_LGL[7] - xi_LGL[2])) \ * ((xi - (xi_LGL[3])) / (xi_LGL[7] - xi_LGL[3])) \ * ((xi - (xi_LGL[4])) / (xi_LGL[7] - xi_LGL[4])) \ * ((xi - (xi_LGL[5])) / (xi_LGL[7] - xi_LGL[5])) \ * ((xi - (xi_LGL[6])) / (xi_LGL[7] - xi_LGL[6])) \ # L_i(\eta) L_0_eta = ((eta - (xi_LGL[1])) / (xi_LGL[0] - eta_LGL[1]))\ * ((eta - (xi_LGL[2])) / (xi_LGL[0] - eta_LGL[2]))\ * ((eta - (xi_LGL[3])) / (xi_LGL[0] - eta_LGL[3]))\ * ((eta - (xi_LGL[4])) / (xi_LGL[0] - eta_LGL[4]))\ * ((eta - (xi_LGL[5])) / (xi_LGL[0] - eta_LGL[5]))\ * ((eta - (xi_LGL[6])) / (xi_LGL[0] - eta_LGL[6]))\ * ((eta - (xi_LGL[7])) / (xi_LGL[0] - eta_LGL[7])) L_1_eta = ((eta - (xi_LGL[0])) / (xi_LGL[1] - eta_LGL[0]))\ * ((eta - (xi_LGL[2])) / (xi_LGL[1] - eta_LGL[2]))\ * ((eta - (xi_LGL[3])) / (xi_LGL[1] - eta_LGL[3]))\ * ((eta - (xi_LGL[4])) / (xi_LGL[1] - eta_LGL[4]))\ * ((eta - (xi_LGL[5])) / (xi_LGL[1] - eta_LGL[5]))\ * ((eta - (xi_LGL[6])) / (xi_LGL[1] - eta_LGL[6]))\ * ((eta - (xi_LGL[7])) / (xi_LGL[1] - eta_LGL[7])) L_2_eta = ((eta - (xi_LGL[0])) / (xi_LGL[2] - eta_LGL[0])) \ * ((eta - (xi_LGL[1])) / (xi_LGL[2] - eta_LGL[1])) \ * ((eta - (xi_LGL[3])) / (xi_LGL[2] - eta_LGL[3])) \ * ((eta - (xi_LGL[4])) / (xi_LGL[2] - eta_LGL[4])) \ * ((eta - (xi_LGL[5])) / (xi_LGL[2] - eta_LGL[5])) \ * ((eta - (xi_LGL[6])) / (xi_LGL[2] - eta_LGL[6])) \ * ((eta - (xi_LGL[7])) / (xi_LGL[2] - eta_LGL[7])) L_3_eta = ((eta - (xi_LGL[0])) / (xi_LGL[3] - eta_LGL[0])) \ * ((eta - (xi_LGL[1])) / (xi_LGL[3] - eta_LGL[1])) \ * ((eta - (xi_LGL[2])) / (xi_LGL[3] - eta_LGL[2])) \ * ((eta - (xi_LGL[4])) / (xi_LGL[3] - eta_LGL[4])) \ * ((eta - (xi_LGL[5])) / (xi_LGL[3] - eta_LGL[5])) \ * ((eta - (xi_LGL[6])) / (xi_LGL[3] - eta_LGL[6])) \ * ((eta - (xi_LGL[7])) / (xi_LGL[3] - eta_LGL[7])) L_4_eta = ((eta - (xi_LGL[0])) / (xi_LGL[4] - eta_LGL[0])) \ * ((eta - (xi_LGL[1])) / (xi_LGL[4] - eta_LGL[1])) \ * ((eta - (xi_LGL[2])) / (xi_LGL[4] - eta_LGL[2])) \ * ((eta - (xi_LGL[3])) / (xi_LGL[4] - eta_LGL[3])) \ * ((eta - (xi_LGL[5])) / (xi_LGL[4] - eta_LGL[5])) \ * ((eta - (xi_LGL[6])) / (xi_LGL[4] - eta_LGL[6])) \ * ((eta - (xi_LGL[7])) / (xi_LGL[4] - eta_LGL[7])) L_5_eta = ((eta - (xi_LGL[0])) / (xi_LGL[5] - eta_LGL[0])) \ * ((eta - (xi_LGL[1])) / (xi_LGL[5] - eta_LGL[1])) \ * ((eta - (xi_LGL[2])) / (xi_LGL[5] - eta_LGL[2])) \ * ((eta - (xi_LGL[3])) / (xi_LGL[5] - eta_LGL[3])) \ * ((eta - (xi_LGL[4])) / (xi_LGL[5] - eta_LGL[4])) \ * ((eta - (xi_LGL[6])) / (xi_LGL[5] - eta_LGL[6])) \ * ((eta - (xi_LGL[7])) / (xi_LGL[5] - eta_LGL[7])) L_6_eta = ((eta - (xi_LGL[0])) / (xi_LGL[6] - eta_LGL[0])) \ * ((eta - (xi_LGL[1])) / (xi_LGL[6] - eta_LGL[1])) \ * ((eta - (xi_LGL[2])) / (xi_LGL[6] - eta_LGL[2])) \ * ((eta - (xi_LGL[3])) / (xi_LGL[6] - eta_LGL[3])) \ * ((eta - (xi_LGL[4])) / (xi_LGL[6] - eta_LGL[4])) \ * ((eta - (xi_LGL[5])) / (xi_LGL[6] - eta_LGL[5])) \ * ((eta - (xi_LGL[7])) / (xi_LGL[6] - eta_LGL[7])) L_7_eta = ((eta - (xi_LGL[0])) / (xi_LGL[7] - eta_LGL[0])) \ * ((eta - (xi_LGL[1])) / (xi_LGL[7] - eta_LGL[1])) \ * ((eta - (xi_LGL[2])) / (xi_LGL[7] - eta_LGL[2])) \ * ((eta - (xi_LGL[3])) / (xi_LGL[7] - eta_LGL[3])) \ * ((eta - (xi_LGL[4])) / (xi_LGL[7] - eta_LGL[4])) \ * ((eta - (xi_LGL[5])) / (xi_LGL[7] - eta_LGL[5])) \ * ((eta - (xi_LGL[6])) / (xi_LGL[7] - eta_LGL[6])) \
L_0_xi.simplify_full() L_1_xi.simplify_full() L_2_xi.simplify_full() L_3_xi.simplify_full() L_4_xi.simplify_full() L_5_xi.simplify_full() L_6_xi.simplify_full() L_7_xi.simplify_full() L_0_eta.simplify_full() L_1_eta.simplify_full() L_2_eta.simplify_full() L_3_eta.simplify_full() L_4_eta.simplify_full() L_5_eta.simplify_full() L_6_eta.simplify_full() L_7_eta.simplify_full()
−3.35156250001ξ7+3.35156250001ξ6+3.86718750001ξ5−3.86718750001ξ4−1.0546875ξ3+1.0546875ξ2+0.0390624999999ξ−0.0390624999999
8.14072271825ξ7−7.09659483138ξ6−11.347477684ξ5+9.89205188146ξ4+3.33160871212ξ3−2.90429707348ξ2−0.124853746365ξ+0.108840023398
−10.3581368289ξ7+6.12891144098ξ6+18.6833551584ξ5−11.054944637ξ4−8.6700371412ξ3+5.13006254948ξ2+0.34481881174ξ−0.204029353468
11.3898137485ξ7−2.38387910961ξ6−24.032962502ξ5+5.03008025554ξ4+15.6735080469ξ3−3.280452976ξ2−3.0303592934ξ+0.63425183007
−11.3898137485ξ7−2.38387910961ξ6+24.032962502ξ5+5.03008025554ξ4−15.6735080469ξ3−3.280452976ξ2+3.0303592934ξ+0.63425183007
10.3581368289ξ7+6.12891144098ξ6−18.6833551584ξ5−11.054944637ξ4+8.6700371412ξ3+5.13006254948ξ2−0.34481881174ξ−0.204029353468
−8.14072271825ξ7−7.09659483138ξ6+11.347477684ξ5+9.89205188146ξ4−3.33160871212ξ3−2.90429707348ξ2+0.124853746365ξ+0.108840023398
3.35156250001ξ7+3.35156250001ξ6−3.86718750001ξ5−3.86718750001ξ4+1.0546875ξ3+1.0546875ξ2−0.0390624999999ξ−0.0390624999999
−3.35156250001η7+3.35156250001η6+3.86718750001η5−3.86718750001η4−1.0546875η3+1.0546875η2+0.0390624999999η−0.0390624999999
8.14072271825η7−7.09659483138η6−11.347477684η5+9.89205188146η4+3.33160871212η3−2.90429707348η2−0.124853746365η+0.108840023398
−10.3581368289η7+6.12891144098η6+18.6833551584η5−11.054944637η4−8.6700371412η3+5.13006254948η2+0.34481881174η−0.204029353468
11.3898137485η7−2.38387910961η6−24.032962502η5+5.03008025554η4+15.6735080469η3−3.280452976η2−3.0303592934η+0.63425183007
−11.3898137485η7−2.38387910961η6+24.032962502η5+5.03008025554η4−15.6735080469η3−3.280452976η2+3.0303592934η+0.63425183007
10.3581368289η7+6.12891144098η6−18.6833551584η5−11.054944637η4+8.6700371412η3+5.13006254948η2−0.34481881174η−0.204029353468
−8.14072271825η7−7.09659483138η6+11.347477684η5+9.89205188146η4−3.33160871212η3−2.90429707348η2+0.124853746365η+0.108840023398
3.35156250001η7+3.35156250001η6−3.86718750001η5−3.86718750001η4+1.0546875η3+1.0546875η2−0.0390624999999η−0.0390624999999
L_0_xi_eta = L_0_xi * L_0_eta L_0_xi_eta.simplify_full()
(11.2329711915η7−11.2329711915η6−12.9611206055η5+12.9611206055η4+3.53485107423η3−3.53485107423η2−0.130920410156η+0.130920410156)ξ7+0.130920410156η7+(−11.2329711915η7+11.2329711915η6+12.9611206055η5−12.9611206055η4−3.53485107423η3+3.53485107423η2+0.130920410156η−0.130920410156)ξ6−0.130920410156η6+(−12.9611206055η7+12.9611206055η6+14.9551391602η5−14.9551391602η4−4.07867431643η3+4.07867431643η2+0.151062011719η−0.151062011719)ξ5−0.151062011719η5+(12.9611206055η7−12.9611206055η6−14.9551391602η5+14.9551391602η4+4.07867431643η3−4.07867431643η2−0.151062011719η+0.151062011719)ξ4+0.151062011719η4+(3.53485107423η7−3.53485107423η6−4.07867431643η5+4.07867431643η4+1.11236572266η3−1.11236572266η2−0.0411987304688η+0.0411987304688)ξ3+0.0411987304688η3+(−3.53485107423η7+3.53485107423η6+4.07867431643η5−4.07867431643η4−1.11236572266η3+1.11236572266η2+0.0411987304688η−0.0411987304688)ξ2−0.0411987304688η2+(−0.130920410156η7+0.130920410156η6+0.151062011719η5−0.151062011719η4−0.0411987304688η3+0.0411987304688η2+0.00152587890624η−0.00152587890624)ξ−0.00152587890624η+0.00152587890624