login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit

Submitted by : (unknown) at: 2007-11-17T22:23:33-08:00 (9 years ago)
Name :
Axiom Version :
Category : Severity : Status :
Optional subject :  
Optional comment :

fricas
a := operator 'a

\label{eq1}a(1)
fricas
seriesSolve((2*x-2)*D(a(x),x)+4*a(x),a,x=0,[1])
fricas
Compiling function %D with type List(UnivariateTaylorSeries(
      Expression(Integer),x,0)) -> UnivariateTaylorSeries(Expression(
      Integer),x,0)

\label{eq2}\begin{array}{@{}l}
\displaystyle
1 +{2 \  x}+{3 \ {{x}^{2}}}+{4 \ {{x}^{3}}}+{5 \ {{x}^{4}}}+{6 \ {{x}^{5}}}+{7 \ {{x}^{6}}}+{8 \ {{x}^{7}}}+{9 \ {{x}^{8}}}+ 
\
\
\displaystyle
{{10}\ {{x}^{9}}}+{{11}\ {{x}^{10}}}+{O \left({{x}^{11}}\right)}
(2)
Type: UnivariateTaylorSeries?(Expression(Integer),x,0)

does not work. A workaround within the interpreter is to provide a function like:

  fixedPointExquo(s, p) == fixedPointExquo(s, p::UTS(EXPR INT, x, 0))$UTSODE(EXPR INT,UTS(EXPR INT, x, 0))

but that won't work in spad.

fricas
a := operator 'a

\label{eq3}a(3)
fricas
fixedPointExquo(s, p) == fixedPointExquo(s, p::UTS(EXPR INT, x, 0))$UTSODE(EXPR INT,UTS(EXPR INT, x, 0))
Type: Void
fricas
p := seriesSolve((2*x-2)*D(a(x),x)+4*a(x),a,x=0,[1])
fricas
Compiling function %H with type List(UnivariateTaylorSeries(
      Expression(Integer),x,0)) -> UnivariateTaylorSeries(Expression(
      Integer),x,0)

\label{eq4}\begin{array}{@{}l}
\displaystyle
1 +{2 \  x}+{3 \ {{x}^{2}}}+{4 \ {{x}^{3}}}+{5 \ {{x}^{4}}}+{6 \ {{x}^{5}}}+{7 \ {{x}^{6}}}+{8 \ {{x}^{7}}}+{9 \ {{x}^{8}}}+ 
\
\
\displaystyle
{{10}\ {{x}^{9}}}+{{11}\ {{x}^{10}}}+{O \left({{x}^{11}}\right)}
(4)
Type: UnivariateTaylorSeries?(Expression(Integer),x,0)
fricas
(2*x-2)*D(p,x)+4*p

\label{eq5}\begin{array}{@{}l}
\displaystyle
{4 \  x}+{{\left({{12}\  x}- 4 \right)}\  x}+{{\left({{24}\  x}-{12}\right)}\ {{x}^{2}}}+{{\left({{40}\  x}-{24}\right)}\ {{x}^{3}}}+ 
\
\
\displaystyle
{{\left({{60}\  x}-{40}\right)}\ {{x}^{4}}}+{{\left({{84}\  x}-{60}\right)}\ {{x}^{5}}}+{{\left({{112}\  x}-{84}\right)}\ {{x}^{6}}}+ 
\
\
\displaystyle
{{\left({{144}\  x}-{112}\right)}\ {{x}^{7}}}+{{\left({{180}\  x}-{144}\right)}\ {{x}^{8}}}+{{\left({{220}\  x}-{180}\right)}\ {{x}^{9}}}+ 
\
\
\displaystyle
{{\left({{264}\  x}-{220}\right)}\ {{x}^{10}}}+{O \left({{x}^{11}}\right)}
(5)
Type: UnivariateTaylorSeries?(Expression(Integer),x,0)

Looks like the case of two different x's again.

I have a better fix, but I'm still not happy. And seriesSolve still crashes quite often...

I'll send it per mail,

Martin

UnivariateTaylorSeriesODESolver is part of the Hidden group in exposed.lsp.pamphlet. Exposing this package prevents the error:
fricas
)clear completely
All user variables and function definitions have been cleared. All )browse facility databases have been cleared. Internally cached functions and constructors have been cleared. )clear completely is finished.
fricas
)expose UnivariateTaylorSeriesODESolver
UnivariateTaylorSeriesODESolver is now explicitly exposed in frame initial

fricas
a := operator 'a

\label{eq6}a(6)
fricas
seriesSolve((2*x-2)*D(a(x),x)+4*a(x),a,x=0,[1])
fricas
Compiling function %D with type List(UnivariateTaylorSeries(
      Expression(Integer),x,0)) -> UnivariateTaylorSeries(Expression(
      Integer),x,0)

\label{eq7}\begin{array}{@{}l}
\displaystyle
1 +{2 \  x}+{3 \ {{x}^{2}}}+{4 \ {{x}^{3}}}+{5 \ {{x}^{4}}}+{6 \ {{x}^{5}}}+{7 \ {{x}^{6}}}+{8 \ {{x}^{7}}}+{9 \ {{x}^{8}}}+ 
\
\
\displaystyle
{{10}\ {{x}^{9}}}+{{11}\ {{x}^{10}}}+{O \left({{x}^{11}}\right)}
(7)
Type: UnivariateTaylorSeries?(Expression(Integer),x,0)

fixed in FriCAS revision 423 --kratt6, Tue, 13 Jan 2009 02:05:24 -0800 reply
Status: open => fixed somewhere

Patch From Fricas --japp, Wed, 04 Feb 2009 17:16:40 -0800 reply
http://fricas.svn.sourceforge.net/viewvc/fricas?view=rev&revision=423

Status: fixed somewhere => fix proposed




  Subject:   Be Bold !!
  ( 13 subscribers )  
Please rate this page: