Arnold, OpenAxiom is a fork of the Axiom project. The OpenAxiom mailing list is open-axiom-devel@lists.sourceforge.net If you are using OpenAxiom you might expect to get more specific answers there.
Regards, Bill Page. On Tue, Jun 16, 2009 at 4:20 AM, Arnold Doray wrote: > I'm trying to create a set of orthogonal polynomials using the > Gram-Schmidt process, but I hit a problem: > > dot(f,g) == integrate(f*g*x^2,x=-1..1) > proj(f,g) == dot(f,g)*f/dot(f,f) > > p0 := 1 > p1 := x - proj(p0,x) > p2 := x^2 - proj(p0,x^2) - proj(p1,x^2) > > p0 and p1 evaluate correctly, but p2 hits this error: > > ; (DEFUN |*2;dot;5;frame1394| ...) is being compiled. > ;; The variable |*2;dot;5;frame1394;MV| is undefined. > ;; The compiler will assume this variable is a global. > Internal Error > The function coerce with signature OrderedCompletion Integer -> > Expression OrderedCompletion Integer is missing from domain > Expression(OrderedCompletion (Integer)) > > I've also tried to define p(n) as a recursive function like so: > > p(0) == 1 > p(n | n > 0) == x^n - sum(proj(p(k),x^n),k=0..n-1) > > But get a stack overflow: (eg, evaluating p(1)) > > OpenAxiom will attempt to step through and interpret the code. > Cannot compile map: proj > We will attempt to interpret the code. > Cannot compile map: p > We will attempt to interpret the code. > > >> System error: > Bind stack overflow. > > I would appreciate any pointers where I am going wrong. I'm an Axiom > newbie. These calculations were done using open-axiom 1.2.1 > > Thanks > Arnold > > > > _______________________________________________ > Axiom-mail mailing list > axiom-m...@nongnu.org > http://lists.nongnu.org/mailman/listinfo/axiom-mail > ------------------------------------------------------------------------------ Crystal Reports - New Free Runtime and 30 Day Trial Check out the new simplified licensing option that enables unlimited royalty-free distribution of the report engine for externally facing server and web deployment. http://p.sf.net/sfu/businessobjects _______________________________________________ open-axiom-devel mailing list open-axiom-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/open-axiom-devel
