(D18)CSM$用户:[鉴于。程序。NPAINLEVE.SINGLE]P_KDV。OUT(输出);(C19)批次(“np_exec.max”);(C20)/*****************//*批处理文件NP_EXEC。最大*//* ************************************************************************* */exec_painleve(eq、alpha、do_resonances、max_resonance、do_simplifiation)$----------------------------------------------------------------方程的PAINLEVE分析,F+F F+F=0X X X X T----------------------------------------------------------------阿尔法在原始方程中用U G代替f。0g的最小幂为[2α-1,α-3]2α-1 2*G的系数是UαG0倍注:对于α=0,此术语无效,验证α=0是否会导致主导行为,如果确实如此,则与该用户一起再次运行该程序α的供应价值。亨斯,把贝塔系数设为0。α-3 3*G的系数是U(α-2)(α-1)α(G)0倍注:此术语不适用于[α=0,α=1,α=2],验证阿尔法的任何这些值是否会导致占主导地位行为,如果确实如此,则使用此值再次运行程序作为用户提供的ALPHA,称为BETA。----------------------------------------------------------------对于g的指数(2α-1)和(α-3),当α=-2时,g的幂为-5-->求U02 1学期-2 U G(12(G)+U)--占主导地位0 X X 0 5G公司等式中。----------------------------------------------------------------21)U=-12(G)-->找到共振0倍阿尔法R+阿尔法用U G+U G替换方程式中的f0卢比3转-5术语((G)(R-6)(R-4)(R+1))U G占主导地位X R(X R)在方程式中。2个非负整数根是[R=4,R=6]最大共振=6-->检查共振。6====\K-2将方程中的f替换为>G U的幂级数。/K(K)====K=02U=-12(G)0倍1 2 2*系数--IS 6(G)((-12(G))-36 G G+5 U G)4 X X X X X 1 XG X公司U=12克1 X X1 2 2*--的系数为24 G(4 G G-3(G)+U(G)3 X X X X X X X X 2 XG公司+G G)T X(T X)24 G G-3(G)+G GX X X X X X X X T XU=--------------------------------2 2(G)X1 2*--IS-12((G)G-4 G G的系数2 X X X X X X X X X X X XG公司3 4 2+3(G)-G G-U(G)+G(G))/GX X T X X X X 3 X T X X X X XU=2 3 2(G)G-4 G G G+3(G)-G G+G(G)X X X X X X X X X-X X X×X X X T X X X-----------------------------------------------------------------------4(G)X1*系数-IS 0G公司U是任意的!4兼容条件已满足!4 3*1的系数是-(G)G-9(G)GX X X X X X X X X-X X X XX X X3 2 2 3-17(G)G G+48(G)(G)G-2 G(G)X X X X X X X X X-X X X×X X X T X X X2 2 3 2+70(G)G(G)-174(G)G+17(G)X X X X X X X X X-X X X T X X X×X X X3 5 3 2 2-8G(G)G+81(G)-21G(GT X X X X X X X X T X X T X X X6 3 2 2 8+6个U(G)G-9个G(G)G+(G)G+6个U(G)4 X X X T X X X X X X T X X X X X 5 X7 4 4 3 5+6个U(G)+G(G)+2个G(G4 X T T X T X X X X T X X X X XX4 3U=-((G)G-9(G)GG5 X X X X X X X X X X X X X X X X X X X X X X X3 2 2 3-17(G)G G+48(G)(G)G-2 G(G)X X X X X X X X X-X X X×X X X T X X X2 2 3 2+70(G)G(G)-174(G)G+17(G)X X X X X X X X X-X X X T X X X×X X X3 5 3 2 2-8G(G)G+81(G)-21G(GT X X X X X X X X T X X T X X X6 3 2 2 7+6个U(G)G-9个G(G)G+(G)G+6个U(G)4 X X X T X X X X X X T X X X X X 4 XX4 4 3 8+G(G)+2 G(GT T X T X X X X T X X X X X*G的系数为0你是任意的!6兼容条件已满足!----------------------------------------------------------------(C21)输出()$----------------------------------------------------------------计算结束时,可以使用以下内容:*U值(uval[j,k,l]型,其中1<=j<=1和0<=k<=[6]和1<=l<=[1])j代表j_thα,k代表u[k],l代表u[0]的第l个解集*字母值(类型alpha[j],其中1<=j<=1)j代表j_thα*相容性条件(类型compcond[j,k],其中1<=j<=1和1<=k<=[1])j代表j_thα,k代表u[0]的k_th解集*共振(类型res[j,k],其中1<=j<=1和1<=k<=[1])j代表j_thα,k代表u[0]的k_th解集----------------------------------------------------------------要再次查看此菜单,只需键入<output()>----------------------------------------------------------------(C22)/****************NP_EXEC结束。最大********************/(D22)完成(C23)关闭文件();