| | 836 | |
| | 837 | (rt:deftest oct.carlson-rf.1d |
| | 838 | ;; Rf(0,2,1) = integrate(1/sqrt(1-s^4), s, 0 ,1) |
| | 839 | ;; = 1/4*beta(1/2,1/2) |
| | 840 | ;; = sqrt(%pi)/4*gamma(1/4)/gamma(3/4) |
| | 841 | (let ((rf (carlson-rf 0d0 2d0 1d0)) |
| | 842 | (true 1.31102877714605990523241979494d0)) |
| | 843 | (check-accuracy 53 rf true)) |
| | 844 | nil) |
| | 845 | |
| | 846 | (rt:deftest oct.carlson-rf.1q |
| | 847 | ;; Rf(0,2,1) = integrate(1/sqrt(1-s^4), s, 0 ,1) |
| | 848 | (let ((rf (carlson-rf #q0 #q2 #q1)) |
| | 849 | (true #q1.311028777146059905232419794945559706841377475715811581408410851900395q0)) |
| | 850 | (check-accuracy 212 rf true)) |
| | 851 | nil) |
| | 852 | |
| | 853 | (rt:deftest oct.carlson-rd.1d |
| | 854 | ;; Rd(0,2,1) = 3*integrate(s^2/sqrt(1-s^4), s, 0 ,1) |
| | 855 | ;; = 3*beta(3/4,1/2)/4 |
| | 856 | ;; = 3*sqrt(%pi)*gamma(3/4)/gamma(1/4) |
| | 857 | (let ((rd (carlson-rd 0d0 2d0 1d0)) |
| | 858 | (true 1.7972103521033883d0)) |
| | 859 | (check-accuracy 51 rd true)) |
| | 860 | nil) |
| | 861 | |
| | 862 | (rt:deftest oct.carlson-rd.1q |
| | 863 | (let ((rd (carlson-rd #q0 #q2 #q1)) |
| | 864 | (true #q1.797210352103388311159883738420485817340818994823477337395512429419599q0)) |
| | 865 | (check-accuracy 212 rd true)) |
| | 866 | nil) |
