{VERSION 6 0 "IBM INTEL NT" "6.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 25 "restart: with(algcurves):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 264 "The puiseux program can give unnecessarily lar ge answers. Reading the following code into Maple before using the pui seux command will help to prevent one (but not all) of the causes of t hese large answers, namely it will help to prevent some unnecessary ex panding." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 2224 "`algcurves/li ft_exp` := proc(v, f, x, y)\nlocal i, ii, r, res, v7, vv7, v3, ext, a, j, n, np, ram, j3;\n if v[5] = 1 then return \{v\} end if;\n v3 := degree(v[3], x);\n res := \{\};\n r := v[1] + y*x^v[2];\n \+ vv7 := v[7]*v3 + v[2] - 1; vv7 := vv7 + v[5];\n ii := `algcurve s/truncate_subs`(subs(x = v[3], f), x, y, r, vv7 + 1,\n v[4]); \n if ii = 0 then error \"degree estimate was wrong\" end if;\n \+ v7 := (ldegree(ii, x) - v[2])/v3;\n r := `algcurves/v_ext_m`(\n \+ `algcurves/g_factors`(tcoeff(ii, x), y, v[4]), y);\n for i in r do res := res union `algcurves/lift_exp`([\n v[1] + x^v[2]*i[1 ], v[2] + 1, v[3], [op(i[3]), op(v[4])], i[2],\n v[6]*i[4], v7, [op(v[8]), [op(1 .. 4, v)]]], f, x, y)\n end do;\n if add(i[5]* i[6]*degree(i[3], x)/(v[6]*v3), i = res) <>\n degree(tcoeff(ii, x), y) then error \"found wrong number of expansions\"\n end if;\n \+ if v[5] = degree(tcoeff(ii, x), y) then\n if ldegree(ii, x) <> \+ vv7 then error \"degree estimate was wrong\"\n end if;\n \+ return res\n end if;\n ii := collect(ii, y);\n ii := add(`al gcurves/normal_tcoeff`(coeff(ii, y, i), x)*y^i,\n i = 0 .. degr ee(ii, y));\n np := `algcurves/Newtonpolygon`(ii, x, y);\n if no ps(np) = 2 and np[1][3] = 0 then\n error \"found wrong number o f expansions\"\n end if;\n for j in np do\n if 2 < nops(j ) and 0 < j[3] and j[3] < 1 then\n r := `algcurves/g_factor s`(j[4], x, v[4]);\n r := `algcurves/v_ext_m`(r, x);\n \+ for i in r do\n j3 := j[3] - v[2];\n \+ ext := [op(i[3]), op(v[4])];\n n := mods(1/numer(j3 ), denom(j3));\n ram := i[1]^n*x^denom(j3);\n \+ a := v[2]*denom(j3) - numer(j[3]);\n res := res u nion `algcurves/lift_exp`([collect(\n subs(x = ram, v[1])\n + x^a*i[1]^((1 - n*numer(j3))/denom(j3)), x, normal),\n a + 1, normal(subs(x = ram, v[3])), \+ ext, i[2],\n v[6]*i[4],\n (j[2] \+ - j[1]*j[3] - a/degree(ram, x))/degree(v[3], x),\n \+ [op(v[8]), [op(1 .. 4, v)]]], f, x, y)\n end do\n en d if\n end do;\n res\nend proc:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 107 "A:=y^2 +2*y*x^2+x^4+0*y^2*x+0*x^3*y+b*x^2*y^2+d*x*y^3+e*y^4+a*x^4*y+h*x^3*y^2 +j*x^2*y^3+k*x*y^4+l*y^5+c*y^3;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% \"AG,:*$)%\"yG\"\"#\"\"\"F**(F)F*F(F*)%\"xGF)F*F**$)F-\"\"%F*F**(%\"bG F*F,F*F'F*F**(%\"dGF*F-F*)F(\"\"$F*F**&%\"eGF*)F(F0F*F**(%\"aGF*F/F*F( F*F**(%\"hGF*)F-F6F*F'F*F**(%\"jGF*F,F*F5F*F**(%\"kGF*F-F*F9F*F**&%\"l GF*)F(\"\"&F*F**&%\"cGF*F5F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 21 "Define polynomial [1]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "p uiseux(A,x=0,y,0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$-%'RootOfG6#,* \"\"\"F(*&%\"eGF()%#_ZG\"\"#F(F(*&%\"lGF()F,\"\"$F(F(*&%\"cGF(F,F(F(,& *$)%\"xGF-F(!\"\"*&)F7F1F(-F%6#,*%\"aGF8*$F+F(F(%\"bGF(F3F8F(F(" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 13 "discrim=0 [4]" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 18 "A1:=eval(A,c=b-a);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A1G,:*$)%\"yG\"\"#\"\"\"F**(F)F*F(F*)%\"xGF)F*F**$)F -\"\"%F*F**(%\"bGF*F,F*F'F*F**(%\"dGF*F-F*)F(\"\"$F*F**&%\"eGF*)F(F0F* F**(%\"aGF*F/F*F(F*F**(%\"hGF*)F-F6F*F'F*F**(%\"jGF*F,F*F5F*F**(%\"kGF *F-F*F9F*F**&%\"lGF*)F(\"\"&F*F**&,&F2F*F;!\"\"F*F5F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "reevaluate polynomial [3]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "puiseux(A1,x=0,y,0);" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#<$-%'RootOfG6#,*\"\"\"F(*&%\"lGF()%#_ZG\"\"$F(F( *&%\"eGF()F,\"\"#F(F(*&,&%\"bGF(%\"aG!\"\"F(F,F(F(,&*$)%\"xGF1F(F6*&)* &F:F(,&%\"hGF6%\"dGF(F6#\"\"(F1F()F>\"\"%F(F(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 15 "numerator=0 [2]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "A2:=eval(A1,h=d);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A2G,:* $)%\"yG\"\"#\"\"\"F**(F)F*F(F*)%\"xGF)F*F**$)F-\"\"%F*F**(%\"bGF*F,F*F 'F*F**(%\"dGF*F-F*)F(\"\"$F*F**&%\"eGF*)F(F0F*F**(%\"aGF*F/F*F(F*F**(F 4F*)F-F6F*F'F*F**(%\"jGF*F,F*F5F*F**(%\"kGF*F-F*F9F*F**&%\"lGF*)F(\"\" &F*F**&,&F2F*F;!\"\"F*F5F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "r eevaluate polynomial [3]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "puiseux(A2,x=0,y,0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$-%'RootOfG6 #,*\"\"\"F(*&%\"lGF()%#_ZG\"\"$F(F(*&%\"eGF()F,\"\"#F(F(*&,&%\"bGF(%\" aG!\"\"F(F,F(F(,&*$)%\"xGF1F(F6*&)F:\"\"%F(-F%6#,**$F0F(F(*&,&F4F(*&F1 F(F5F(F6F(F,F(F(F/F(%\"jGF6F(F(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 13 "discrim=0 [4]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "expand ((b-2*a)^2-4*(e-j));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,,*$)%\"bG\"\" #\"\"\"F(*(\"\"%F(F&F(%\"aGF(!\"\"*&F*F()F+F'F(F(*&F*F(%\"eGF(F,*&F*F( %\"jGF(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "A3:=eval(A2,e= 1/4*(b^2-4*b*a+4*a^2+4*j));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A3G, :*$)%\"yG\"\"#\"\"\"F**(F)F*F(F*)%\"xGF)F*F**$)F-\"\"%F*F**(%\"bGF*F,F *F'F*F**(%\"dGF*F-F*)F(\"\"$F*F**&,**&#F*F0F**$)F2F)F*F*F**&F2F*%\"aGF *!\"\"*$)F>F)F*F*%\"jGF*F*)F(F0F*F**(F>F*F/F*F(F*F**(F4F*)F-F6F*F'F*F* *(FBF*F,F*F5F*F**(%\"kGF*F-F*FCF*F**&%\"lGF*)F(\"\"&F*F**&,&F2F*F>F?F* F5F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "reevaluate polynomial [ 3]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "puiseux(A3,x=0,y,0); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#<$,(*(\"#K!\"\",(*&%\"dG\"\"\"%\"b GF+F'*(\"\"#F+F*F+%\"aGF+F+*&F.F+%\"kGF+F+\"\"&*&%\"xGF+,(*(F.F'F*F+F, F+F+*&F*F+F/F+F'F1F'F'#\"\"*F.F'*,F&F'F(\"\"%,&F,F'*&F.F+F/F+F+F+F4F;F 5!\"%F+**F;F'F(F.F4F.F5!\"#F'-%'RootOfG6#,*F;F+*(F;F+%\"lGF+)%#_ZG\"\" $F+F+*&,**$)F,F.F+F+*(F;F+F,F+F/F+F'*&F;F+)F/F.F+F+*&F;F+%\"jGF+F+F+)F HF.F+F+*&,&*&F;F+F,F+F+*&F;F+F/F+F'F+FHF+F+" }}}{EXCHG {PARA 0 "" 0 " " {TEXT -1 15 "numerator=0 [2]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "A4:=eval(A3,k=1/2*d*b-d*a);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%#A4G,:*$)%\"yG\"\"#\"\"\"F**(F)F*F(F*)%\"xGF)F*F**$)F-\"\"%F*F**(% \"bGF*F,F*F'F*F**(%\"dGF*F-F*)F(\"\"$F*F**&,**&#F*F0F**$)F2F)F*F*F**&F 2F*%\"aGF*!\"\"*$)F>F)F*F*%\"jGF*F*)F(F0F*F**(F>F*F/F*F(F*F**(F4F*)F-F 6F*F'F*F**(FBF*F,F*F5F*F**(,&*&#F*F)F**&F4F*F2F*F*F**&F4F*F>F*F?F*F-F* FCF*F**&%\"lGF*)F(\"\"&F*F**&,&F2F*F>F?F*F5F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 26 "reevalutate polynomail [3]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "puiseux(A4,x=0,y,0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$,(*$)%\"xG\"\"#\"\"\"!\"\"*&)F'\"\"%F),&%\"aGF)*&F(F* %\"bGF)F*F)F)*&)F'\"\"&F)-%'RootOfG6#,2*&)F1F(F)F/F)F**(F-F)F1F))F/F(F )F)*&F-F))F/\"\"$F)F**(F-F)%\"jGF)F/F)F**(F(F)FAF)F1F)F)*&F-F)%\"lGF)F **&F-F))%#_ZGF(F)F)*(F-F)%\"dGF)FGF)F)F)F)-F66#,*F-F)*(F-F)FDF))FGF?F) F)*&,**$F:F)F)*(F-F)F1F)F/F)F**&F-F)F " 0 "" {MPLTEXT 1 0 58 "expand((4*d)^2-16*(-b^2*a+4*b*a^2-4*a^3+2*j*b-4*j*a-4*l));" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# ,0*&\"#;\"\"\")%\"dG\"\"#F&F&*(F%F&)%\"bGF)F&%\"aGF&F&*(\"#kF&F,F&)F-F )F&!\"\"*&F/F&)F-\"\"$F&F&*(F/F&%\"jGF&F-F&F&*(\"#KF&F6F&F,F&F1*&F/F&% \"lGF&F&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 13 "discrim=0 [4]" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "A5:=eval(A4,l=-1/64*(16*d^2+ 16*b^2*a-64*b*a^2+64*a^3-32*j*b+64*j*a));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A5G,:*$)%\"yG\"\"#\"\"\"F**(F)F*F(F*)%\"xGF)F*F**$)F -\"\"%F*F**(%\"bGF*F,F*F'F*F**(%\"dGF*F-F*)F(\"\"$F*F**&,**&#F*F0F**$) F2F)F*F*F**&F2F*%\"aGF*!\"\"*$)F>F)F*F*%\"jGF*F*)F(F0F*F**(F>F*F/F*F(F *F**(F4F*)F-F6F*F'F*F**(FBF*F,F*F5F*F**(,&*&#F*F)F**&F4F*F2F*F*F**&F4F *F>F*F?F*F-F*FCF*F**&,.*&F0F?F4F)F?*(F0F?F2F)F>F*F?*&F2F*FAF*F**$)F>F6 F*F?*(F)F?FBF*F2F*F**&FBF*F>F*F?F*)F(\"\"&F*F**&,&F2F*F>F?F*F5F*F*" }} }{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "reevaluate polynomial [3]" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "puiseux(A5,x=0,y,0);" }} {PARA 12 "" 1 "" {XPPMATH 20 "6#<$-%'RootOfG6#,*\"\"%!\"\"*&,.*$)%\"dG \"\"#\"\"\"F0*&)%\"bGF/F0%\"aGF0F0*(F(F0F3F0)F4F/F0F)*&F(F0)F4\"\"$F0F 0*(F/F0%\"jGF0F3F0F)*(F(F0F;F0F4F0F0F0)%#_ZGF9F0F0*&,**&F(F0F;F0F)*$F2 F0F)*(F(F0F3F0F4F0F0*&F(F0F6F0F)F0)F>F/F0F0*&,&*&F(F0F3F0F)*&F(F0F4F0F 0F0F>F0F0,***\"#kF)F.\"\"',(*&F/F0F6F0F0*&F3F0F4F0F)F;F0FM,$**F/F0%\"x GF0F.F)FNF)F)#\"#6F/F0*(F/F)FS\"\"&F.F0F)*(F/F),&F3F)*&F/F0F4F0F0F0FSF (F0*$)FSF/F0F)" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 15 "numerator=0 [2] " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "A6:=eval(A5,j=b*a-2*a^2 );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%#A6G,:*$)%\"yG\"\"#\"\"\"F**(F )F*F(F*)%\"xGF)F*F**$)F-\"\"%F*F**(%\"bGF*F,F*F'F*F**(%\"dGF*F-F*)F(\" \"$F*F**&,&*&F0!\"\"F2F)F**$)%\"aGF)F*F:F*)F(F0F*F**(F=F*F/F*F(F*F**(F 4F*)F-F6F*F'F*F**(,&*&F2F*F=F*F**&F)F*FF*F**&,.*&F0F:F4F)F:*(F0F:F2F)F=F*F:* &F2F*F " 0 "" {MPLTEXT 1 0 20 "puiseux(A6,x =0,y,0);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#<$-%'RootOfG6#,*\"\"%!\"\" *&,**$)%\"dG\"\"#\"\"\"F0*&)%\"bGF/F0%\"aGF0F)*(F(F0F3F0)F4F/F0F0*&F(F 0)F4\"\"$F0F)F0)%#_ZGF9F0F0*&,&*$F2F0F)*&F(F0F6F0F0F0)F;F/F0F0*&,&*&F( F0F3F0F)*&F(F0F4F0F0F0F;F0F0,**$)%\"xGF/F0F)*&)FHF(F0,&F4F0*&F/F)F3F0F )F0F0*(F/F)FH\"\"&F.F0F)*&)FH\"\"'F0-F%6#,4*&\"\")F0F@F0F0*&,(*(\"#KF0 F3F0F4F0F)*&FZF0F6F0F0*&FVF0F2F0F0F0F;F0F0*&F3F0F-F0F0*(\"#kF0F3F0F8F0 F)*(\"#[F0F2F0F6F0F0*(\"#;F0)F3F9F0F4F0F)*(F(F0F4F0F-F0F)*&F/F0)F3F(F0 F0*&FZF0)F4F(F0F0F0F0" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 13 "discrim= 0 [4]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "discrim(8*_Z^2+(- 32*b*a+32*a^2+8*b^2)*_Z+b*d^2-64*b*a^3+48*b^2*a^2-16*b^3*a-4*a*d^2+2*b ^4+32*a^4,_Z);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&*(\"#K\"\"\"%\"bGF &)%\"dG\"\"#F&!\"\"*(\"$G\"F&%\"aGF&F(F&F&" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 29 "solve(-32*b*d^2+128*a*d^2,a);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$*&\"\"%!\"\"%\"bG\"\"\"F(" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 19 "A7:=eval(A6,a=b/4);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A7G,:*$)%\"yG\"\"#\"\"\"F**(F)F*F(F*)%\"xGF)F*F**$)F-\"\"%F*F **(%\"bGF*F,F*F'F*F**(%\"dGF*F-F*)F(\"\"$F*F***F6F*\"#;!\"\"F(F0F2F)F* **F0F9F2F*F-F0F(F*F**(F4F*)F-F6F*F'F*F***\"\")F9F2F)F-F)F(F6F**,F0F9F- F*F(F0F4F*F2F*F**&,&*&F0F9F4F)F9*&\"#kF9F2F6F*F*)F(\"\"&F*F***F6F*F0F9 F(F6F2F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "reevaluate polynomi al [3]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "puiseux(A7,x=0,y, 0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$,,*$)%\"xG\"\"#\"\"\"!\"\"*( \"\"%F*F'F,%\"bGF)F**(F(F*F'\"\"&%\"dGF)F**(\"\")F*F'\"\"'F-F(F**(\"(_ r4#F*,$*(F2F)F'F)F0!\"$F*#\"#8F(F0\"#@F*-%'RootOfG6#,*\"#kF**&,&*&\"#; F))F0F(F)F)*$)F-\"\"$F)F*F))%#_ZGFHF)F)*(\"#7F))FJF(F))F-F(F)F**(\"#[F )F-F)FJF)F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 15 "numerator=0 [2]" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "A8:=eval(A7,d=0);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A8G,4*$)%\"yG\"\"#\"\"\"F**(F)F*F(F *)%\"xGF)F*F**$)F-\"\"%F*F**(%\"bGF*F,F*F'F*F**&#\"\"$\"#;F**&)F(F0F*) F2F)F*F*F**&#F*F0F**(F2F*F/F*F(F*F*F**&#F*\"\")F**(F9F*F,F*)F(F5F*F*F* *&#F*\"#kF**&)F(\"\"&F*)F2F5F*F*F**&#F5F0F**&FAF*F2F*F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "reevaluate polynomial [3]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "factor(A8);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$*(\"#k!\"\",&\"\"%\"\"\"*&%\"bGF)%\"yGF)F)F),(*&F(F)F ,F)F)*&F+F))F,\"\"#F)F)*&F(F))%\"xGF1F)F)F1F)" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 22 "Polynomial factors [5]" }}}}{MARK "44 0 0" 18 } {VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }