{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 }{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 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 51 "macro(WHITE=COLOR(RGB,1,1,1),RED=COLOR(RGB,1,0 ,0)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "PlotTr:=proc(Verts ,color)\n POLYGONS(Verts,color,STYLE(PATCHNOGRID));\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "#PLOT(PlotTr([[1,1],[-1,1],[0,0]],R ED));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 488 "MedTr:=proc(Tr)\n local medtr,tries:\n medtr:=[[(Tr[1][1]+Tr[2][1])/2.,(Tr[1][2]+Tr[2][ 2])/2.],\n [(Tr[2][1]+Tr[3][1])/2.,(Tr[2][2]+Tr[3][2])/2.],\n \+ [(Tr[3][1]+Tr[1][1])/2.,(Tr[3][2]+Tr[1][2])/2.]]:\n tries:=[[T r[1][1],Tr[1][2]],[medtr[1][1],medtr[1][2]],[medtr[3][1],medtr[3][2]]] ,\n [[medtr[1][1],medtr[1][2]],[Tr[2][1],Tr[2][2]],[medtr[2][1] ,medtr[2][2]]],\n [[medtr[3][1],medtr[3][2]],[medtr[2][1],medtr [2][2]],[Tr[3][1],Tr[3][2]]];\n RETURN([tries,medtr]);\nend:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 120 "FStep:=proc(Tr,Tries,G)\n l ocal Res;\n Res:=MedTr(Tr);\n Tries:=Res[1..3]:\n G:=PlotTr(Res[ 4],WHITE);\n RETURN(0):\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "Tri:=[[1,1],[-1,0],[0,-1]];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "#Res:=MedTr(Tri);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "#display(PLOT(PlotTr(Res[4],WHITE)),PLOT(PlotTr(Tri,R ED)));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "g[1]:=PlotTr(Tri, RED):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "FStep(Tri,'Tries', 'g[2]'):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "display(PLOT(g[ 2],g[1]));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "N:=7:" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "GTries[1]:=[PlotTr(Tri,RED)]:" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "LTries[1]:=[Tri]:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 200 "for i from 2 to N do\n LTries[i]:=[];\n GTries[ i]:=[];\n for j from 1 to nops(LTries[i-1]) do\n FStep(LTries[i-1][j] ,'L','G');\n LTries[i]:=[op(LTries[i]),op(L)];\n GTries[i]:=[op(GTri es[i]),G];\n od:\nod:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "PLOT('op(G Tries[N+1-i])'$'i'=1..N):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "displa y(%,scaling=constrained);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "Pic:=NULL:\nfor k from 1 to N do\n Pic:=Pic,PLOT('op(GTries[k+1-i] )'$'i'=1..k):\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "displ ay([Pic],insequence=true,axes=NONE);" }}}{EXCHG }{EXCHG }{EXCHG } {EXCHG }}{MARK "20" 0 }{VIEWOPTS 1 1 1 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }