forked from lgcrego/Dynemol
-
Notifications
You must be signed in to change notification settings - Fork 0
/
ckplm-Coul.F
9163 lines (9140 loc) · 539 KB
/
ckplm-Coul.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
* coeficientes ckplm para obtener los factores radiales
* de las distribuciones bicentricas a partir de los
* factores radiales de las funciones trasladadas
* Subrutinas generadas a partir del notebook integrales3C-5.nb
* Preparadas para valores de l hasta 20
*
* ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO!
*
* Los ckplm definidos aqui difieren de los que aparecen en el
* articulo de J. Comp. Chem. (2005) 26, 846-855
* (Translation of STO Charge Distributions) en in factor (-1)**m
*
* ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO! ¡OJO!
*
*
********************************************************
*
subroutine subckplma(kmax,kdim,ckplm)
implicit real * 8 (a-h,o-z)
* Coeficientes con 0 <= l <= 14
dimension ckplm(3311,0:kdim), den(0:20)
do k = 0, kmax
rk = k
r2k = rk+rk
* ckplm para l = 0
den(0) = 1.d0/((r2k+1.D0))
ckplm(1,k) = 1.D0 * den(0)
* ckplm para l = 1
den(0) = 1.d0/((r2k+1.D0)*(r2k+3.D0))
den(1) = 1.d0/((r2k-1.D0)*(r2k+1.D0))
ckplm(2,k) = 3.D0*(rk+1.D0) * den(0)
ckplm(3,k) = 3.D0*rk * den(1)
ckplm(4,k) = -3.D0 * den(0)
ckplm(5,k) = 3.D0 * den(1)
* ckplm para l = 2
den(0) = 1.d0/((r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0))
den(1) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k+3.D0))
den(2) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0))
ckplm(6,k) = 7.5D0*(rk+1.D0)*(rk+2.D0) * den(0)
ckplm(7,k) = 5.D0*(rk+1.D0)*rk * den(1)
ckplm(8,k) = 7.5D0*(rk-1.D0)*rk * den(2)
ckplm(9,k) = -5.D0*(rk+2.D0) * den(0)
ckplm(10,k) = 5.D0 * den(1)
ckplm(11,k) = 5.D0*(rk-1.D0) * den(2)
ckplm(12,k) = 1.25D0 * den(0)
ckplm(13,k) = -2.5D0 * den(1)
ckplm(14,k) = 1.25D0 * den(2)
* ckplm para l = 3
den(0) = 1.d0/((r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0))
den(1) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0))
den(2) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0))
den(3) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k-5.D0))
ckplm(15,k) = 17.5D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0) * den(0)
ckplm(16,k) = 10.5D0*(rk+1.D0)*(rk+2.D0)*rk * den(1)
ckplm(17,k) = 10.5D0*(rk-1.D0)*(rk+1.D0)*rk * den(2)
ckplm(18,k) = 17.5D0*(rk-1.D0)*(rk-2.D0)*rk * den(3)
ckplm(19,k) = -8.75D0*(rk+2.D0)*(rk+3.D0) * den(0)
ckplm(20,k) = -1.75D0*(rk+2.D0)*(rk-5.D0) * den(1)
ckplm(21,k) = 1.75D0*(rk-1.D0)*(rk+6.D0) * den(2)
ckplm(22,k) = 8.75D0*(rk-1.D0)*(rk-2.D0) * den(3)
ckplm(23,k) = 1.75D0*(rk+3.D0) * den(0)
ckplm(24,k) = -1.75D0*(rk+4.D0) * den(1)
ckplm(25,k) = -1.75D0*(rk-3.D0) * den(2)
ckplm(26,k) = 1.75D0*(rk-2.D0) * den(3)
ckplm(27,k) = -0.291666666666667D0 * den(0)
ckplm(28,k) = 0.875D0 * den(1)
ckplm(29,k) = -0.875D0 * den(2)
ckplm(30,k) = 0.291666666666667D0 * den(3)
* ckplm para l = 4
den(0) = 1.d0/((r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.
&D0))
den(1) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.
&D0))
den(2) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k+5.
&D0))
den(3) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0))
den(4) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k-5.D0)*(r2k-7.
&D0))
ckplm(31,k) = 39.375D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0) * d
&en(0)
ckplm(32,k) = 22.5D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*rk * den(1)
ckplm(33,k) = 20.25D0*(rk-1.D0)*(rk+1.D0)*(rk+2.D0)*rk * den(2)
ckplm(34,k) = 22.5D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*rk * den(3)
ckplm(35,k) = 39.375D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*rk * den(4)
ckplm(36,k) = -15.75D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0) * den(0)
ckplm(37,k) = -2.25D0*(rk+2.D0)*(rk+3.D0)*(-7.D0+2.D0*rk) * den(1)
ckplm(38,k) = 20.25D0*(rk-1.D0)*(rk+2.D0) * den(2)
ckplm(39,k) = 2.25D0*(rk-1.D0)*(rk-2.D0)*(9.D0+2.D0*rk) * den(3)
ckplm(40,k) = 15.75D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0) * den(4)
ckplm(41,k) = 2.625D0*(rk+3.D0)*(rk+4.D0) * den(0)
ckplm(42,k) = -1.5D0*(rk+3.D0)*(rk+7.D0) * den(1)
ckplm(43,k) = -2.25D0*(rk*rk+rk-9.D0) * den(2)
ckplm(44,k) = -1.5D0*(rk-2.D0)*(rk-6.D0) * den(3)
ckplm(45,k) = 2.625D0*(rk-2.D0)*(rk-3.D0) * den(4)
ckplm(46,k) = -0.375D0*(rk+4.D0) * den(0)
ckplm(47,k) = 0.375D0*(9.D0+2.D0*rk) * den(1)
ckplm(48,k) = -3.375D0 * den(2)
ckplm(49,k) = -0.375D0*(-7.D0+2.D0*rk) * den(3)
ckplm(50,k) = 0.375D0*(rk-3.D0) * den(4)
ckplm(51,k) = 0.046875D0 * den(0)
ckplm(52,k) = -0.1875D0 * den(1)
ckplm(53,k) = 0.28125D0 * den(2)
ckplm(54,k) = -0.1875D0 * den(3)
ckplm(55,k) = 0.046875D0 * den(4)
* ckplm para l = 5
den(0) = 1.d0/((r2k+11.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7
&.D0)*(r2k+9.D0))
den(1) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.
&D0)*(r2k+9.D0))
den(2) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k+5.
&D0)*(r2k+7.D0))
den(3) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0))
den(4) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k-7.D0))
den(5) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k-5.D0)*(r2k-7.
&D0)*(r2k-9.D0))
ckplm(56,k) = 86.625D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk
&+5.D0) * den(0)
ckplm(57,k) = 48.125D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*rk
&* den(1)
ckplm(58,k) = 41.25D0*(rk-1.D0)*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*rk *
& den(2)
ckplm(59,k) = 41.25D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*rk *
& den(3)
ckplm(60,k) = 48.125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk-3.D0)*rk
&* den(4)
ckplm(61,k) = 86.625D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*rk
&* den(5)
ckplm(62,k) = -28.875D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0) *
&den(0)
ckplm(63,k) = -9.625D0*(rk+2.D0)*(rk-3.D0)*(rk+3.D0)*(rk+4.D0) * d
&en(1)
ckplm(64,k) = -2.75D0*(rk-14.D0)*(rk-1.D0)*(rk+2.D0)*(rk+3.D0) * d
&en(2)
ckplm(65,k) = 2.75D0*(rk+15.D0)*(rk-1.D0)*(rk-2.D0)*(rk+2.D0) * de
&n(3)
ckplm(66,k) = 9.625D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk+4.D0) * de
&n(4)
ckplm(67,k) = 28.875D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0) * d
&en(5)
ckplm(68,k) = 4.125D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0) * den(0)
ckplm(69,k) = -1.375D0*(rk+12.D0)*(rk+3.D0)*(rk+4.D0) * den(1)
ckplm(70,k) = -2.75D0*(rk+3.D0)*(-13.D0+(rk+3.D0)*rk) * den(2)
ckplm(71,k) = -2.75D0*(rk-2.D0)*(-15.D0+(rk-1.D0)*rk) * den(3)
ckplm(72,k) = -1.375D0*(rk-11.D0)*(rk-2.D0)*(rk-3.D0) * den(4)
ckplm(73,k) = 4.125D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0) * den(5)
ckplm(74,k) = -0.515625D0*(rk+4.D0)*(rk+5.D0) * den(0)
ckplm(75,k) = 0.0572916666666667D0*(rk+4.D0)*(81.D0+13.D0*rk) * de
&n(1)
ckplm(76,k) = 0.34375D0*(-48.D0+(rk-7.D0)*rk) * den(2)
ckplm(77,k) = -0.34375D0*(-40.D0+(rk+9.D0)*rk) * den(3)
ckplm(78,k) = -0.0572916666666667D0*(rk-3.D0)*(-68.D0+13.D0*rk) *
&den(4)
ckplm(79,k) = 0.515625D0*(rk-3.D0)*(rk-4.D0) * den(5)
ckplm(80,k) = 0.0572916666666667D0*(rk+5.D0) * den(0)
ckplm(81,k) = -0.0572916666666667D0*(16.D0+3.D0*rk) * den(1)
ckplm(82,k) = 0.114583333333333D0*(rk+11.D0) * den(2)
ckplm(83,k) = 0.114583333333333D0*(rk-10.D0) * den(3)
ckplm(84,k) = -0.0572916666666667D0*(-13.D0+3.D0*rk) * den(4)
ckplm(85,k) = 0.0572916666666667D0*(rk-4.D0) * den(5)
ckplm(86,k) = -0.00572916666666667D0 * den(0)
ckplm(87,k) = 0.0286458333333333D0 * den(1)
ckplm(88,k) = -0.0572916666666667D0 * den(2)
ckplm(89,k) = 0.0572916666666667D0 * den(3)
ckplm(90,k) = -0.0286458333333333D0 * den(4)
ckplm(91,k) = 0.00572916666666667D0 * den(5)
* ckplm para l = 6
den(0) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+
&5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(1) = 1.d0/((r2k+11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5
&.D0)*(r2k+7.D0)*(r2k+9.D0))
den(2) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k+5.
&D0)*(r2k+7.D0)*(r2k+9.D0))
den(3) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k+7.D0))
den(4) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k-7.D0))
den(5) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k-7.D0)*(r2k-9.D0))
den(6) = 1.d0/((r2k-11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k-5
&.D0)*(r2k-7.D0)*(r2k-9.D0))
ckplm(92,k) = 187.6875D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(
&rk+5.D0)*(rk+6.D0) * den(0)
ckplm(93,k) = 102.375D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(r
&k+5.D0)*rk * den(1)
ckplm(94,k) = 85.3125D0*(rk-1.D0)*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(r
&k+4.D0)*rk * den(2)
ckplm(95,k) = 81.25D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*(rk+
&3.D0)*rk * den(3)
ckplm(96,k) = 85.3125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*(r
&k-3.D0)*rk * den(4)
ckplm(97,k) = 102.375D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk-3.D0)*(r
&k-4.D0)*rk * den(5)
ckplm(98,k) = 187.6875D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(
&rk-5.D0)*rk * den(6)
ckplm(99,k) = -53.625D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(r
&k+6.D0) * den(0)
ckplm(100,k) = -4.875D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(-
&11.D0+4.D0*rk) * den(1)
ckplm(101,k) = -8.125D0*(rk-1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(r
&k-9.D0) * den(2)
ckplm(102,k) = 81.25D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk+3.D0) * d
&en(3)
ckplm(103,k) = 8.125D0*(rk+10.D0)*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(r
&k-3.D0) * den(4)
ckplm(104,k) = 4.875D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(15
&.D0+4.D0*rk) * den(5)
ckplm(105,k) = 53.625D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(r
&k-5.D0) * den(6)
ckplm(106,k) = 6.703125D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)
&* den(0)
ckplm(107,k) = -1.21875D0*(rk+22.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)
& * den(1)
ckplm(108,k) = -0.203125D0*(rk+3.D0)*(rk+4.D0)*(-306.D0+rk*(91.D0+
&17.D0*rk)) * den(2)
ckplm(109,k) = -4.0625D0*(rk-2.D0)*(rk+3.D0)*(rk-4.D0)*(rk+5.D0) *
& den(3)
ckplm(110,k) = -0.203125D0*(rk-2.D0)*(rk-3.D0)*(-380.D0+rk*(-57.D0
&+17.D0*rk)) * den(4)
ckplm(111,k) = -1.21875D0*(rk-21.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)
& * den(5)
ckplm(112,k) = 6.703125D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)
&* den(6)
ckplm(113,k) = -0.744791666666667D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)
&* den(0)
ckplm(114,k) = 0.203125D0*(rk+4.D0)*(rk+5.D0)*(33.D0+4.D0*rk) * de
&n(1)
ckplm(115,k) = 0.609375D0*(rk+4.D0)*(-43.D0+(rk-2.D0)*rk) * den(2)
ckplm(116,k) = -2.03125D0*(-40.D0+3.D0*(rk+1.D0)*rk) * den(3)
ckplm(117,k) = -0.609375D0*(rk-3.D0)*(-40.D0+(rk+4.D0)*rk) * den(4
&)
ckplm(118,k) = -0.203125D0*(rk-3.D0)*(rk-4.D0)*(-29.D0+4.D0*rk) *
&den(5)
ckplm(119,k) = 0.744791666666667D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0) *
& den(6)
ckplm(120,k) = 0.0744791666666667D0*(rk+5.D0)*(rk+6.D0) * den(0)
ckplm(121,k) = -0.0135416666666667D0*(rk+5.D0)*(88.D0+13.D0*rk) *
&den(1)
ckplm(122,k) = 0.0338541666666667D0*(216.D0+(rk+47.D0)*rk) * den(2
&)
ckplm(123,k) = 0.135416666666667D0*(rk*rk+rk-50.D0) * den(3)
ckplm(124,k) = 0.0338541666666667D0*(170.D0+(rk-45.D0)*rk) * den(4
&)
ckplm(125,k) = -0.0135416666666667D0*(rk-4.D0)*(-75.D0+13.D0*rk) *
& den(5)
ckplm(126,k) = 0.0744791666666667D0*(rk-4.D0)*(rk-5.D0) * den(6)
ckplm(127,k) = -0.00677083333333334D0*(rk+6.D0) * den(0)
ckplm(128,k) = 0.00677083333333334D0*(25.D0+4.D0*rk) * den(1)
ckplm(129,k) = -0.0338541666666667D0*(rk+9.D0) * den(2)
ckplm(130,k) = 0.338541666666667D0 * den(3)
ckplm(131,k) = 0.0338541666666667D0*(rk-8.D0) * den(4)
ckplm(132,k) = -0.00677083333333334D0*(-21.D0+4.D0*rk) * den(5)
ckplm(133,k) = 0.00677083333333334D0*(rk-5.D0) * den(6)
ckplm(134,k) = 0.00056423611111111D0 * den(0)
ckplm(135,k) = -0.00338541666666667D0 * den(1)
ckplm(136,k) = 0.00846354166666666D0 * den(2)
ckplm(137,k) = -0.0112847222222222D0 * den(3)
ckplm(138,k) = 0.00846354166666666D0 * den(4)
ckplm(139,k) = -0.00338541666666667D0 * den(5)
ckplm(140,k) = 0.00056423611111111D0 * den(6)
* ckplm para l = 7
den(0) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k+1.D0)*(r2k
&+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(1) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k+
&3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(2) = 1.d0/((r2k+11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(3) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(4) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0))
den(5) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k-9.D0))
den(6) = 1.d0/((r2k-11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k-5.D0)*(r2k-7.D0)*(r2k-9.D0))
den(7) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-
&3.D0)*(r2k-5.D0)*(r2k-7.D0)*(r2k-9.D0))
ckplm(141,k) = 402.1875D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*
&(rk+5.D0)*(rk+6.D0)*(rk+7.D0) * den(0)
ckplm(142,k) = 216.5625D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*
&(rk+5.D0)*(rk+6.D0)*rk * den(1)
ckplm(143,k) = 177.1875D0*(rk-1.D0)*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*
&(rk+4.D0)*(rk+5.D0)*rk * den(2)
ckplm(144,k) = 164.0625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*
&(rk+3.D0)*(rk+4.D0)*rk * den(3)
ckplm(145,k) = 164.0625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*
&(rk-3.D0)*(rk+3.D0)*rk * den(4)
ckplm(146,k) = 177.1875D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*
&(rk-3.D0)*(rk-4.D0)*rk * den(5)
ckplm(147,k) = 216.5625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk-3.D0)*
&(rk-4.D0)*(rk-5.D0)*rk * den(6)
ckplm(148,k) = 402.1875D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*
&(rk-5.D0)*(rk-6.D0)*rk * den(7)
ckplm(149,k) = -100.546875D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D
&0)*(rk+6.D0)*(rk+7.D0) * den(0)
ckplm(150,k) = -7.734375D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)
&*(rk+6.D0)*(-13.D0+5.D0*rk) * den(1)
ckplm(151,k) = -6.328125D0*(rk-1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)
&*(rk+5.D0)*(-22.D0+3.D0*rk) * den(2)
ckplm(152,k) = -5.859375D0*(rk-1.D0)*(rk-27.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk+3.D0)*(rk+4.D0) * den(3)
ckplm(153,k) = 5.859375D0*(rk-1.D0)*(rk+28.D0)*(rk-2.D0)*(rk+2.D0)
&*(rk-3.D0)*(rk+3.D0) * den(4)
ckplm(154,k) = 6.328125D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0)*
&(rk-4.D0)*(25.D0+3.D0*rk) * den(5)
ckplm(155,k) = 7.734375D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*
&(rk-5.D0)*(18.D0+5.D0*rk) * den(6)
ckplm(156,k) = 100.546875D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0
&)*(rk-5.D0)*(rk-6.D0) * den(7)
ckplm(157,k) = 11.171875D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)
&*(rk+7.D0) * den(0)
ckplm(158,k) = -0.859375D0*(rk+3.D0)*(rk+4.D0)*(rk+52.D0)*(rk+5.D0
&)*(rk+6.D0) * den(1)
ckplm(159,k) = -0.234375D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(-462.D0+
&rk*(157.D0+19.D0*rk)) * den(2)
ckplm(160,k) = -5.859375D0*(rk-2.D0)*(rk+3.D0)*(rk+4.D0)*(-26.D0+(
&rk+3.D0)*rk) * den(3)
ckplm(161,k) = -5.859375D0*(rk-2.D0)*(rk-3.D0)*(rk+3.D0)*(-28.D0+(
&rk-1.D0)*rk) * den(4)
ckplm(162,k) = -0.234375D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(-600.D0+
&rk*(-119.D0+19.D0*rk)) * den(5)
ckplm(163,k) = -0.859375D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-51.D0
&)*(rk-5.D0) * den(6)
ckplm(164,k) = 11.171875D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)
&*(rk-6.D0) * den(7)
ckplm(165,k) = -1.1171875D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0
&) * den(0)
ckplm(166,k) = 0.0859375D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(117.D0+1
&1.D0*rk) * den(1)
ckplm(167,k) = 0.0703125D0*(rk+4.D0)*(rk+5.D0)*(-594.D0+rk*(9.D0+1
&3.D0*rk)) * den(2)
ckplm(168,k) = 0.1171875D0*(rk+4.D0)*(1250.D0+3.D0*rk*(-73.D0+(rk-
&22.D0)*rk)) * den(3)
ckplm(169,k) = -0.1171875D0*(rk-3.D0)*(-1400.D0+3.D0*(rk-1.D0)*(rk
&+26.D0)*rk) * den(4)
ckplm(170,k) = -0.0703125D0*(rk-3.D0)*(rk-4.D0)*(-590.D0+rk*(17.D0
&+13.D0*rk)) * den(5)
ckplm(171,k) = -0.0859375D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(-106.D0
&+11.D0*rk) * den(6)
ckplm(172,k) = 1.1171875D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)
& * den(7)
ckplm(173,k) = 0.1015625D0*(rk+5.D0)*(rk+6.D0)*(rk+7.D0) * den(0)
ckplm(174,k) = -0.0078125D0*(rk+5.D0)*(rk+6.D0)*(208.D0+25.D0*rk)
&* den(1)
ckplm(175,k) = -0.0234375D0*(rk+5.D0)*(-456.D0+(rk-65.D0)*rk) * de
&n(2)
ckplm(176,k) = 0.1171875D0*(-400.D0+rk*(-38.D0+(rk+13.D0)*rk)) * d
&en(3)
ckplm(177,k) = 0.1171875D0*(350.D0+rk*(-61.D0+(rk-10.D0)*rk)) * de
&n(4)
ckplm(178,k) = -0.0234375D0*(rk-4.D0)*(-390.D0+(rk+67.D0)*rk) * de
&n(5)
ckplm(179,k) = 0.0078125D0*(rk-4.D0)*(rk-5.D0)*(183.D0-25.D0*rk) *
& den(6)
ckplm(180,k) = 0.1015625D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0) * den(7)
ckplm(181,k) = -0.00846354166666666D0*(rk+6.D0)*(rk+7.D0) * den(0)
ckplm(182,k) = 0.00065104166666667D0*(rk+6.D0)*(325.D0+43.D0*rk) *
& den(1)
ckplm(183,k) = 0.001953125D0*(-1050.D0-rk*(239.D0+11.D0*rk)) * den
&(2)
ckplm(184,k) = -0.0162760416666667D0*(-138.D0+(rk-11.D0)*rk) * den
&(3)
ckplm(185,k) = 0.0162760416666667D0*(-126.D0+(rk+13.D0)*rk) * den(
&4)
ckplm(186,k) = 0.001953125D0*(822.D0+rk*(-217.D0+11.D0*rk)) * den(
&5)
ckplm(187,k) = 0.00065104166666667D0*(rk-5.D0)*(282.D0-43.D0*rk) *
& den(6)
ckplm(188,k) = 0.00846354166666666D0*(rk-5.D0)*(rk-6.D0) * den(7)
ckplm(189,k) = 0.00065104166666667D0*(rk+7.D0) * den(0)
ckplm(190,k) = 0.00065104166666667D0*(-36.D0-5.D0*rk) * den(1)
ckplm(191,k) = 0.005859375D0*(rk+9.D0) * den(2)
ckplm(192,k) = -0.00325520833333333D0*(rk+22.D0) * den(3)
ckplm(193,k) = -0.00325520833333333D0*(rk-21.D0) * den(4)
ckplm(194,k) = 0.005859375D0*(rk-8.D0) * den(5)
ckplm(195,k) = 0.00065104166666667D0*(31.D0-5.D0*rk) * den(6)
ckplm(196,k) = 0.00065104166666667D0*(rk-6.D0) * den(7)
ckplm(197,k) = -0.00004650297619048D0 * den(0)
ckplm(198,k) = 0.00032552083333333D0 * den(1)
ckplm(199,k) = -0.0009765625D0 * den(2)
ckplm(200,k) = 0.00162760416666667D0 * den(3)
ckplm(201,k) = -0.00162760416666667D0 * den(4)
ckplm(202,k) = 0.0009765625D0 * den(5)
ckplm(203,k) = -0.00032552083333333D0 * den(6)
ckplm(204,k) = 0.00004650297619048D0 * den(7)
* ckplm para l = 8
den(0) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k+17.D0)*(r2
&k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(1) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k-1.D0)*(r2k
&+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(2) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-
&3.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(3) = 1.d0/((r2k+11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(4) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k+9.D0))
den(5) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k-9.D0))
den(6) = 1.d0/((r2k-11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k-9.D0))
den(7) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-
&3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k-7.D0)*(r2k-9.D0))
den(8) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-15.D0)*(r2k-1.D0)*(r2k
&+1.D0)*(r2k-3.D0)*(r2k-5.D0)*(r2k-7.D0)*(r2k-9.D0))
ckplm(205,k) = 854.6484375D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D
&0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0) * den(0)
ckplm(206,k) = 455.8125D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*
&(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*rk * den(1)
ckplm(207,k) = 368.15625D0*(rk-1.D0)*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)
&*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*rk * den(2)
ckplm(208,k) = 334.6875D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*
&(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*rk * den(3)
ckplm(209,k) = 325.390625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk-3.D0)*(rk+3.D0)*(rk+4.D0)*rk * den(4)
ckplm(210,k) = 334.6875D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)*
&(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*rk * den(5)
ckplm(211,k) = 368.15625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)
&*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*rk * den(6)
ckplm(212,k) = 455.8125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk-3.D0)*
&(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*rk * den(7)
ckplm(213,k) = 854.6484375D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D
&0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*rk * den(8)
ckplm(214,k) = -189.921875D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D
&0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0) * den(0)
ckplm(215,k) = -37.984375D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0
&)*(rk+6.D0)*(rk+7.D0)*(-5.D0+2.D0*rk) * den(1)
ckplm(216,k) = -20.453125D0*(rk-1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0
&)*(rk+5.D0)*(rk+6.D0)*(-13.D0+2.D0*rk) * den(2)
ckplm(217,k) = -9.296875D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk+3.D0)
&*(rk+4.D0)*(rk+5.D0)*(-33.D0+2.D0*rk) * den(3)
ckplm(218,k) = 325.390625D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0
&)*(rk+3.D0)*(rk+4.D0) * den(4)
ckplm(219,k) = 9.296875D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0)*
&(rk+3.D0)*(rk-4.D0)*(35.D0+2.D0*rk) * den(5)
ckplm(220,k) = 20.453125D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0)
&*(rk-4.D0)*(rk-5.D0)*(15.D0+2.D0*rk) * den(6)
ckplm(221,k) = 37.984375D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)
&*(rk-5.D0)*(rk-6.D0)*(7.D0+2.D0*rk) * den(7)
ckplm(222,k) = 189.921875D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0
&)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0) * den(8)
ckplm(223,k) = 18.9921875D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0
&)*(rk+7.D0)*(rk+8.D0) * den(0)
ckplm(224,k) = -75.96875D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)
&*(rk+7.D0) * den(1)
ckplm(225,k) = -2.921875D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)
&*(-65.D0+2.D0*(rk+12.D0)*rk) * den(2)
ckplm(226,k) = -0.53125D0*(rk-2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*
&(-528.D0+rk*(83.D0+16.D0*rk)) * den(3)
ckplm(227,k) = -9.296875D0*(rk*rk+rk-35.D0)*(rk-2.D0)*(rk-3.D0)*(r
&k+3.D0)*(rk+4.D0) * den(4)
ckplm(228,k) = -0.53125D0*(rk-2.D0)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*
&(-595.D0+rk*(-51.D0+16.D0*rk)) * den(5)
ckplm(229,k) = -2.921875D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)
&*(-87.D0+2.D0*(rk-10.D0)*rk) * den(6)
ckplm(230,k) = 75.96875D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*
&(rk-6.D0) * den(7)
ckplm(231,k) = 18.9921875D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0
&)*(rk-6.D0)*(rk-7.D0) * den(8)
ckplm(232,k) = -1.7265625D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0
&)*(rk+8.D0) * den(0)
ckplm(233,k) = 0.575520833333333D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(
&rk+7.D0)*(27.D0+2.D0*rk) * den(1)
ckplm(234,k) = 0.1328125D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(-507.D0+
&rk*(29.D0+10.D0*rk)) * den(2)
ckplm(235,k) = 0.3984375D0*(rk+4.D0)*(rk+5.D0)*(642.D0+rk*(-153.D0
&+rk*(-23.D0+2.D0*rk))) * den(3)
ckplm(236,k) = -4.6484375D0*(rk-3.D0)*(rk+4.D0)*(-70.D0+3.D0*(rk+1
&.D0)*rk) * den(4)
ckplm(237,k) = -0.3984375D0*(rk-3.D0)*(rk-4.D0)*(-770.D0+rk*(-101.
&D0+rk*(29.D0+2.D0*rk))) * den(5)
ckplm(238,k) = -0.1328125D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(-526.D0
&+rk*(-9.D0+10.D0*rk)) * den(6)
ckplm(239,k) = -0.575520833333333D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*
&(rk-6.D0)*(-25.D0+2.D0*rk) * den(7)
ckplm(240,k) = 1.7265625D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)
&*(rk-7.D0) * den(8)
ckplm(241,k) = 0.143880208333333D0*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(
&rk+8.D0) * den(0)
ckplm(242,k) = -0.230208333333333D0*(rk+10.D0)*(rk+5.D0)*(rk+6.D0)
&*(rk+7.D0) * den(1)
ckplm(243,k) = -0.0265625D0*(rk-26.D0)*(rk+5.D0)*(rk+6.D0)*(23.D0+
&3.D0*rk) * den(2)
ckplm(244,k) = 0.0885416666666667D0*(rk+5.D0)*(-876.D0+rk*(-10.D0+
&(rk+27.D0)*rk)) * den(3)
ckplm(245,k) = 0.154947916666667D0*(2100.D0+(rk*rk+rk-122.D0)*(rk+
&1.D0)*rk) * den(4)
ckplm(246,k) = 0.0885416666666667D0*(rk-4.D0)*(840.D0+rk*(-61.D0+(
&rk-24.D0)*rk)) * den(5)
ckplm(247,k) = -0.0265625D0*(rk+27.D0)*(rk-4.D0)*(rk-5.D0)*(-20.D0
&+3.D0*rk) * den(6)
ckplm(248,k) = -0.230208333333333D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*
&(rk-9.D0) * den(7)
ckplm(249,k) = 0.143880208333333D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(
&rk-7.D0) * den(8)
ckplm(250,k) = -0.0110677083333333D0*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)
& * den(0)
ckplm(251,k) = 0.00221354166666667D0*(rk+6.D0)*(rk+7.D0)*(125.D0+1
&4.D0*rk) * den(1)
ckplm(252,k) = -0.006640625D0*(rk+6.D0)*(425.D0+rk*(73.D0+2.D0*rk)
&) * den(2)
ckplm(253,k) = -0.0110677083333333D0*(-1590.D0+rk*(-281.D0+rk*(9.D
&0+2.D0*rk))) * den(3)
ckplm(254,k) = 0.387369791666667D0*(rk-6.D0)*(rk+7.D0) * den(4)
ckplm(255,k) = 0.0110677083333333D0*(1302.D0+rk*(-293.D0+rk*(-3.D0
&+2.D0*rk))) * den(5)
ckplm(256,k) = 0.006640625D0*(rk-5.D0)*(354.D0+rk*(-69.D0+2.D0*rk)
&) * den(6)
ckplm(257,k) = -0.00221354166666667D0*(rk-5.D0)*(rk-6.D0)*(-111.D0
&+14.D0*rk) * den(7)
ckplm(258,k) = 0.0110677083333333D0*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)
&* den(8)
ckplm(259,k) = 0.0007905505952381D0*(rk+7.D0)*(rk+8.D0) * den(0)
ckplm(260,k) = -0.00021081349206349D0*(rk+7.D0)*(135.D0+16.D0*rk)
&* den(1)
ckplm(261,k) = 0.00221354166666667D0*(183.D0+2.D0*(rk+20.D0)*rk) *
& den(2)
ckplm(262,k) = -0.06640625D0*(rk+8.D0) * den(3)
ckplm(263,k) = -0.00368923611111111D0*(rk*rk+rk-147.D0) * den(4)
ckplm(264,k) = 0.06640625D0*(rk-7.D0) * den(5)
ckplm(265,k) = 0.00221354166666667D0*(145.D0+2.D0*(rk-18.D0)*rk) *
& den(6)
ckplm(266,k) = -0.00021081349206349D0*(rk-6.D0)*(-119.D0+16.D0*rk)
& * den(7)
ckplm(267,k) = 0.0007905505952381D0*(rk-6.D0)*(rk-7.D0) * den(8)
ckplm(268,k) = -0.00005270337301587D0*(rk+8.D0) * den(0)
ckplm(269,k) = 0.00005270337301587D0*(49.D0+6.D0*rk) * den(1)
ckplm(270,k) = -0.00036892361111111D0*(19.D0+2.D0*rk) * den(2)
ckplm(271,k) = 0.00036892361111111D0*(31.D0+2.D0*rk) * den(3)
ckplm(272,k) = -0.0129123263888889D0 * den(4)
ckplm(273,k) = -0.00036892361111111D0*(-29.D0+2.D0*rk) * den(5)
ckplm(274,k) = 0.00036892361111111D0*(-17.D0+2.D0*rk) * den(6)
ckplm(275,k) = -0.00005270337301587D0*(-43.D0+6.D0*rk) * den(7)
ckplm(276,k) = 0.00005270337301587D0*(rk-7.D0) * den(8)
ckplm(277,k) = 3.29396081349206D-6 * den(0)
ckplm(278,k) = -0.00002635168650794D0 * den(1)
ckplm(279,k) = 0.00009223090277778D0 * den(2)
ckplm(280,k) = -0.00018446180555556D0 * den(3)
ckplm(281,k) = 0.00023057725694444D0 * den(4)
ckplm(282,k) = -0.00018446180555556D0 * den(5)
ckplm(283,k) = 0.00009223090277778D0 * den(6)
ckplm(284,k) = -0.00002635168650794D0 * den(7)
ckplm(285,k) = 3.29396081349206D-6 * den(8)
* ckplm para l = 9
den(0) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k+17.D0)*(r2
&k+19.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(1) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k+17.D0)*(r2
&k-1.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(2) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k-1.D0)*(r2k
&+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(3) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-
&3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k+9.D0))
den(4) = 1.d0/((r2k+11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k+9.D0))
den(5) = 1.d0/((r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.
&D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k-9.D0)*(r2k+9.D0))
den(6) = 1.d0/((r2k-11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k-9.D0))
den(7) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-
&3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k-9.D0))
den(8) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-15.D0)*(r2k-1.D0)*(r2k
&+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k-7.D0)*(r2k-9.D0))
den(9) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-15.D0)*(r2k-17.D0)*(r2
&k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k-5.D0)*(r2k-7.D0)*(r2k-9.D0))
ckplm(286,k) = 1804.2578125D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.
&D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(287,k) = 955.1953125D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D
&0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*rk * den(1)
ckplm(288,k) = 764.15625D0*(rk-1.D0)*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)
&*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*rk * den(2)
ckplm(289,k) = 685.78125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)
&*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*rk * den(3)
ckplm(290,k) = 654.609375D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk-3.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*rk * den(4)
ckplm(291,k) = 654.609375D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk+4.D0)*rk * den(5)
ckplm(292,k) = 685.78125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)
&*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk-5.D0)*rk * den(6)
ckplm(293,k) = 764.15625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0)
&*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*rk * den(7)
ckplm(294,k) = 955.1953125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk-3.D
&0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*rk * den(8)
ckplm(295,k) = 1804.2578125D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.
&D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)*rk * den(9)
ckplm(296,k) = -360.8515625D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.
&D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(297,k) = -21.2265625D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D
&0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(-17.D0+7.D0*rk) * den(1)
ckplm(298,k) = -84.90625D0*(rk-1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)
&*(rk+5.D0)*(rk-6.D0)*(rk+6.D0)*(rk+7.D0) * den(2)
ckplm(299,k) = -45.71875D0*(rk-13.D0)*(rk-1.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0) * den(3)
ckplm(300,k) = -14.546875D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0
&)*(rk+3.D0)*(rk-44.D0)*(rk+4.D0)*(rk+5.D0) * den(4)
ckplm(301,k) = 14.546875D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0)
&*(rk+3.D0)*(rk+45.D0)*(rk-4.D0)*(rk+4.D0) * den(5)
ckplm(302,k) = 45.71875D0*(rk+14.D0)*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)
&*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk-5.D0) * den(6)
ckplm(303,k) = 84.90625D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0)*
&(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk+7.D0) * den(7)
ckplm(304,k) = 21.2265625D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0
&)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(24.D0+7.D0*rk) * den(8)
ckplm(305,k) = 360.8515625D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D
&0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0) * den(9)
ckplm(306,k) = 32.8046875D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0
&)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(307,k) = 1.9296875D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk-68.D0
&)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0) * den(1)
ckplm(308,k) = -3.859375D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)
&*(rk+7.D0)*(-87.D0+2.D0*(rk+17.D0)*rk) * den(2)
ckplm(309,k) = -2.078125D0*(6.D0*rk*rk+46.D0*rk-247.D0)*(rk-2.D0)*
&(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0) * den(3)
ckplm(310,k) = -14.546875D0*(rk-2.D0)*(rk-3.D0)*(rk+3.D0)*(rk+4.D0
&)*(rk+5.D0)*(-43.D0+(rk+3.D0)*rk) * den(4)
ckplm(311,k) = -14.546875D0*(rk-2.D0)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0
&)*(rk+4.D0)*(-45.D0+(rk-1.D0)*rk) * den(5)
ckplm(312,k) = -2.078125D0*(6.D0*rk*rk-34.D0*rk-287.D0)*(rk-2.D0)*
&(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk-5.D0) * den(6)
ckplm(313,k) = -3.859375D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)
&*(rk-6.D0)*(-119.D0+2.D0*(rk-15.D0)*rk) * den(7)
ckplm(314,k) = 1.9296875D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)
&*(rk+69.D0)*(rk-6.D0)*(rk-7.D0) * den(8)
ckplm(315,k) = 32.8046875D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0
&)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0) * den(9)
ckplm(316,k) = -2.73372395833333D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(
&rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(317,k) = 1.447265625D0*(rk+17.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.
&D0)*(rk+7.D0)*(rk+8.D0) * den(1)
ckplm(318,k) = 1.9296875D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)
&*(-57.D0+(rk+5.D0)*rk) * den(2)
ckplm(319,k) = 0.0494791666666667D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*
&(8931.D0+rk*(-2516.D0+rk*(-204.D0+29.D0*rk))) * den(3)
ckplm(320,k) = 0.51953125D0*(rk-3.D0)*(rk+4.D0)*(rk+5.D0)*(1176.D0
&+rk*(-124.D0+(rk-39.D0)*rk)) * den(4)
ckplm(321,k) = -0.51953125D0*(rk-3.D0)*(rk-4.D0)*(rk+4.D0)*(-1260.
&D0+(rk-1.D0)*(rk+43.D0)*rk) * den(5)
ckplm(322,k) = -0.0494791666666667D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)
&*(-11214.D0+rk*(-2021.D0+rk*(291.D0+29.D0*rk))) * den(6)
ckplm(323,k) = -1.9296875D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0
&)*(-61.D0+(rk-3.D0)*rk) * den(7)
ckplm(324,k) = -1.447265625D0*(rk-16.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5
&.D0)*(rk-6.D0)*(rk-7.D0) * den(8)
ckplm(325,k) = 2.73372395833333D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(r
&k-6.D0)*(rk-7.D0)*(rk-8.D0) * den(9)
ckplm(326,k) = 0.210286458333333D0*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(
&rk+8.D0)*(rk+9.D0) * den(0)
ckplm(327,k) = -0.0123697916666667D0*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)
&*(rk+8.D0)*(272.D0+23.D0*rk) * den(1)
ckplm(328,k) = -0.1484375D0*(rk-18.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D
&0)*(rk+9.D0) * den(2)
ckplm(329,k) = 0.0494791666666667D0*(rk+5.D0)*(rk+6.D0)*(-2556.D0+
&rk*(116.D0+(rk+69.D0)*rk)) * den(3)
ckplm(330,k) = 0.173177083333333D0*(rk+5.D0)*(3444.D0+rk*(-490.D0+
&rk*(-145.D0+(rk+10.D0)*rk))) * den(4)
ckplm(331,k) = 0.173177083333333D0*(rk-4.D0)*(3780.D0+(rk-1.D0)*rk
&*(-174.D0+(rk-5.D0)*rk)) * den(5)
ckplm(332,k) = 0.0494791666666667D0*(rk-4.D0)*(rk-5.D0)*(2604.D0+r
&k*(-19.D0+(rk-66.D0)*rk)) * den(6)
ckplm(333,k) = -0.1484375D0*(rk+19.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D
&0)*(rk-8.D0) * den(7)
ckplm(334,k) = -0.0123697916666667D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)
&*(rk-7.D0)*(-249.D0+23.D0*rk) * den(8)
ckplm(335,k) = 0.210286458333333D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(
&rk-7.D0)*(rk-8.D0) * den(9)
ckplm(336,k) = -0.0150204613095238D0*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)
&*(rk+9.D0) * den(0)
ckplm(337,k) = 0.00088355654761905D0*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)
&*(425.D0+41.D0*rk) * den(1)
ckplm(338,k) = -0.00212053571428571D0*(rk+6.D0)*(rk+7.D0)*(1875.D0
&+rk*(247.D0+3.D0*rk)) * den(2)
ckplm(339,k) = -0.0247395833333333D0*(rk+15.D0)*(rk+6.D0)*(-73.D0+
&(rk-3.D0)*rk) * den(3)
ckplm(340,k) = -0.0123697916666667D0*(11760.D0+rk*(824.D0+rk*(-379
&.D0+(rk-26.D0)*rk))) * den(4)
ckplm(341,k) = 0.0123697916666667D0*(10584.D0+rk*(-1500.D0+rk*(-29
&5.D0+(rk+30.D0)*rk))) * den(5)
ckplm(342,k) = 0.0247395833333333D0*(rk-14.D0)*(rk-5.D0)*(-69.D0+(
&rk+5.D0)*rk) * den(6)
ckplm(343,k) = 0.00212053571428571D0*(rk-5.D0)*(rk-6.D0)*(1631.D0+
&rk*(-241.D0+3.D0*rk)) * den(7)
ckplm(344,k) = -0.00088355654761905D0*(rk-5.D0)*(rk-6.D0)*(rk-7.D0
&)*(-384.D0+41.D0*rk) * den(8)
ckplm(345,k) = 0.0150204613095238D0*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*
&(rk-8.D0) * den(9)
ckplm(346,k) = 0.00100136408730159D0*(rk+7.D0)*(rk+8.D0)*(rk+9.D0)
& * den(0)
ckplm(347,k) = -0.00053013392857143D0*(rk+7.D0)*(rk+8.D0)*(68.D0+7
&.D0*rk) * den(1)
ckplm(348,k) = 0.00035342261904762D0*(rk+7.D0)*(1509.D0+2.D0*rk*(1
&33.D0+5.D0*rk)) * den(2)
ckplm(349,k) = 0.00082465277777778D0*(-5622.D0+rk*(-1193.D0+2.D0*(
&rk-21.D0)*rk)) * den(3)
ckplm(350,k) = -0.00247395833333333D0*(-1911.D0+rk*(-124.D0+(rk+24
&.D0)*rk)) * den(4)
ckplm(351,k) = -0.00247395833333333D0*(1764.D0+rk*(-169.D0+(rk-21.
&D0)*rk)) * den(5)
ckplm(352,k) = 0.00082465277777778D0*(4473.D0+rk*(-1103.D0+2.D0*(r
&k+24.D0)*rk)) * den(6)
ckplm(353,k) = 0.00035342261904762D0*(rk-6.D0)*(1253.D0+2.D0*rk*(-
&123.D0+5.D0*rk)) * den(7)
ckplm(354,k) = -0.00053013392857143D0*(rk-6.D0)*(rk-7.D0)*(-61.D0+
&7.D0*rk) * den(8)
ckplm(355,k) = 0.00100136408730159D0*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)
& * den(9)
ckplm(356,k) = -0.00006258525545635D0*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(357,k) = 3.68148561507937D-6*(rk+8.D0)*(833.D0+89.D0*rk) * d
&en(1)
ckplm(358,k) = -0.00001472594246032D0*(4116.D0+rk*(841.D0+41.D0*rk
&)) * den(2)
ckplm(359,k) = 0.00010308159722222D0*(916.D0+rk*(131.D0+3.D0*rk))
&* den(3)
ckplm(360,k) = 0.00036078559027778D0*(-304.D0+(rk-15.D0)*rk) * den
&(4)
ckplm(361,k) = -0.00036078559027778D0*(-288.D0+(rk+17.D0)*rk) * de
&n(5)
ckplm(362,k) = -0.00010308159722222D0*(788.D0+rk*(-125.D0+3.D0*rk)
&) * den(6)
ckplm(363,k) = 0.00001472594246032D0*(3316.D0+rk*(-759.D0+41.D0*rk
&)) * den(7)
ckplm(364,k) = -3.68148561507937D-6*(rk-7.D0)*(-744.D0+89.D0*rk) *
& den(8)
ckplm(365,k) = 0.00006258525545635D0*(rk-7.D0)*(rk-8.D0) * den(9)
ckplm(366,k) = 3.68148561507937D-6*(rk+9.D0) * den(0)
ckplm(367,k) = -3.68148561507937D-6*(64.D0+7.D0*rk) * den(1)
ckplm(368,k) = 0.00001472594246032D0*(51.D0+5.D0*rk) * den(2)
ckplm(369,k) = -0.00010308159722222D0*(rk+14.D0) * den(3)
ckplm(370,k) = 0.00005154079861111D0*(rk+37.D0) * den(4)
ckplm(371,k) = 0.00005154079861111D0*(rk-36.D0) * den(5)
ckplm(372,k) = -0.00010308159722222D0*(rk-13.D0) * den(6)
ckplm(373,k) = 0.00001472594246032D0*(-46.D0+5.D0*rk) * den(7)
ckplm(374,k) = -3.68148561507937D-6*(-57.D0+7.D0*rk) * den(8)
ckplm(375,k) = 3.68148561507937D-6*(rk-8.D0) * den(9)
ckplm(376,k) = -2.0452697861552D-7 * den(0)
ckplm(377,k) = 1.84074280753968D-6 * den(1)
ckplm(378,k) = -7.36297123015873D-6 * den(2)
ckplm(379,k) = 0.0000171802662037D0 * den(3)
ckplm(380,k) = -0.00002577039930556D0 * den(4)
ckplm(381,k) = 0.00002577039930556D0 * den(5)
ckplm(382,k) = -0.0000171802662037D0 * den(6)
ckplm(383,k) = 7.36297123015873D-6 * den(7)
ckplm(384,k) = -1.84074280753968D-6 * den(8)
ckplm(385,k) = 2.0452697861552D-7 * den(9)
* ckplm para l = 10
den(0) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k+17.D0)*(r2
&k+19.D0)*(r2k+1.D0)*(r2k+21.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(
&r2k+9.D0))
den(1) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k+17.D0)*(r2
&k+19.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r
&2k+9.D0))
den(2) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k+17.D0)*(r2
&k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2
&k+9.D0))
den(3) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k+15.D0)*(r2k-1.D0)*(r2k
&+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k+7.D0)*(r2k
&+9.D0))
den(4) = 1.d0/((r2k+11.D0)*(r2k+13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-
&3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k+
&9.D0))
den(5) = 1.d0/((r2k+11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k-9.D0)*(r2k+9
&.D0))
den(6) = 1.d0/((r2k-11.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3
&.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k-9.D0)*(r2k+9
&.D0))
den(7) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-
&3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k+7.D0)*(r2k-
&9.D0))
den(8) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-15.D0)*(r2k-1.D0)*(r2k
&+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k+5.D0)*(r2k-7.D0)*(r2k
&-9.D0))
den(9) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-15.D0)*(r2k-17.D0)*(r2
&k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k+3.D0)*(r2k-5.D0)*(r2k-7.D0)*(r2
&k-9.D0))
den(10) = 1.d0/((r2k-11.D0)*(r2k-13.D0)*(r2k-15.D0)*(r2k-17.D0)*(r
&2k-19.D0)*(r2k-1.D0)*(r2k+1.D0)*(r2k-3.D0)*(r2k-5.D0)*(r2k-7.D0)*(
&r2k-9.D0))
ckplm(386,k) = 3788.94140625D0*(rk+10.D0)*(rk+1.D0)*(rk+2.D0)*(rk+
&3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0)
&* den(0)
ckplm(387,k) = 1994.1796875D0*(rk+1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.
&D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0)*rk * den(1)
ckplm(388,k) = 1583.61328125D0*(rk-1.D0)*(rk+1.D0)*(rk+2.D0)*(rk+3
&.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*rk * den(2)
ckplm(389,k) = 1407.65625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*rk * den(3)
ckplm(390,k) = 1326.4453125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.
&D0)*(rk-3.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*rk * den(4)
ckplm(391,k) = 1302.328125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D
&0)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk+4.D0)*(rk+5.D0)*rk * den(5)
ckplm(392,k) = 1326.4453125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.
&D0)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk+4.D0)*(rk-5.D0)*rk * den(6)
ckplm(393,k) = 1407.65625D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*rk * den(7)
ckplm(394,k) = 1583.61328125D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk+2
&.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*rk * den(8)
ckplm(395,k) = 1994.1796875D0*(rk-1.D0)*(rk+1.D0)*(rk-2.D0)*(rk-3.
&D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)*rk * den(9)
ckplm(396,k) = 3788.94140625D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4
&.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)*(rk-9.D0)*rk * den(10
&)
ckplm(397,k) = -688.8984375D0*(rk+10.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4
&.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(398,k) = -36.2578125D0*(rk+2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D
&0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0)*(-19.D0+8.D0*rk) * den(
&1)
ckplm(399,k) = -57.5859375D0*(rk-1.D0)*(rk+2.D0)*(rk+3.D0)*(rk+4.D
&0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(-17.D0+3.D0*rk) * den(
&2)
ckplm(400,k) = -25.59375D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk+3.D0)
&*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(-45.D0+4.D0*rk) * den(3)
ckplm(401,k) = -48.234375D0*(rk-1.D0)*(rk-26.D0)*(rk-2.D0)*(rk+2.D
&0)*(rk-3.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0) * den(4)
ckplm(402,k) = 1302.328125D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D
&0)*(rk+3.D0)*(rk-4.D0)*(rk+4.D0)*(rk+5.D0) * den(5)
ckplm(403,k) = 48.234375D0*(rk-1.D0)*(rk+27.D0)*(rk-2.D0)*(rk+2.D0
&)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)*(rk+4.D0)*(rk-5.D0) * den(6)
ckplm(404,k) = 25.59375D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0)*
&(rk+3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(49.D0+4.D0*rk) * den(7)
ckplm(405,k) = 57.5859375D0*(rk-1.D0)*(rk-2.D0)*(rk+2.D0)*(rk-3.D0
&)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(20.D0+3.D0*rk) * den(8)
ckplm(406,k) = 36.2578125D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D0
&)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)*(27.D0+8.D0*rk) * den(9)
ckplm(407,k) = 688.8984375D0*(rk-1.D0)*(rk-2.D0)*(rk-3.D0)*(rk-4.D
&0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)*(rk-9.D0) * den(10)
ckplm(408,k) = 57.408203125D0*(rk+10.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5
&.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(409,k) = 6.04296875D0*(rk-38.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D
&0)*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(1)
ckplm(410,k) = -0.533203125D0*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.
&D0)*(rk+7.D0)*(rk+8.D0)*(-1122.D0+rk*(457.D0+19.D0*rk)) * den(2)
ckplm(411,k) = -1.421875D0*(rk-2.D0)*(rk+3.D0)*(rk+4.D0)*(rk+5.D0)
&*(rk+6.D0)*(rk+7.D0)*(-660.D0+rk*(137.D0+13.D0*rk)) * den(3)
ckplm(412,k) = -1.33984375D0*(rk-2.D0)*(rk-3.D0)*(rk+3.D0)*(rk+4.D
&0)*(rk+5.D0)*(rk+6.D0)*(-884.D0+rk*(87.D0+17.D0*rk)) * den(4)
ckplm(413,k) = -24.1171875D0*(rk*rk+rk-54.D0)*(rk-2.D0)*(rk-3.D0)*
&(rk+3.D0)*(rk-4.D0)*(rk+4.D0)*(rk+5.D0) * den(5)
ckplm(414,k) = -1.33984375D0*(rk-2.D0)*(rk-3.D0)*(rk+3.D0)*(rk-4.D
&0)*(rk+4.D0)*(rk-5.D0)*(-954.D0+rk*(-53.D0+17.D0*rk)) * den(6)
ckplm(415,k) = -1.421875D0*(rk-2.D0)*(rk-3.D0)*(rk+3.D0)*(rk-4.D0)
&*(rk-5.D0)*(rk-6.D0)*(-784.D0+rk*(-111.D0+13.D0*rk)) * den(7)
ckplm(416,k) = -0.533203125D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.
&D0)*(rk-6.D0)*(rk-7.D0)*(-1560.D0+rk*(-419.D0+19.D0*rk)) * den(8)
ckplm(417,k) = 6.04296875D0*(rk-2.D0)*(rk+39.D0)*(rk-3.D0)*(rk-4.D
&0)*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0) * den(9)
ckplm(418,k) = 57.408203125D0*(rk-2.D0)*(rk-3.D0)*(rk-4.D0)*(rk-5.
&D0)*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)*(rk-9.D0) * den(10)
ckplm(419,k) = -4.416015625D0*(rk+10.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6
&.D0)*(rk+7.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(420,k) = 0.232421875D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D
&0)*(rk+8.D0)*(rk+9.D0)*(171.D0+8.D0*rk) * den(1)
ckplm(421,k) = 0.123046875D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D
&0)*(rk+8.D0)*(-1479.D0+rk*(164.D0+23.D0*rk)) * den(2)
ckplm(422,k) = 0.0546875D0*(rk+4.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)
&*(13935.D0+rk*(-4367.D0+rk*(-177.D0+44.D0*rk))) * den(3)
ckplm(423,k) = 1.33984375D0*(rk-3.D0)*(rk+4.D0)*(rk+5.D0)*(rk+6.D0
&)*(834.D0+rk*(-124.D0+(rk-21.D0)*rk)) * den(4)
ckplm(424,k) = -36.17578125D0*(rk*rk+rk-36.D0)*(rk-3.D0)*(rk-4.D0)
&*(rk+4.D0)*(rk+5.D0) * den(5)
ckplm(425,k) = -1.33984375D0*(rk-3.D0)*(rk-4.D0)*(rk+4.D0)*(rk-5.D
&0)*(-936.D0+rk*(-79.D0+(rk+24.D0)*rk)) * den(6)
ckplm(426,k) = -0.0546875D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D0
&)*(-18081.D0+rk*(-3881.D0+rk*(309.D0+44.D0*rk))) * den(7)
ckplm(427,k) = -0.123046875D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.
&D0)*(rk-7.D0)*(-1620.D0+rk*(-118.D0+23.D0*rk)) * den(8)
ckplm(428,k) = -0.232421875D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.
&D0)*(rk-7.D0)*(rk-8.D0)*(-163.D0+8.D0*rk) * den(9)
ckplm(429,k) = 4.416015625D0*(rk-3.D0)*(rk-4.D0)*(rk-5.D0)*(rk-6.D
&0)*(rk-7.D0)*(rk-8.D0)*(rk-9.D0) * den(10)
ckplm(430,k) = 0.3154296875D0*(rk+10.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7
&.D0)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(431,k) = -0.033203125D0*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8.
&D0)*(rk+9.D0)*(152.D0+11.D0*rk) * den(1)
ckplm(432,k) = -0.0029296875D0*(rk+5.D0)*(rk+6.D0)*(rk+7.D0)*(rk+8
&.D0)*(-12648.D0+rk*(-355.D0+83.D0*rk)) * den(2)
ckplm(433,k) = -0.0078125D0*(rk-5.D0)*(rk+5.D0)*(rk+6.D0)*(rk+7.D0
&)*(-5268.D0+(rk-598.D0)*rk) * den(3)
ckplm(434,k) = 0.013671875D0*(rk+5.D0)*(rk+6.D0)*(77028.D0+rk*(-15
&406.D0+rk*(-2209.D0+rk*(274.D0+13.D0*rk)))) * den(4)
ckplm(435,k) = 0.24609375D0*(rk-4.D0)*(rk+5.D0)*(5292.D0+(rk*rk+rk
&-198.D0)*(rk+1.D0)*rk) * den(5)
ckplm(436,k) = 0.013671875D0*(rk-4.D0)*(rk-5.D0)*(89964.D0+rk*(102
&18.D0+rk*(-2953.D0+rk*(-222.D0+13.D0*rk)))) * den(6)
ckplm(437,k) = -0.0078125D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk+6.D0
&)*(-4669.D0+(rk+600.D0)*rk) * den(7)
ckplm(438,k) = -0.0029296875D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7
&.D0)*(-12210.D0+rk*(521.D0+83.D0*rk)) * den(8)
ckplm(439,k) = -0.033203125D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.
&D0)*(rk-8.D0)*(-141.D0+11.D0*rk) * den(9)
ckplm(440,k) = 0.3154296875D0*(rk-4.D0)*(rk-5.D0)*(rk-6.D0)*(rk-7.
&D0)*(rk-8.D0)*(rk-9.D0) * den(10)
ckplm(441,k) = -0.0210286458333333D0*(rk+10.D0)*(rk+6.D0)*(rk+7.D0
&)*(rk+8.D0)*(rk+9.D0) * den(0)
ckplm(442,k) = 0.00553385416666667D0*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)
&*(rk+9.D0)*(95.D0+8.D0*rk) * den(1)
ckplm(443,k) = 0.0009765625D0*(rk+6.D0)*(rk+7.D0)*(rk+8.D0)*(-5865
&.D0+(rk-596.D0)*rk) * den(2)
ckplm(444,k) = -0.00651041666666667D0*(rk+6.D0)*(rk+7.D0)*(-6393.D
&0+rk*(-373.D0+rk*(81.D0+4.D0*rk))) * den(3)
ckplm(445,k) = -0.0227864583333333D0*(rk+6.D0)*(10878.D0+rk*(-22.D
&0+rk*(-325.D0+(rk-8.D0)*rk))) * den(4)
ckplm(446,k) = 0.123046875D0*(10584.D0+5.D0*(rk*rk+rk-100.D0)*(rk+
&1.D0)*rk) * den(5)
ckplm(447,k) = 0.0227864583333333D0*(rk-5.D0)*(10584.D0+rk*(-600.D
&0+rk*(-295.D0+(rk+12.D0)*rk))) * den(6)
ckplm(448,k) = 0.00651041666666667D0*(rk-5.D0)*(rk-6.D0)*(5943.D0+
&rk*(-523.D0+rk*(-69.D0+4.D0*rk))) * den(7)
ckplm(449,k) = 0.0009765625D0*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*(5268.
&D0-(rk+598.D0)*rk) * den(8)
ckplm(450,k) = -0.00553385416666667D0*(rk-5.D0)*(rk-6.D0)*(rk-7.D0
&)*(rk-8.D0)*(-87.D0+8.D0*rk) * den(9)
ckplm(451,k) = 0.0210286458333333D0*(rk-5.D0)*(rk-6.D0)*(rk-7.D0)*
&(rk-8.D0)*(rk-9.D0) * den(10)
ckplm(452,k) = 0.00131429036458333D0*(rk+10.D0)*(rk+7.D0)*(rk+8.D0
&)*(rk+9.D0) * den(0)
ckplm(453,k) = -0.00013834635416667D0*(rk+7.D0)*(rk+8.D0)*(rk+9.D0
&)*(342.D0+31.D0*rk) * den(1)
ckplm(454,k) = 0.00003662109375D0*(rk+7.D0)*(rk+8.D0)*(19686.D0+rk
&*(2845.D0+79.D0*rk)) * den(2)
ckplm(455,k) = 0.00003255208333333D0*(rk+7.D0)*(-206880.D0+rk*(-32
&234.D0+rk*(-159.D0+83.D0*rk))) * den(3)
ckplm(456,k) = -0.00005696614583333D0*(-818496.D0+rk*(-118998.D0+r
&k*(9533.D0+rk*(1542.D0+19.D0*rk)))) * den(4)
ckplm(457,k) = -0.003076171875D0*(14112.D0+(rk*rk+rk-338.D0)*(rk+1
&.D0)*rk) * den(5)
ckplm(458,k) = -0.00005696614583333D0*(-691488.D0+rk*(133514.D0+rk
&*(5021.D0+rk*(-1466.D0+19.D0*rk)))) * den(6)
ckplm(459,k) = 0.00003255208333333D0*(rk-6.D0)*(174888.D0+rk*(-316
&67.D0+rk*(408.D0+83.D0*rk))) * den(7)
ckplm(460,k) = 0.00003662109375D0*(rk-6.D0)*(rk-7.D0)*(16920.D0+rk
&*(-2687.D0+79.D0*rk)) * den(8)
ckplm(461,k) = -0.00013834635416667D0*(rk-6.D0)*(rk-7.D0)*(rk-8.D0
&)*(-311.D0+31.D0*rk) * den(9)
ckplm(462,k) = 0.00131429036458333D0*(rk-6.D0)*(rk-7.D0)*(rk-8.D0)
&*(rk-9.D0) * den(10)
ckplm(463,k) = -0.00007731119791667D0*(rk+10.D0)*(rk+8.D0)*(rk+9.D
&0) * den(0)
ckplm(464,k) = 4.06901041666667D-6*(rk+8.D0)*(rk+9.D0)*(931.D0+88.
&D0*rk) * den(1)
ckplm(465,k) = -0.00001220703125D0*(rk+8.D0)*(6321.D0+rk*(1084.D0+
&43.D0*rk)) * den(2)
ckplm(466,k) = 0.00001627604166667D0*(55860.D0+rk*(12503.D0+rk*(72
&3.D0+4.D0*rk))) * den(3)
ckplm(467,k) = 0.00039876302083333D0*(-2616.D0+rk*(-304.D0+(rk+9.D
&0)*rk)) * den(4)
ckplm(468,k) = -0.0107666015625D0*(rk*rk+rk-96.D0) * den(5)
ckplm(469,k) = -0.00039876302083333D0*(2304.D0+rk*(-319.D0+(rk-6.D
&0)*rk)) * den(6)
ckplm(470,k) = 0.00001627604166667D0*(44076.D0+rk*(-11069.D0+(711.
&D0-4.D0*rk)*rk)) * den(7)
ckplm(471,k) = 0.00001220703125D0*(rk-7.D0)*(5280.D0+rk*(-998.D0+4
&3.D0*rk)) * den(8)
ckplm(472,k) = 4.06901041666667D-6*(rk-7.D0)*(rk-8.D0)*(843.D0-88.
&D0*rk) * den(9)
ckplm(473,k) = 0.00007731119791667D0*(rk-7.D0)*(rk-8.D0)*(rk-9.D0)
& * den(10)
ckplm(474,k) = 4.29506655092593D-6*(rk+10.D0)*(rk+9.D0) * den(0)
ckplm(475,k) = -4.52112268518519D-7*(rk+9.D0)*(608.D0+59.D0*rk) *
&den(1)
ckplm(476,k) = 2.03450520833333D-6*(3552.D0+rk*(673.D0+31.D0*rk))
&* den(2)
ckplm(477,k) = 5.42534722222222D-6*(-2430.D0-rk*(359.D0+11.D0*rk))
& * den(3)
ckplm(478,k) = -9.49435763888889D-6*(-1842.D0+(rk-149.D0)*rk) * de
&n(4)
ckplm(479,k) = 0.00005696614583333D0*(rk*rk+rk-324.D0) * den(5)
ckplm(480,k) = -9.49435763888889D-6*(-1692.D0+(rk+151.D0)*rk) * de
&n(6)
ckplm(481,k) = 5.42534722222222D-6*(-2082.D0+(337.D0-11.D0*rk)*rk)
& * den(7)
ckplm(482,k) = 2.03450520833333D-6*(2910.D0+rk*(-611.D0+31.D0*rk))
& * den(8)
ckplm(483,k) = 4.52112268518519D-7*(rk-8.D0)*(549.D0-59.D0*rk) * d
&en(9)
ckplm(484,k) = 4.29506655092593D-6*(rk-8.D0)*(rk-9.D0) * den(10)
ckplm(485,k) = 2.26056134259259D-7*(-10.D0-rk) * den(0)
ckplm(486,k) = 2.26056134259259D-7*(81.D0+8.D0*rk) * den(1)
ckplm(487,k) = 6.103515625D-6*(-11.D0-rk) * den(2)
ckplm(488,k) = 2.71267361111111D-6*(55.D0+4.D0*rk) * den(3)
ckplm(489,k) = -9.49435763888889D-6*(rk+24.D0) * den(4)
ckplm(490,k) = 0.00025634765625D0 * den(5)
ckplm(491,k) = 9.49435763888889D-6*(rk-23.D0) * den(6)
ckplm(492,k) = 2.71267361111111D-6*(51.D0-4.D0*rk) * den(7)
ckplm(493,k) = 6.103515625D-6*(rk-10.D0) * den(8)
ckplm(494,k) = 2.26056134259259D-7*(73.D0-8.D0*rk) * den(9)
ckplm(495,k) = 2.26056134259259D-7*(rk-9.D0) * den(10)
ckplm(496,k) = 1.1302806712963D-8 * den(0)
ckplm(497,k) = -1.1302806712963D-7 * den(1)
ckplm(498,k) = 5.08626302083333D-7 * den(2)
ckplm(499,k) = -1.35633680555556D-6 * den(3)
ckplm(500,k) = 2.37358940972222D-6 * den(4)
ckplm(501,k) = -2.84830729166667D-6 * den(5)
ckplm(502,k) = 2.37358940972222D-6 * den(6)
ckplm(503,k) = -1.35633680555556D-6 * den(7)
ckplm(504,k) = 5.08626302083333D-7 * den(8)