MathGroup Archive 1989

[Date Index] [Thread Index] [Author Index]

Search the Archive

Mathematica Problem

  • Subject: Mathematica Problem
  • From: jack at sun.acs.udel.edu
  • Date: Mon, 19 Jun 89 14:51:49 -0500
  • Apparently-to: mathgroup-out at yoda.ncsa.uiuc.edu

Perhaps someone can help solve the following Mathematica problem:

     Given y' = y, 

solve using a series solution and prove that the result is E[x].  

It is proving that the result is E[x] that I have not been able to do.
Here is a sample session:

In[1]:= Derivative[n_][y][x_]:=Derivative[n-1][y][x]

In[2]:= Series[y[x],{x,0,5}]

                              2         3         4         5
                        y[0] x    y[0] x    y[0] x    y[0] x        6
Out[2]= y[0] + y[0] x + ------- + ------- + ------- + ------- + O[x]
                           2         6        24        120

In[3]:= seriesbyy0[x_,n_]:=Normal[Expand[Cancel[Series[y[x],{x,0,n}]/y[0]]]]

In[4]:= expseries[x_,n_]:=Sum[x^i/i!,{i,0,n}]

In[5]:= seriesbyy0[x,6]

                 2    3    4    5     6
                x    x    x    x     x
Out[5]= 1 + x + -- + -- + -- + --- + ---
                2    6    24   120   720

In[6]:= expseries[x,6]

                 2    3    4    5     6
                x    x    x    x     x
Out[6]= 1 + x + -- + -- + -- + --- + ---
                2    6    24   120   720


In[7]:= seriesbyy0[x,100]-expseries[x,100]

Out[7]= 0


As shown, I can prove the results are identical for a constant value of
the argument "n", but not for general n.  I realize that the identity is 
proven by visual inspection, but I would like to prove it by showing the 
identity 

	    seriesbyy0[x,n]  ==  expseries[x,n]

evaluates to true.  I have tried LogicalExpand, but to no avail.


				Jack Seltzer
                                jack at sun.acs.udel.edu
 



  • Prev by Date: Mathematica Postscript Macros
  • Next by Date: Re: Mathematica Problem
  • Previous by thread: Mathematica Postscript Macros
  • Next by thread: Re: Mathematica Problem