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