\mnb150ÿ{\rtf1\ansi\deff0\deftab720{\fonttbl{\f0\fswiss MS Sans Serif;}{\f1\froman\fcharset2 Symbol;}{\f2\fswiss\fprq2 System;}{\f3\fmodern\fprq1 Courier New;}{\f4\froman\fcharset1 Times New Roman;}} {\colortbl\red0\green0\blue0;\red0\green0\blue255;\red255\green0\blue0;} \deflang1031\pard\ri4\plain\f4\fs22\cf0 01. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}a:=poly(x^4+x^2+1,[x]) \par \pard\ri4\plain\f4\fs22\cf0 02. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}punkte:=[-1,0,1] \par \pard\ri4\plain\f4\fs22\cf0 03. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}werte:=[a(punkte[i]) $ i=1..nops(punkte)] \par \pard\ri4\plain\f4\fs22\cf0 04. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}divisors:=proc(z) \par \pard\li600\ri1\fi-300\plain\f3\fs22\cf2 local divi; \par begin \par divi:=numlib::divisors(z); \par [op(divi),op(map(divi,_negate))] \par end_proc: \par \pard\ri4\plain\f4\fs22\cf0 05. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}divisorliste:=[divisors(werte[k]) $ k=1..nops(werte)] \par \pard\ri4\plain\f4\fs22\cf0 06. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}tabelle:=combinat::cartesianProduct(op(divisorliste)) \par \pard\ri4\plain\f4\fs22\cf0 07. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}poltabelle:=map(tabelle,y->interpolate([-1,0,1],y,x)) \par \pard\ri4\plain\f4\fs22\cf0 08. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}reste:=map(poltabelle,y->divide(a,y)[2]) \par \pard\ri4\plain\f4\fs22\cf0 09./10. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}fac:=[]: \par \pard\li600\ri1\fi-300\plain\f3\fs22\cf2 for i from 1 to nops(reste) do \par if reste[i]=poly(0,[x]) \par then fac:=[op(fac),poltabelle[i]] \par end_if: \par end_for: \par fac \par \pard\ri4\plain\f4\fs22\cf0 11. \par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f3\fs22\cf2 {\pntext\f1\'b7\tab}factor(a) \par }