(D13)CSM$用户:[WHEREMAN.PROGRAMS.NPAINLEVE.SINGLE]P_SG.OUT;4(C14)批次(“np_exec.max”)$(C15)/*****************//*批处理文件NP_EXEC。最大*//* ************************************************************************* */exec_painleve(eq、alpha、do_resonances、max_resonance、do_simplifiation)$----------------------------------------------------------------方程的PAINLEVE分析,-2 F F+2 F F-F+F=0T X T X----------------------------------------------------------------阿尔法在原始方程中用U G代替f。0g的最小幂为[2α-2,3α,α]2α-2 2*G的系数是-2 UαG0吨X注:对于α=0,此术语无效,验证α=0是否会导致主导行为,如果确实如此,则使用此用户再次运行程序α的供应价值。亨斯,把贝塔系数设为0。3α3*G的系数IS-U0阿尔法*G的系数是U0----------------------------------------------------------------对于g的指数(2α-2)和(3α),当alpha=-2时,g的幂为-6---->求解U02个1术语U(4 G G-U)——占主导地位0 T X 0 6G公司等式中。----------------------------------------------------------------1)当U=4 G G时--->找到共振0吨X阿尔法R+阿尔法用U G+U G替换方程式中的f0卢比2 2右-6第(8(G)(G)项(R-2)(R+1))U G占主导地位T X R(T X R)等式中。唯一的非负整数根是[R=2]最大共振=2-->检查共振。2====\K-2将方程式中的f替换为>G U的幂级数。/K(K)====K=0U=4 G G0吨X1 2 2*--IS-16(G)(4 G+U)(G)的系数5吨X 1倍G公司U=-4克1吨X1*--的系数为04G公司你是任意的!2兼容条件已满足!----------------------------------------------------------------对于g的指数(2α-2)和(α),α=2时,g的幂为2-->求U02术语-U(4 U G G-1)G占主导地位0 0 T X等式中。----------------------------------------------------------------11)使用U=--------->找到共振0 4克T X(T X)阿尔法R+阿尔法用U G+U G替换方程式中的f0卢比R+2级术语((R-2)(R+1))U G占主导地位等式中。唯一的非负整数根是[R=2]最大共振=2-->检查共振。2====\K+2公司将方程式中的f替换为>G U的幂级数。/K(K)====K=01带U=-------0 4克T X(T X)2 24个U(G)(G)-G3 1吨X吨X*G的系数为------------------------2 24(G)(G)T X(T X)G公司T X(T X)U=-------------1 2 24(G)(G)T X(T X)4*G的系数为0你是任意的!2兼容条件已满足!----------------------------------------------------------------对于g的指数(3α)和(α),g的幂不是最小值--跳过此字母值。----------------------------------------------------------------(C16)输出()$----------------------------------------------------------------计算结束时,可以使用以下内容:*U值(uval[j,k,l]型,其中1<=j<=2和0<=k<=[2,2]和1<=l<=[1,1])j代表j_thα,k代表u[k],l代表u[0]的第l个解集*字母值(类型alpha[j],其中1<=j<=2)j代表j_thα*相容性条件(类型compcond[j,k],其中1<=j<=2和1<=k<=[1,1])j代表j_thα,k代表u[0]的k_th解集*共振(类型res[j,k],其中1<=j<=2和1<=k<=[1,1])j代表j_thα,k代表u[0]的k_th解集----------------------------------------------------------------要再次查看此菜单,只需键入<output()>----------------------------------------------------------------(C17)/****************NP_EXEC结束。最大***************************/(C18)关闭文件();