Funciones bessel

Código fuente para generar la serie de potencias de J0 y J1

Program Js
Implicit none
Real, allocatable:: JJ0(:), JJ1(:), xx(:)
real x, b, h, J0, fact,a, factt, J1, aa
integer i,j ,m,n,k
m=21
a=0
b=7.7
h=0.1
n=int((b-a)/h)+2
fact=1
factt=1
x=0
allocate(JJ0(0:n), JJ1(0:n), xx(0:n))

do j=1,n
xx(0)=x
J0=((((-1)**0)*(x**(0)))/((2**(0))*(1**2)))J1=((((1)**0)*(x**(2*0+1)))/((2**(2*0+1))*fact*factt))
do i=1,m
aa=i
fact=fact*i
factt=factt*(i+1)

J0=J0+ ((((-1)**aa)*(x**(2*aa)))/((2**(2*aa))*(fact**2)))

J1=J1+ ((((-1)**aa)*(x**(2*aa+1)))/((2**(2*aa+1))*fact*factt))
enddo
xx(j)=x
fact=1
factt=1
Print*, x, J0, J1
x=h*j+a
enddo
endprogram Js

******X*******J0********J1**********
0.0000000 1.000000 0.0000000
0.1000000 0.99750164.993752E-02
0.2000000 0.9900250 9.950083E-02
0.3000000 0.9776263 0.1483188
0.4000000 0.9603982 0.1960266
0.5000000 0.9384698 0.2422685
0.6000000 0.9120048 0.2867010
0.7000000 0.8812009 0.3289957
0.8000000 0.8462873 0.3688420
0.9000000 0.8075238 0.4059496
1.0000000.7651977 0.4400506
1.100000 0.7196220 0.4709024
1.200000 0.6711327 0.4982891
1.300000 0.6200860 0.5220232
1.400000 0.5668551 0.5419477
1.500000 0.5118276 0.5579365
1.600000 0.4554021 0.5698959
1.700000 0.3979848 0.5777653
1.800000 0.3399864 0.58151691.900000 0.2818186 0.5811570
2.000000 0.2238908 0.5767248
2.100000 0.1666068 0.5682922
2.200000 0.1103623 0.5559630
2.300000 5.5539776E-02 0.5398725
2.400000 2.5076792E-03 0.5201852
2.500000 -4.8383772E-02 0.4970941
2.600000 -9.6805029E-02 0.4708183
2.700000-0.1424494 0.4416014
2.800000 -0.1850360 0.4097092
2.900000 -0.2243116 0.3754275
3.000000 -0.2600519 0.3390590
3.100000 -0.2920642 0.3009211
3.200000 -0.3201883 0.2613432
3.300000 -0.3442961 0.2206634
3.400000 -0.3642955 0.1792257
3.500000 -0.3801278 0.13737753.600000 -0.3917691 9.5465444E-02
3.700000 -0.3992302 5.3833924E-02
3.800000 -0.4025564 1.2821008E-02
3.900000 -0.4018263 -2.724397E-02
4.000000 -0.3971498 -6.604331E-02
4.100000 -0.3886696 -0.1032733
4.200000 -0.3765574 -0.1386470
4.300000 -0.3610112 -0.1718967
4.400000-0.3422567 -0.2027754
4.500000 -0.3205426 -0.2310604
4.600000 -0.2961379 -0.2565522
4.700000 -0.2693305 -0.2790812
4.800000 -0.2404253 -0.2984994
4.900000 -0.2097381 -0.3146951
5.000000 -0.1775967 -0.3275793
5.100000 -0.1443349 -0.3370963
5.200000 -0.1102904 -0.34322295.300000 -7.5803004E-02 -0.3459610
5.400000 -4.1210510E-02 -0.3453438
5.500000 -6.8437452E-03 -0.3414388
5.600000 2.6971286E-02 -0.3343331
5.700000 5.9921224E-02 -0.3241474
5.800000 9.1703020E-02 -0.3110269
5.900000 0.1220335 -0.2951427
6.000000 0.1506453 -0.2766838
6.1000000.1772905 -0.2558658
6.200000 0.2017479 -0.2329163
6.300000 0.2238129 -0.2080870
6.400000 0.2433075 -0.1816381
6.500000 0.2600951 -0.1538405
6.600000 0.2740444 -0.1249794
6.700000 0.2850636 -9.5340960E-02
6.800000 0.2930970 -6.5217510E-02
6.900000 0.2981021 -3.4900077E-02...