I. INTRODUCTION

It is well known that the expression ${\mathrm{J}}_0$ of Bessel equations is the graphic most widely published in Differential Equation textbooks. Every time we talk about the vibration of a circular membrane the graphics of Figure 1is always referred to. Mathematical software packages such as MapleTM always includes this curve in their libraries. If you hit once a circular membrane, it is easy to assume that the first undulation will have a higher amplitude than the following undulations which will decrease gradually as shown also in Figure 1.

In this work we proceed to find the $\mathrm{Jv}(\mathbf{x})$ functions of the first kind and order v to satisfy the second order differential equation of Bessel.

To study acoustical radiation and vibration analysis one usually solves problems in cylindrical coordinates that are associated with Bessel functions of integer order. We can solve also spherical problems with half-integer order Bessel Function.

II. BESSEL FUNCTIONS

a) Differential Equation of Oder $\nu$

Wefocus this work on Besselfunctions of the first kind of positive order and real arguments[1]for

Assuming

which leads to

Differentiating and substituting

Writing out the differential equation give us

Substituting this back to the differential equation gives

Leading to

If ${\mathbf{r}}_1=\mathbf{v}$

giving

Writing this out

So that

Or

$c_1=0$ in (1) leads ${\mathrm{toc}}_3=c_5=c_7=\dots =0$ so, for $k=0,2,4$. after making $k+2=2n$, $n=1$, 2, 3.. we will have that

The coefficients with even index are determined by the following formula:

i. The Gamma Function

Since we can select ${\mathbf{a}}_0$, we take it to be [2]

The Gamma function $\mathrm{\Gamma }$ is defined as

Where

Then we can write

The solution for this $c_n$ is

b) Half Integer Order [3]

c) Spherical Bessel [4]

If we solve the Helmholtz's equation we obtain the following differential

One of the solutions of (3) is

The Bessel function of half-integral order is used to define the important function:

$J_n$ is called the spherical Bessel function of the first kind. For $n=0$ we can see that (1) becomes [5]

III. PLOTTING WITH MAPLE

a) Plotting $\sin \left(x\right)/x$

Figure 2: Plot of $\sin \left(\mathbf{x}\right)/\mathbf{x}$

b) Plotting of $J(-1/2)$

Figure 2: Plot of $\mathrm{Sin}(2/\pi x{)}^2\mathrm{Sin}(x)$

c) Plotting $J(1/2)$ with variations

Figure 4: Plot of Two variations of $\mathrm{J}(1/2)$

IV. COMPARING BESSEL WITH A LINEAR MODEL UNDERDAMPED

Given that

V. ANIMATION OF BESSEL ON MAPLE SOFTWARE

animate(wave, $[uC(0,3)],t=0..P(0,3))$;

1*1234567890

1

VI. CONCLUSION

A study of Bessel equations can lead us to new applications in acoustics. By changing their parameters, we can determine the ratio between a membrane dimension and the acoustic power it can radiate or the way the membrane will vibrate. For example, we can find analogies between the elasticity of a trampoline jumper and a loudspeaker membrane. We can also study suspension, material thickness, and surfaces and compare them to computer animations which allow us to explore these analogies. Thus, educational institutions with low budgets can make physics experiments on vibration and mechanics at more affordable prices using computer animations. Knowing Bessel equation and its relationship with Helmholtz's formula we can extend the application to the resonance phenomena too.

ACKNOWLEDGMENT

The authors thank Dr. James S. Sochacki for his invaluable assistance and his thorough review of this paper. Dr. Sochacki is a James Madison University Emeritus

Professor in the department of Mathematics and Statistics. Dr. Sochacki is also coauthor of the Parker-Sochacki's algorithm widely for solving systems of ordinary Differential Equations (https://en.wikipedia.org/wiki/Parker%E2%80%93Sochacki_method).