Contact
CoCalc Logo Icon
StoreFeaturesDocsShareSupport News AboutSign UpSign In
| Download
Views: 1021
%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]
(ξ\displaystyle \xi, η\displaystyle \eta)
#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ξ53.86718750001ξ41.0546875ξ3+1.0546875ξ2+0.0390624999999ξ0.0390624999999\displaystyle -3.35156250001 \, \xi^{7} + 3.35156250001 \, \xi^{6} + 3.86718750001 \, \xi^{5} - 3.86718750001 \, \xi^{4} - 1.0546875 \, \xi^{3} + 1.0546875 \, \xi^{2} + 0.0390624999999 \, \xi - 0.0390624999999
8.14072271825ξ77.09659483138ξ611.347477684ξ5+9.89205188146ξ4+3.33160871212ξ32.90429707348ξ20.124853746365ξ+0.108840023398\displaystyle 8.14072271825 \, \xi^{7} - 7.09659483138 \, \xi^{6} - 11.347477684 \, \xi^{5} + 9.89205188146 \, \xi^{4} + 3.33160871212 \, \xi^{3} - 2.90429707348 \, \xi^{2} - 0.124853746365 \, \xi + 0.108840023398
10.3581368289ξ7+6.12891144098ξ6+18.6833551584ξ511.054944637ξ48.6700371412ξ3+5.13006254948ξ2+0.34481881174ξ0.204029353468\displaystyle -10.3581368289 \, \xi^{7} + 6.12891144098 \, \xi^{6} + 18.6833551584 \, \xi^{5} - 11.054944637 \, \xi^{4} - 8.6700371412 \, \xi^{3} + 5.13006254948 \, \xi^{2} + 0.34481881174 \, \xi - 0.204029353468
11.3898137485ξ72.38387910961ξ624.032962502ξ5+5.03008025554ξ4+15.6735080469ξ33.280452976ξ23.0303592934ξ+0.63425183007\displaystyle 11.3898137485 \, \xi^{7} - 2.38387910961 \, \xi^{6} - 24.032962502 \, \xi^{5} + 5.03008025554 \, \xi^{4} + 15.6735080469 \, \xi^{3} - 3.280452976 \, \xi^{2} - 3.0303592934 \, \xi + 0.63425183007
11.3898137485ξ72.38387910961ξ6+24.032962502ξ5+5.03008025554ξ415.6735080469ξ33.280452976ξ2+3.0303592934ξ+0.63425183007\displaystyle -11.3898137485 \, \xi^{7} - 2.38387910961 \, \xi^{6} + 24.032962502 \, \xi^{5} + 5.03008025554 \, \xi^{4} - 15.6735080469 \, \xi^{3} - 3.280452976 \, \xi^{2} + 3.0303592934 \, \xi + 0.63425183007
10.3581368289ξ7+6.12891144098ξ618.6833551584ξ511.054944637ξ4+8.6700371412ξ3+5.13006254948ξ20.34481881174ξ0.204029353468\displaystyle 10.3581368289 \, \xi^{7} + 6.12891144098 \, \xi^{6} - 18.6833551584 \, \xi^{5} - 11.054944637 \, \xi^{4} + 8.6700371412 \, \xi^{3} + 5.13006254948 \, \xi^{2} - 0.34481881174 \, \xi - 0.204029353468
8.14072271825ξ77.09659483138ξ6+11.347477684ξ5+9.89205188146ξ43.33160871212ξ32.90429707348ξ2+0.124853746365ξ+0.108840023398\displaystyle -8.14072271825 \, \xi^{7} - 7.09659483138 \, \xi^{6} + 11.347477684 \, \xi^{5} + 9.89205188146 \, \xi^{4} - 3.33160871212 \, \xi^{3} - 2.90429707348 \, \xi^{2} + 0.124853746365 \, \xi + 0.108840023398
3.35156250001ξ7+3.35156250001ξ63.86718750001ξ53.86718750001ξ4+1.0546875ξ3+1.0546875ξ20.0390624999999ξ0.0390624999999\displaystyle 3.35156250001 \, \xi^{7} + 3.35156250001 \, \xi^{6} - 3.86718750001 \, \xi^{5} - 3.86718750001 \, \xi^{4} + 1.0546875 \, \xi^{3} + 1.0546875 \, \xi^{2} - 0.0390624999999 \, \xi - 0.0390624999999
3.35156250001η7+3.35156250001η6+3.86718750001η53.86718750001η41.0546875η3+1.0546875η2+0.0390624999999η0.0390624999999\displaystyle -3.35156250001 \, \eta^{7} + 3.35156250001 \, \eta^{6} + 3.86718750001 \, \eta^{5} - 3.86718750001 \, \eta^{4} - 1.0546875 \, \eta^{3} + 1.0546875 \, \eta^{2} + 0.0390624999999 \, \eta - 0.0390624999999
8.14072271825η77.09659483138η611.347477684η5+9.89205188146η4+3.33160871212η32.90429707348η20.124853746365η+0.108840023398\displaystyle 8.14072271825 \, \eta^{7} - 7.09659483138 \, \eta^{6} - 11.347477684 \, \eta^{5} + 9.89205188146 \, \eta^{4} + 3.33160871212 \, \eta^{3} - 2.90429707348 \, \eta^{2} - 0.124853746365 \, \eta + 0.108840023398
10.3581368289η7+6.12891144098η6+18.6833551584η511.054944637η48.6700371412η3+5.13006254948η2+0.34481881174η0.204029353468\displaystyle -10.3581368289 \, \eta^{7} + 6.12891144098 \, \eta^{6} + 18.6833551584 \, \eta^{5} - 11.054944637 \, \eta^{4} - 8.6700371412 \, \eta^{3} + 5.13006254948 \, \eta^{2} + 0.34481881174 \, \eta - 0.204029353468
11.3898137485η72.38387910961η624.032962502η5+5.03008025554η4+15.6735080469η33.280452976η23.0303592934η+0.63425183007\displaystyle 11.3898137485 \, \eta^{7} - 2.38387910961 \, \eta^{6} - 24.032962502 \, \eta^{5} + 5.03008025554 \, \eta^{4} + 15.6735080469 \, \eta^{3} - 3.280452976 \, \eta^{2} - 3.0303592934 \, \eta + 0.63425183007
11.3898137485η72.38387910961η6+24.032962502η5+5.03008025554η415.6735080469η33.280452976η2+3.0303592934η+0.63425183007\displaystyle -11.3898137485 \, \eta^{7} - 2.38387910961 \, \eta^{6} + 24.032962502 \, \eta^{5} + 5.03008025554 \, \eta^{4} - 15.6735080469 \, \eta^{3} - 3.280452976 \, \eta^{2} + 3.0303592934 \, \eta + 0.63425183007
10.3581368289η7+6.12891144098η618.6833551584η511.054944637η4+8.6700371412η3+5.13006254948η20.34481881174η0.204029353468\displaystyle 10.3581368289 \, \eta^{7} + 6.12891144098 \, \eta^{6} - 18.6833551584 \, \eta^{5} - 11.054944637 \, \eta^{4} + 8.6700371412 \, \eta^{3} + 5.13006254948 \, \eta^{2} - 0.34481881174 \, \eta - 0.204029353468
8.14072271825η77.09659483138η6+11.347477684η5+9.89205188146η43.33160871212η32.90429707348η2+0.124853746365η+0.108840023398\displaystyle -8.14072271825 \, \eta^{7} - 7.09659483138 \, \eta^{6} + 11.347477684 \, \eta^{5} + 9.89205188146 \, \eta^{4} - 3.33160871212 \, \eta^{3} - 2.90429707348 \, \eta^{2} + 0.124853746365 \, \eta + 0.108840023398
3.35156250001η7+3.35156250001η63.86718750001η53.86718750001η4+1.0546875η3+1.0546875η20.0390624999999η0.0390624999999\displaystyle 3.35156250001 \, \eta^{7} + 3.35156250001 \, \eta^{6} - 3.86718750001 \, \eta^{5} - 3.86718750001 \, \eta^{4} + 1.0546875 \, \eta^{3} + 1.0546875 \, \eta^{2} - 0.0390624999999 \, \eta - 0.0390624999999
L_0_xi_eta = L_0_xi * L_0_eta L_0_xi_eta.simplify_full()
(11.2329711915η711.2329711915η612.9611206055η5+12.9611206055η4+3.53485107423η33.53485107423η20.130920410156η+0.130920410156)ξ7+0.130920410156η7+(11.2329711915η7+11.2329711915η6+12.9611206055η512.9611206055η43.53485107423η3+3.53485107423η2+0.130920410156η0.130920410156)ξ60.130920410156η6+(12.9611206055η7+12.9611206055η6+14.9551391602η514.9551391602η44.07867431643η3+4.07867431643η2+0.151062011719η0.151062011719)ξ50.151062011719η5+(12.9611206055η712.9611206055η614.9551391602η5+14.9551391602η4+4.07867431643η34.07867431643η20.151062011719η+0.151062011719)ξ4+0.151062011719η4+(3.53485107423η73.53485107423η64.07867431643η5+4.07867431643η4+1.11236572266η31.11236572266η20.0411987304688η+0.0411987304688)ξ3+0.0411987304688η3+(3.53485107423η7+3.53485107423η6+4.07867431643η54.07867431643η41.11236572266η3+1.11236572266η2+0.0411987304688η0.0411987304688)ξ20.0411987304688η2+(0.130920410156η7+0.130920410156η6+0.151062011719η50.151062011719η40.0411987304688η3+0.0411987304688η2+0.00152587890624η0.00152587890624)ξ0.00152587890624η+0.00152587890624\displaystyle {\left(11.2329711915 \, \eta^{7} - 11.2329711915 \, \eta^{6} - 12.9611206055 \, \eta^{5} + 12.9611206055 \, \eta^{4} + 3.53485107423 \, \eta^{3} - 3.53485107423 \, \eta^{2} - 0.130920410156 \, \eta + 0.130920410156\right)} \xi^{7} + 0.130920410156 \, \eta^{7} + {\left(-11.2329711915 \, \eta^{7} + 11.2329711915 \, \eta^{6} + 12.9611206055 \, \eta^{5} - 12.9611206055 \, \eta^{4} - 3.53485107423 \, \eta^{3} + 3.53485107423 \, \eta^{2} + 0.130920410156 \, \eta - 0.130920410156\right)} \xi^{6} - 0.130920410156 \, \eta^{6} + {\left(-12.9611206055 \, \eta^{7} + 12.9611206055 \, \eta^{6} + 14.9551391602 \, \eta^{5} - 14.9551391602 \, \eta^{4} - 4.07867431643 \, \eta^{3} + 4.07867431643 \, \eta^{2} + 0.151062011719 \, \eta - 0.151062011719\right)} \xi^{5} - 0.151062011719 \, \eta^{5} + {\left(12.9611206055 \, \eta^{7} - 12.9611206055 \, \eta^{6} - 14.9551391602 \, \eta^{5} + 14.9551391602 \, \eta^{4} + 4.07867431643 \, \eta^{3} - 4.07867431643 \, \eta^{2} - 0.151062011719 \, \eta + 0.151062011719\right)} \xi^{4} + 0.151062011719 \, \eta^{4} + {\left(3.53485107423 \, \eta^{7} - 3.53485107423 \, \eta^{6} - 4.07867431643 \, \eta^{5} + 4.07867431643 \, \eta^{4} + 1.11236572266 \, \eta^{3} - 1.11236572266 \, \eta^{2} - 0.0411987304688 \, \eta + 0.0411987304688\right)} \xi^{3} + 0.0411987304688 \, \eta^{3} + {\left(-3.53485107423 \, \eta^{7} + 3.53485107423 \, \eta^{6} + 4.07867431643 \, \eta^{5} - 4.07867431643 \, \eta^{4} - 1.11236572266 \, \eta^{3} + 1.11236572266 \, \eta^{2} + 0.0411987304688 \, \eta - 0.0411987304688\right)} \xi^{2} - 0.0411987304688 \, \eta^{2} + {\left(-0.130920410156 \, \eta^{7} + 0.130920410156 \, \eta^{6} + 0.151062011719 \, \eta^{5} - 0.151062011719 \, \eta^{4} - 0.0411987304688 \, \eta^{3} + 0.0411987304688 \, \eta^{2} + 0.00152587890624 \, \eta - 0.00152587890624\right)} \xi - 0.00152587890624 \, \eta + 0.00152587890624