Re: [Axiom-mail] Coersion error
Tim Daly <[email protected]> Tue, 16 Jun 2009 06:19:33 -0400
| Newsgroups | gmane.comp.mathematics.axiom.devel,gmane.comp.mathematics.axiom.user |
|---|---|
| Message-ID | <[email protected]> |
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 > [email protected] > http://lists.nongnu.org/mailman/listinfo/axiom-mail > > You are not doing anything wrong. The error message says that there is a missing function (coerce) from one domain to another. I will look at the problem and get back to you. Tim