{VERSION 4 0 "IBM INTEL LINUX22" "4.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 "" 0 21 "" 0 1 0 0 0 1 0 0 0 0 2 0 0 0 0 1 } {CSTYLE "" -1 256 "" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 0 }{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 "Heading 1" 0 3 1 {CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }1 0 0 0 8 4 0 0 0 0 0 0 -1 0 } {PSTYLE "Heading 2" -1 4 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 2 1 0 1 0 2 2 0 1 }{PSTYLE "Title" -1 18 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 1 2 2 2 1 1 1 1 }3 1 0 0 12 12 1 0 1 0 2 2 19 1 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Co urier" 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 "R3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Times" 1 12 128 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 "Normal" -1 258 1 {CSTYLE "" -1 -1 "Times" 1 16 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Title" -1 259 1 {CSTYLE "" -1 -1 "Ti mes" 1 24 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }3 1 0 0 12 12 3 0 3 0 2 2 19 1 }{PSTYLE "Maple Output" -1 260 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 8 "restart;" }}}{EXCHG {PARA 18 "" 0 "" {TEXT 256 10 "PLANX - 08" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 350 "dir:= \"C:/MapleV6/Users/Douillet/\";\ndir:= \"/ho me/douillet/docs/Ensait/planx/\";\n# ksv is convert \" \" -> \"_\" and \",\" -> \".\"\n# nam, ix, ksv:= \"experimentique.txt\", 5, true;\n# \+ nam, ix, ksv:= \"nmt-edu.txt\", 8, false;\n# nam, ix, ksv:= \"nist-cer amic.txt\", 5, false;\nnam, ix, ksv:= \"nist-ceramic.cut.txt\", 4, fal se;\nnam, ix, ksv:= \"nist-light.txt\", 7, false;" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 15 "Initialisations" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "kernelopts(ASSERT=true);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 46 "with(pldx); ever_loaded; with(plots, display):" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "with(simul); ever_loaded; with(Line arAlgebra); rand(); xmacro();" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "fich:= cat(dir, nam); fclose(%):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "interface(rtablesize=20);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 80 "(cut4):= proc(x)\nif type(x,numeric) then Float(rou nd(x*10000),-4)\nelse x fi end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 7 "Lecture" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 385 "kill_space s_virg:= proc(txt, ksv) local line, ouca; \nif ksv = false then return txt fi;\nline:= txt: do ouca:= searchtext(\" \", line);\nif ouca=0 th en break fi;\nline:= cat(substring(line,1..ouca-1), \"_\", substring(l ine,ouca+1..-1));\nod;\ndo ouca:= searchtext(\",\", line);\nif ouca=0 \+ then break fi;\nline:= cat(substring(line,1..ouca-1), \".\", substring (line,ouca+1..-1));\nod;\nreturn line; end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 312 "fclose(fich): hnd:= fopen(fich,READ):\nfrm:= cat( \"%s\"$(ix+1)): line:= readline(hnd): kill_spaces_virg(%, ksv): titres := sscanf(%, frm);\nfrm:= cat(\"%s\"$(ix), \"%f\"):\nj:=1: do\n line:= readline(hnd); if line=0 then break fi;\n kill_spaces_virg(line, ksv) :\n data[j]:= sscanf(%, frm); j:= j+1;\nod: jx:=j-1; fclose(fich);" }} }{EXCHG {PARA 0 "" 0 "" {MPLTEXT 0 21 245 "frm:= cat(\"%s\\t\"$(ix-1), \"%s\\n\");\nmm:= :\nmmm:= Matrix([Column(mm,[2..8, 1])]);\nhnd:= fopen(fich,WRITE ):\nfor j to jx+1 do fprintf(hnd, frm, op(convert(Row(mmm,j), list))); od:\nfclose(fich):" }}}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 42 "Lecture s peciale de \"nist-ceramic.cut.txt\"" }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 42 "u est le num\351ro du param\350tre \340 ne pas lire" }}{PARA 0 "" 0 "" {MPLTEXT 0 21 478 "ix:=4: u:=1: valu:=\"1\"; notu:= subs(u=NUL L, [$1..ix+1]);\nfclose(fich): hnd:= fopen(fich,READ):\nfrm:= cat(\"%s \"$(ix+1)): line:= readline(hnd): kill_spaces_virg(%, ksv): \ntmp:= ss canf(%, frm): titres:= [seq(tmp[k], k=notu)];\nfrm:= cat(\"%s\"$(ix), \+ \"%f\"):\nj:=1: do\n line:= readline(hnd); if line=0 then break fi;\n \+ kill_spaces_virg(line, ksv):\n tmp:= sscanf(%, frm); if tmp[u]=valu th en \n data[j]:= [seq(tmp[k], k=notu)]; j:= j+1; fi; \nod: jx:=j-1; i x:= nops(titres)-1; fclose(fich);" }}}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 19 "Modalit\351s et codage" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 179 "zone:= table( ): zone[0]:=1..1: deb:=2:\nfor m to ix do\nmodalite[m]:= sort(convert( \{seq(data[j][m], j=1..jx)\}, list));\nmaxm[m]:= nops(%);\nzone[m]:= d eb..deb+%-2; deb:= rhs(%)+1:\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 322 "coda:= table():\nmkcoda:= proc(moda) local lesmod, n, nn;\nglobal coda; \nle smod:= modalite[moda]; n:= nops(%)-1;\ncoda[moda]:= proc() \n error s printf(\"code inattendu, moda=%a, valeur=%a\", moda, args) \nend;\ncod a[moda](lesmod[1]):= (-1)$ n;\nfor nn from 1 to n do\ncoda[moda](lesmo d[nn+1]):= 0$ (nn-1),1,0$(n-nn);\nod: NULL;\nend:" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 76 "seq(mkcoda(m), m=1..ix); # xprint(coda[4]);x print(coda[7]); xprint(coda[1]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "codeur: = proc(dataline)[1,seq(coda[m](dataline[m]), m=1..ix)] end;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "titres; data[1];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 83 "mma, mmb:= Matrix([seq(codeur(data[ j]), j=1..jx)]),\n ;" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 0 "" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 15 "Moindres carr\351s" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "ma, mb:= eval(mma), eval(mmb );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "1/(Transpose(ma).ma). Transpose(ma).mb: mf:= evalf(%);" }}}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 3 "FRV" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "jjx,iix:= Dimensio n(ma); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 169 "mr:= ma.mf-mb: \+ [moy, var](convert(%,list)); var_res_exp:= %[2];\nlib:= convert(evalf( mb), list): moy_b, var_b:= (moy, var)(%);\n(op@xab)(lib): %[2]-%[1], \+ sqrt(var(lib)); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "evalf(T r(mb).mb/jjx)= moy_b^2+ var_b;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "Tr(mr).mr/jjx; " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 " frv_exp:= var_b/var_res_exp; sqrt(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "var_res_th:= var_res_exp*jjx/(jjx-iix); frv_th:= var_ b/%;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "sqrt(var_res_exp), \+ sqrt(var_res_th);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 167 "pl1:= plots[pointplot]([seq ([mb[j], mr[j]], j=1..jx-1)]):\nrerange:= proc(quoi) local a,b;\n a,b: = op(1,quoi), op(2,quoi); (a+b)/2+(a-b)/2*1.1..(a+b)/2+(b-a)/2*1.1; en d:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 155 "[([min, max]@op@ con vert)(mb, list), ([min, max]@op@ convert)(mr, list)]:\nevalf(%,2); dis playg(pl1, tickmarks=%, labels=[exper, ``], view=map(rerange, %));" }} }{EXCHG {PARA 0 "" 0 "" {MPLTEXT 0 21 157 "cestac:= < ('2*ra()-1'$16) \+ >:\n1/(Transpose(ma).ma).Transpose(ma).(mb+cestac): mff:= evalf(%):\n[ seq(1-mff[j]/mf[j], j=1..12)]:evalf(%,3);\nma.(mff-mf), cestac;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 10 "Influences " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "convert(mf,list):\nli:= \+ map(Z-> Float(round(100*Z),-2), %);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "seq(li[zone[i]], i=0..ix);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 87 "cave r:= : \{Entries(evalf(ma.caver))\};\nASSERT(nops(%) =1); r0:= op(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 236 "for i \+ to ix do ; zoni:= zone[i];\n# modalite[i], mf[zoni]; \nc||(titres[i]): = Vector([(0$(lhs(zoni)-1)), mf[zoni],(0$(iix-rhs(zoni)))]): \nr||(tit res[i]):= evalf(ma.%);\nsort(convert(\{Entries(%)\}, list)); print( ca t(` r`, titres[i])=%);\nod:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 26 "Co mpensation des d\351calages" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 115 "[se q( add(k, k=convert(r||(titres[i]), list)), i=1..ix)];\nadd(k, k=%); ( moy(convert(mb, list)) - r0)*Dimension(mb);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 298 "graphe:= NULL; x:= 1; # x est l'abscisse du point \nfor i to ix do\ntmp:= convert(mf[zone[i]], list); how:= nops(%);\n[ [x, -add(k, k=tmp)], seq([x+k, tmp[k]], k=1..how)];\ngraphe:= graphe, \+ plots[pointplot](%, style='line'), plots[pointplot](%);\nx:=x+how+2 \n od: graphe:= graphe, plot(0, y=0..x): x:='x':" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "zone[1];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 252 "xticks:= [seq( -3+(lhs(zone[i])+rhs(zone[i])+1.)/2+2*i=titres[i ], i=1..ix)]:\nselect(has, map(op, [graphe]), POINTS): \nyticks:= eval f( ([min, max]@op@map2)(op, 2, map(op, %)) ,2):\ndisplayg(graphe, axes =boxed, tickmarks=[xticks, yticks]); influences:= %:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 187 "xres_ini(ix, [` #`, ` min`, ` max`, ` de lta`, `sd`]):\nxres_run(proc(i) convert(r||(titres[i]), list); evalf(x ab(%)): \n[titres[i], op(%), rhs(%)-lhs(%), (sqrt@var)(%%)] end): map (cut4,%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 234 "optimal:= tab le(): for i to ix do\n convert(mf[zone[i]], list); [-add(k,k=%), op(%) ];\n member((max@op)(%), %, 'ou'); optimal[i]:= modalite[i][ou];\nod: \+ \noptimal:= [seq(optimal[i], i=1..ix)];\ncodeur(optimal); evalf(Vector [row](%).mf, 6);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 112 "conver t(mb, list): member((max@op)(%), %, 'ou'):\nchallenger:= data[ou]; cod eur(%); evalf(Vector[row](%).mf, 6);" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 14 "Visualisatio ns" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "item:= proc(a,b) local ou; member(a,b,'ou'); ou-1; end;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "data[2]; seq( item(%[i], modalite[i]), i=1..ix);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 231 "[$1..ix]; grx,gry:= select( type,%,odd), select(type,%,even);\nif nam=\"nist-ceramic.txt\" then gr x, gry:=[1],[3,4,5,2]; fi;\nxxx, yyy:= mul(maxm[k], k=grx), mul(maxm[k ], k=gry);\ninterface(rtablesize=max(xxx,yyy,interface(rtablesize))); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 221 "tmpx:=0: for k from no ps(grx) by -1 to 1 do tmpx:= tmpx*maxm[grx[k]]+x[grx[k]] od:\ntmpy:=0: for k from nops(gry) by -1 to 1 do tmpy:= tmpy*maxm[gry[k]]+x[gry[k]] od:\nmarqueur:= unapply([tmpx+1,tmpy+1], seq(x[i],i=1..ix));" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "plan:= Matrix(xxx,yyy, (a,b) ->0): " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 153 "for j from 1 to jx do da ta[j]; seq( item(%[i], modalite[i]), i=1..ix);\n(op@marqueur)(%); plan [%]:=plan[%]+1; od: subs(0=`.`, 1=`*`,plan); # Transpose(%);" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 15 "Valeurs propres" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "ms:= Matrix(Transpose(ma).ma, shape=symmetric);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 403 "if nam=\"nmt-edu.txt\" then \n mv,mp:= E igenvectors(ms): [Column(mp,1..-1)]: sort(%, \n proc(a,b) local u: if \+ a[1]<>b[1] then return evalb(a[1]>b[1]) fi ;\n if a[2]<>b[2] then retu rn evalb(a[2]>b[2]) fi ;\n for u from 3 by 2 to iix do \n if a[u]<>b[u ] then return evalb(abs(a[u])>abs(b[u])) fi ;\n od; true end): mp:= Ma trix(%);\n 1/mp.ms.mp; mv:= Vector([seq(%[j,j], j=1..iix)]); \n print( 'mp, mv'=mp, mv); \nfi: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 " " }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 1 {PARA 3 " " 0 "" {TEXT -1 15 "Special ceramic" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "procol: = proc(u,v); zip((a,b)->a*b, Column(mA,u), Column(mA,v)) end:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "lots:=[[2,3],[2,4],[3,4]],[[ 5,4]];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 165 "Titres:= Cte, op (titres[1..ix]): \nTitres:= %, seq( cat(seq(Titres[j], j=u)), u=lots[1 ]):\nTitres:= [%, seq( cat(seq(Titres[j], j=u)), u=lots[2])]; ASSERT(n ops(%)=8); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 115 "mA:= copy(m ma);\nmA:= Matrix([mA, seq(procol(op(u)), u=lots[1])]);\nmA:= Matrix([ mA, seq(procol(op(u)), u=lots[2])]);" }}}{EXCHG {PARA 0 "" 0 "" {MPLTEXT 0 21 322 "Titres:= [Cte, \"B\",\"1\",\"2\",\"3\", \"B1\",\"B2 \",\"B3\",\"12\",\"13\",\"23\",\"B12\",\"B13\",\"B23\",\"123\",\"B123 \"];ASSERT(nops(%)=16);\nmA:= copy(mma);\nmA:= Matrix([mA, procol(2,3) ,procol(2,4),procol(2,5),procol(3,4),procol(3,5),procol(4,5)]);\nmA:= \+ Matrix([mA, procol(6,4),procol(6,5),procol(7,5),procol(9,5)]);\nmA:= M atrix([mA, procol(12,5)]);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 43 "mF:= 1/(Transpose(mA).mA).Transpose(mA).mb;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "convert(mF,list):\nli:= map( Z-> Float(round(100*Z),-2), %);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "jjx, iix:= Dimension(mA);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 208 "for i to iix do ; zoni:= i..i;\nc||(Titres[i]):= Vec tor([(0$(lhs(zoni)-1)), mF[zoni],(0$(iix-rhs(zoni)))]): \nr||(Titres[i ]):= evalf(mA.%);\nsort(convert(\{Entries(%)\}, list)); print( cat(` r `, Titres[i])=%);\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 317 " graphe:= NULL: x:= 1: # x est l'abscisse du point \nfor i from 2 to ii x do\nzonei:=i..i; tmp:= convert(mF[zonei], list); how:= nops(%);\n[[x , -add(k, k=tmp)], seq([x+k, tmp[k]], k=1..how)];\ngraphe:= graphe, pl ots[pointplot](%, style='line'), plots[pointplot](%);\nx:=x+how+2 \nod : graphe:= graphe, plot(0, y=0..x): x:='x':" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 205 "select(has, map(op, [graphe]), POINTS): \nticks:= ([min, max]@op@map2)(op, 2, map(op, %));\ndisplayg(graphe, tickmarks= [[seq(3*i-4.5=Titres[i], i=2..iix)], evalf(ticks,2)], axes=boxed); inf luences_inter:= %:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 188 "xres _ini(iix, [` #`, ` min`, ` max`, ` delta`, `sd`]):\nxres_run(proc(i) c onvert(r||(Titres[i]), list); evalf(xab(%)): \n[Titres[i], op(%), rhs( %)-lhs(%), (sqrt@var)(%%)] end): map(cut4,%);" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 0 "" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 " " }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 13 "Special light" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "ix; jx; Dimension(ma);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 134 "procol:= proc(u) local tmp, i; tmp := Column(mA,u);\nfor i from 2 to nargs do \ntmp:= zip((a,b)->a*b, t mp, Column(mA,args[i])) od; end:" }}}{EXCHG {PARA 0 "" 0 "" {MPLTEXT 0 21 139 "8=234, 7=245, 6=345;\n28=34, 38=24, 48=23, 27=45, 47=25, 57= 24, 36=45, 46=35, 56=34,\n 67=23, 68=25, \+ 78=35 " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 94 "mA:= copy(mm a): mA:= Matrix([mA, procol(3,4,5)]): \nTranspose(op(NullSpace(%))); \+ Dimension(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}}{MARK "22" 0 }{VIEWOPTS 1 1 0 3 2 1804 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }