• To: mathgroup at smc.vnet.net
• From: Chris Chiasson <chris.chiasson at gmail.com>
• Date: Tue, 19 Apr 2005 04:55:12 -0400 (EDT)
• References: <200504180708.DAA20024@smc.vnet.net>
• Reply-to: Chris Chiasson <chris.chiasson at gmail.com>
• Sender: owner-wri-mathgroup at wolfram.com

```Try this instead:
Integrate[(hillb*((f0 + hilla)/(2*
Pi*f*a*Cos[ArcSin[x/a]] + hillb))) - hilla, {x, -a, a}]
% /. {hilla -> 3, hillb -> 50, f0 -> 8, a -> 1, f -> 1}
% // N

On 4/18/05, Chris Chiasson <chris.chiasson at gmail.com> wrote:
> Of course, that might be expected when none of the fuctions have their
> arguments specified with brackets, as they are supposed to...
>
> On 4/18/05, Chris Chiasson <chris.chiasson at gmail.com> wrote:
> > In Mathematica 5.0, if I evaluate the else branch (not shown above) of
> > the output from the above mentioned integration command, the kernel
> > crashes.
> >
> > On 4/18/05, Jim Martin <jim.martin at utah.edu> wrote:
> > > Hello Mathematica Experts:
> > >
> > > I am a biomechanist and work mostly in the area of muscle contraction. I
> > > do a lot of numerical computations using excel, but right now I need an
> > > analytical solution that represents force as a function of position
> > > integrated over a shortening amplitude. I downloaded a trial version of
> > > Mathematica and have tried to obtain a solution for this:
> > >
> > > Integrate[(hillb*((f0 + hilla)/(2*pi*f*a*Cos(ArcSin(x/a)) + hillb))) -
> > > hilla, {x, -a, a}]
> > >
> > > Mathematica returns this:
> > > (-4 a ArcSin Cos f hilla pi + (f0 + hilla) hillb (-Log[hillb - 2 a
> > >      ArcSin Cos f pi] + Log[hillb + 2 a ArcSin Cos f pi]))/(4 a ArcSin
> > > Cos f pi)
> > >
> > > I know the line wrap makes this hard to read so please feel free to
> > > email me and I can send you the output as a picture.
> > >
> > > In a sample data set, hilla=3, hillb=50, f0=8, a=1, f=1
> > >
> > > I can numerically integrate this function and obtain a value for that
> > > sample data set of 14.04. When I put those sample values into the
> > > solution that Mathematica produces, I get 10.01.
> > >
> > > Can any of you please give a hand here? I must be making some simple
> > > Mathematica-beginner error but I just can't see it.
> > >
> > > In Mathematica, Log is Log to base e, right (LN in excel)? Did I use
> > > variables that have intrinsic functions in Mathematica? Maybe I am
> > > misunderstanding the output with regard to implicit parentheses etc.
> > > Any help appreciated!
> > >
> > > Thanks,
> > >
> > > Jim
> > >
> > >
> >
> > --
> > Chris Chiasson
> > Kettering University
> > Mechanical Engineering
> > 1 810 265 3161
> >
>
> --
> Chris Chiasson
> Kettering University
> Mechanical Engineering
> 1 810 265 3161
>

--
Chris Chiasson
Kettering University
Mechanical Engineering