(D13)CSM$用户:[WHEREMAN.PROGRAMS.NPAINLEVE.SINGLE]P_RAM3.OUT;1(C14)批次(“np_exec.max”);(C15)/*****************//*批处理文件NP_EXEC。最大*//* ************************************************************************* */exec_painleve(eq、alpha、do_resonances、max_resonance、do_simplifiation)$替换X-->G+X0----------------------------------------------------------------2方程的PAINLEVE分析,F+F=0----------------------------------------------------------------阿尔法用U G代替原始方程中的f。0g的最小幂为[2α,α-1]2α2*G的系数是U0α-1*G是Uα的系数0注:对于α=0,此术语无效,验证α=0是否会导致主导行为,如果确实如此,则使用此用户再次运行程序α的供应价值。亨斯,把贝塔系数设为0。----------------------------------------------------------------对于g的指数(2α)和(α-1),当α=-1时,g的幂为-2-->求U01术语(U-1)U——占主导地位0 0 2等式中。----------------------------------------------------------------1)U=1-->找到共振0阿尔法R+阿尔法用U G+U G替换方程式中的f0卢比右-2术语(R+1)U G占主导地位R(右)等式中。r不存在非负整数根。----------------------------------------------------------------(C16)输出()$----------------------------------------------------------------计算结束时,可以使用以下内容:*U值(uval[j,k,l]型,其中1<=j<=1和0<=k<=[0]和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()>----------------------------------------------------------------(C17)/***************************NP_EXEC结束。最大********************/(D17)完成(C18)关闭文件();