subroutine design_00_001_3d ( x ) !*****************************************************************************80 ! !! design_00_001_3d() returns a 0-design of 1 point in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 07 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Parameters: ! ! Output, real ( kind = rk8 ) X(3,1), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,1) x(1:3,1:1) = reshape ( (/ & 1.0D+00, 0.0D+00, 0.0D+00 /), (/ 3, 1 /) ) return end subroutine design_01_002_3d ( x ) !*****************************************************************************80 ! !! DESIGN_01_002_3D returns a 1-design of 2 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 07 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,2), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,2) x(1:3,1:2) = reshape ( (/ & 1.0D+00, 0.0D+00, 0.0D+00, & -1.0D+00, 0.0D+00, 0.0D+00 /), (/ 3, 2 /) ) return end subroutine design_02_004_3d ( x ) !*****************************************************************************80 ! !! DESIGN_02_004_3D returns a 2-design of 4 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 07 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,4), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ), parameter :: a = 0.577350269189625763D+00 real ( kind = rk8 ) x(3,4) x(1:3,1:4) = reshape ( (/ & a, a, a, & a, -a, -a, & -a, a, -a, & -a, -a, a /), (/ 3, 4 /) ) return end subroutine design_03_006_3d ( x ) !*****************************************************************************80 ! !! DESIGN_03_006_3D returns a 3-design of 6 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 07 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,6), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,6) x(1:3,1:6) = reshape ( (/ & 1.0D+00, 0.0D+00, 0.0D+00, & -1.0D+00, 0.0D+00, 0.0D+00, & 0.0D+00, 1.0D+00, 0.0D+00, & 0.0D+00,-1.0D+00, 0.0D+00, & 0.0D+00, 0.0D+00, 1.0D+00, & 0.0D+00, 0.0D+00, -1.0D+00 /), (/ 3, 6 /) ) return end subroutine design_04_014_3d ( x ) !*****************************************************************************80 ! !! DESIGN_04_014_3D returns a 4-design of 14 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 07 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,14), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,14) x(1:3,1:14) = reshape ( (/ & 1.0D+00, 0.0D+00, 0.0D+00, & 0.5947189772040725D+00, 0.1776539926025823D+00, 0.7840589244857197D+00, & 0.5947189772040725D+00, -0.7678419429698292D+00, -0.2381765915652909D+00, & 0.5947189772040725D+00, 0.590187950367247D+00, -0.5458823329204288D+00, & -0.5947189772040725D+00, 0.1776539926025823D+00, -0.7840589244857197D+00, & -0.5947189772040725D+00, 0.590187950367247D+00, 0.5458823329204288D+00, & -0.5947189772040725D+00, -0.7678419429698292D+00, 0.2381765915652909D+00, & 0.3012536847870683D+00, 0.879474443923065D+00, 0.3684710570566285D+00, & 0.3012536847870683D+00, -0.7588425179318781D+00, 0.5774116818882528D+00, & 0.3012536847870683D+00, -0.1206319259911869D+00, -0.9458827389448813D+00, & -0.3012536847870683D+00, 0.879474443923065D+00, -0.3684710570566285D+00, & -0.3012536847870683D+00, -0.1206319259911869D+00, 0.9458827389448813D+00, & -0.3012536847870683D+00, -0.7588425179318781D+00, -0.5774116818882528D+00, & -1.0D+00, 0.0D+00, 0.0D+00 /), & (/ 3, 14 /) ) return end subroutine design_05_012_3d ( x ) !*****************************************************************************80 ! !! DESIGN_05_012_3D returns a 5-design of 12 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 07 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,12), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,12) x(1:3,1:12) = reshape ( (/ & 0.850650808352D+00, 0.0D+00, -0.525731112119D+00, & 0.525731112119D+00, -0.850650808352D+00, 0.000000000000D+00, & 0.0D+00, -0.525731112119D+00, 0.850650808352D+00, & 0.850650808352D+00, 0.0D+00, 0.525731112119D+00, & -0.525731112119D+00, -0.850650808352D+00, 0.0D+00, & 0.0D+00, 0.525731112119D+00, -0.850650808352D+00, & -0.850650808352D+00, 0.0D+00, -0.525731112119D+00, & -0.525731112119D+00, 0.850650808352D+00, 0.0D+00, & 0.0D+00, 0.525731112119D+00, 0.850650808352D+00, & -0.850650808352D+00, 0.0D+00, 0.525731112119D+00, & 0.525731112119D+00, 0.850650808352D+00, 0.0D+00, & 0.0D+00, -0.525731112119D+00, -0.850650808352D+00 /), & (/ 3, 12 /) ) return end subroutine design_06_026_3d ( x ) !*****************************************************************************80 ! !! DESIGN_06_026_3D returns a 6-design of 26 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,26), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,26) x(1:3,1:26) = reshape ( (/ & -.6156289423339210D+00, & .1008195633550640D+00, & .7815602478412730D+00, & .3951267761987170D+00, & .4827405208135130D+00, & .7815602473869220D+00, & .2205021660506020D+00, & -.5835600880219350D+00, & .7815602461965580D+00, & -.4164742585347590D+00, & .6658869344429950D+00, & .6189860923445970D+00, & -.3684378761516490D+00, & -.6936207556288800D+00, & .6189860893248600D+00, & .7849121331270050D+00, & .2773382046045800D-01, & .6189860890784880D+00, & -.9753444939429849D+00, & .2176679667586930D+00, & .3638371039680900D-01, & .6761782314481880D+00, & .7358391298492950D+00, & .3638370924238100D-01, & .2991662611656320D+00, & -.9535070916939770D+00, & .3638370884709700D-01, & .7506752044088399D+00, & -.5662857931564620D+00, & .3403044783056840D+00, & .1150802794034220D+00, & .9332466938773990D+00, & .3403044778714140D+00, & -.8657554869972310D+00, & -.3669608990989080D+00, & .3403044743559930D+00, & .2031526000000000D-08, & -.1957257000000000D-08, & .1000000000000000D+01, & -.2031525000000000D-08, & -.1957257000000000D-08, & -.1000000000000000D+01, & .8657554869972310D+00, & -.3669608990989070D+00, & -.3403044743559920D+00, & -.1150802794034210D+00, & .9332466938773980D+00, & -.3403044778714150D+00, & -.7506752044088399D+00, & -.5662857931564620D+00, & -.3403044783056840D+00, & -.2991662611656320D+00, & -.9535070916939770D+00, & -.3638370884709700D-01, & -.6761782314481870D+00, & .7358391298492950D+00, & -.3638370924238000D-01, & .9753444939429849D+00, & .2176679667586940D+00, & -.3638371039681000D-01, & -.7849121331270040D+00, & .2773382046045900D-01, & -.6189860890784890D+00, & .3684378761516490D+00, & -.6936207556288810D+00, & -.6189860893248590D+00, & .4164742585347590D+00, & .6658869344429950D+00, & -.6189860923445970D+00, & -.2205021660506030D+00, & -.5835600880219350D+00, & -.7815602461965580D+00, & -.3951267761987180D+00, & .4827405208135130D+00, & -.7815602473869220D+00, & .6156289423339210D+00, & .1008195633550640D+00, & -.7815602478412730D+00 /), (/ 3, 26 /) ) return end subroutine design_07_024_3d ( x ) !*****************************************************************************80 ! !! DESIGN_07_024_3D returns a 7-design of 24 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,24), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,24) x(1:3,1:24) = reshape ( (/ & 0.8662468181078205913835980D+00, 0.4225186537611115291185464D+00, & 0.2666354015167047203315344D+00, & 0.8662468181078205913835980D+00, -0.4225186537611115291185464D+00, & -0.2666354015167047203315344D+00, & 0.8662468181078205913835980D+00, 0.2666354015167047203315344D+00, & -0.4225186537611115291185464D+00, & 0.8662468181078205913835980D+00, -0.2666354015167047203315344D+00, & 0.4225186537611115291185464D+00, & -0.8662468181078205913835980D+00, 0.4225186537611115291185464D+00, & -0.2666354015167047203315344D+00, & -0.8662468181078205913835980D+00, -0.4225186537611115291185464D+00, & 0.2666354015167047203315344D+00, & -0.8662468181078205913835980D+00, 0.2666354015167047203315344D+00, & 0.4225186537611115291185464D+00, & -0.8662468181078205913835980D+00, -0.2666354015167047203315344D+00, & -0.4225186537611115291185464D+00, & 0.2666354015167047203315344D+00, 0.8662468181078205913835980D+00, & 0.4225186537611115291185464D+00, & -0.2666354015167047203315344D+00, 0.8662468181078205913835980D+00, & -0.4225186537611115291185464D+00, & -0.4225186537611115291185464D+00, 0.8662468181078205913835980D+00, & 0.2666354015167047203315344D+00, & 0.4225186537611115291185464D+00, 0.8662468181078205913835980D+00, & -0.2666354015167047203315344D+00, & -0.2666354015167047203315344D+00, -0.8662468181078205913835980D+00, & 0.4225186537611115291185464D+00, & 0.2666354015167047203315344D+00, -0.8662468181078205913835980D+00, & -0.4225186537611115291185464D+00, & 0.4225186537611115291185464D+00, -0.8662468181078205913835980D+00, & 0.2666354015167047203315344D+00, & -0.4225186537611115291185464D+00, -0.8662468181078205913835980D+00, & -0.2666354015167047203315344D+00, & 0.4225186537611115291185464D+00, 0.2666354015167047203315344D+00, & 0.8662468181078205913835980D+00, & -0.4225186537611115291185464D+00, -0.2666354015167047203315344D+00, & 0.8662468181078205913835980D+00, & 0.2666354015167047203315344D+00, -0.4225186537611115291185464D+00, & 0.8662468181078205913835980D+00, & -0.2666354015167047203315344D+00, 0.4225186537611115291185464D+00, & 0.8662468181078205913835980D+00, & 0.4225186537611115291185464D+00, -0.2666354015167047203315344D+00, & -0.8662468181078205913835980D+00, & -0.4225186537611115291185464D+00, 0.2666354015167047203315344D+00, & -0.8662468181078205913835980D+00, & 0.2666354015167047203315344D+00, 0.4225186537611115291185464D+00, & -0.8662468181078205913835980D+00, & -0.2666354015167047203315344D+00, -0.4225186537611115291185464D+00, & -0.8662468181078205913835980D+00 & /), (/ 3, 24 /) ) return end subroutine design_08_036_3d ( x ) !*****************************************************************************80 ! !! DESIGN_08_036_3D returns an 8-design of 36 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,36), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,36) x(1:3,1:36) = reshape ( (/ & 0.5074754464108170D+00, -0.3062000132395710D+00, 0.8054254920116630D+00, & -0.3062000132395690D+00, 0.8054254920116630D+00, 0.5074754464108170D+00, & -0.5074754464108170D+00, 0.3062000132395700D+00, 0.8054254920116630D+00, & 0.8054254920116630D+00, 0.5074754464108170D+00, -0.3062000132395690D+00, & 0.3062000132395690D+00, 0.8054254920116640D+00, -0.5074754464108170D+00, & 0.8054254920116630D+00, -0.5074754464108170D+00, 0.3062000132395690D+00, & 0.3062000132395690D+00, -0.8054254920116630D+00, 0.5074754464108160D+00, & -0.8054254920116630D+00, -0.5074754464108170D+00, -0.3062000132395690D+00, & -0.3062000132395700D+00, -0.8054254920116640D+00, -0.5074754464108160D+00, & -0.8054254920116630D+00, 0.5074754464108180D+00, 0.3062000132395690D+00, & 0.5074754464108170D+00, 0.3062000132395700D+00, -0.8054254920116630D+00, & -0.5074754464108170D+00, -0.3062000132395700D+00, -0.8054254920116630D+00, & 0.6263636702652710D+00, -0.2435277754091940D+00, -0.7405152092807200D+00, & -0.2435277754091950D+00, -0.7405152092807200D+00, 0.6263636702652710D+00, & -0.6263636702652710D+00, 0.2435277754091940D+00, -0.7405152092807200D+00, & -0.7405152092807200D+00, 0.6263636702652700D+00, -0.2435277754091950D+00, & 0.2435277754091950D+00, -0.7405152092807190D+00, -0.6263636702652710D+00, & -0.7405152092807200D+00, -0.6263636702652700D+00, 0.2435277754091950D+00, & 0.2435277754091950D+00, 0.7405152092807190D+00, 0.6263636702652710D+00, & 0.7405152092807200D+00, -0.6263636702652700D+00, -0.2435277754091950D+00, & -0.2435277754091950D+00, 0.7405152092807200D+00, -0.6263636702652710D+00, & 0.7405152092807200D+00, 0.6263636702652700D+00, 0.2435277754091950D+00, & 0.6263636702652710D+00, 0.2435277754091940D+00, 0.7405152092807200D+00, & -0.6263636702652710D+00, -0.2435277754091940D+00, 0.7405152092807200D+00, & -0.2862487234260350D+00, 0.9571203270924580D+00, -0.4452356458542100D-01, & 0.9571203270924580D+00, -0.4452356458542000D-01, -0.2862487234260350D+00, & 0.2862487234260350D+00, -0.9571203270924580D+00, -0.4452356458542100D-01, & -0.4452356458542000D-01, -0.2862487234260350D+00, 0.9571203270924580D+00, & -0.9571203270924580D+00, -0.4452356458541900D-01, 0.2862487234260350D+00, & -0.4452356458542100D-01, 0.2862487234260340D+00, -0.9571203270924580D+00, & -0.9571203270924580D+00, 0.4452356458542000D-01, -0.2862487234260340D+00, & 0.4452356458542100D-01, 0.2862487234260340D+00, 0.9571203270924580D+00, & 0.9571203270924580D+00, 0.4452356458542000D-01, 0.2862487234260340D+00, & 0.4452356458542100D-01, -0.2862487234260340D+00, -0.9571203270924580D+00, & -0.2862487234260340D+00, -0.9571203270924580D+00, 0.4452356458542100D-01, & 0.2862487234260350D+00, 0.9571203270924580D+00, 0.4452356458542100D-01 & /), (/ 3, 36 /) ) return end subroutine design_09_048_3d ( x ) !*****************************************************************************80 ! !! DESIGN_09_048_3D returns a 9-design of 48 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,48), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,48) x(1:3,1:48) = reshape ( (/ & 0.9333646931993099D+00, 0.3535421889214720D+00, -0.06195377423185971D+00, & 0.9333646931993099D+00, -0.3535421889214720D+00, 0.06195377423185971D+00, & 0.9333646931993099D+00, -0.06195377423185971D+00, -0.3535421889214720D+00, & 0.9333646931993099D+00, 0.06195377423185971D+00, 0.3535421889214720D+00, & -0.9333646931993099D+00, 0.3535421889214720D+00, 0.06195377423185971D+00, & -0.9333646931993099D+00, -0.3535421889214720D+00, -0.06195377423185971D+00, & -0.9333646931993099D+00, -0.06195377423185971D+00, 0.3535421889214720D+00, & -0.9333646931993099D+00, 0.06195377423185971D+00, -0.3535421889214720D+00, & -0.06195377423185971D+00, 0.9333646931993099D+00, 0.3535421889214720D+00, & 0.06195377423185971D+00, 0.9333646931993099D+00, -0.3535421889214720D+00, & -0.3535421889214720D+00, 0.9333646931993099D+00, -0.06195377423185971D+00, & 0.3535421889214720D+00, 0.9333646931993099D+00, 0.06195377423185971D+00, & 0.06195377423185971D+00, -0.9333646931993099D+00, 0.3535421889214720D+00, & -0.06195377423185971D+00, -0.9333646931993099D+00, -0.3535421889214720D+00, & 0.3535421889214720D+00, -0.9333646931993099D+00, -0.06195377423185971D+00, & -0.3535421889214720D+00, -0.9333646931993099D+00, 0.06195377423185971D+00, & 0.3535421889214720D+00, -0.06195377423185971D+00, 0.9333646931993099D+00, & -0.3535421889214720D+00, 0.06195377423185971D+00, 0.9333646931993099D+00, & -0.06195377423185971D+00, -0.3535421889214720D+00, 0.9333646931993099D+00, & 0.06195377423185971D+00, 0.3535421889214720D+00, 0.9333646931993099D+00, & 0.3535421889214720D+00, 0.06195377423185971D+00, -0.9333646931993099D+00, & -0.3535421889214720D+00, -0.06195377423185971D+00, -0.9333646931993099D+00, & -0.06195377423185971D+00, 0.3535421889214720D+00, -0.9333646931993099D+00, & 0.06195377423185971D+00, -0.3535421889214720D+00, -0.9333646931993099D+00, & 0.7068416977125503D+00, 0.6397400986197920D+00, 0.3018400579657685D+00, & 0.7068416977125503D+00, -0.6397400986197920D+00, -0.3018400579657685D+00, & 0.7068416977125503D+00, 0.3018400579657685D+00, -0.6397400986197920D+00, & 0.7068416977125503D+00, -0.3018400579657685D+00, 0.6397400986197920D+00, & -0.7068416977125503D+00, 0.6397400986197920D+00, -0.3018400579657685D+00, & -0.7068416977125503D+00, -0.6397400986197920D+00, 0.3018400579657685D+00, & -0.7068416977125503D+00, 0.3018400579657685D+00, 0.6397400986197920D+00, & -0.7068416977125503D+00, -0.3018400579657685D+00, -0.6397400986197920D+00, & 0.3018400579657685D+00, 0.7068416977125503D+00, 0.6397400986197920D+00, & -0.3018400579657685D+00, 0.7068416977125503D+00, -0.6397400986197920D+00, & -0.6397400986197920D+00, 0.7068416977125503D+00, 0.3018400579657685D+00, & 0.6397400986197920D+00, 0.7068416977125503D+00, -0.3018400579657685D+00, & -0.3018400579657685D+00, -0.7068416977125503D+00, 0.6397400986197920D+00, & 0.3018400579657685D+00, -0.7068416977125503D+00, -0.6397400986197920D+00, & 0.6397400986197920D+00, -0.7068416977125503D+00, 0.3018400579657685D+00, & -0.6397400986197920D+00, -0.7068416977125503D+00, -0.3018400579657685D+00, & 0.6397400986197920D+00, 0.3018400579657685D+00, 0.7068416977125503D+00, & -0.6397400986197920D+00, -0.3018400579657685D+00, 0.7068416977125503D+00, & 0.3018400579657685D+00, -0.6397400986197920D+00, 0.7068416977125503D+00, & -0.3018400579657685D+00, 0.6397400986197920D+00, 0.7068416977125503D+00, & 0.6397400986197920D+00, -0.3018400579657685D+00, -0.7068416977125503D+00, & -0.6397400986197920D+00, 0.3018400579657685D+00, -0.7068416977125503D+00, & 0.3018400579657685D+00, 0.6397400986197920D+00, -0.7068416977125503D+00, & -0.3018400579657685D+00, -0.6397400986197920D+00, -0.7068416977125503D+00 & /), (/ 3, 48 /) ) return end subroutine design_10_060_3d ( x ) !*****************************************************************************80 ! !! DESIGN_10_060_3D returns a 10-design of 60 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,60), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,60) x(1:3,1:60) = reshape ( (/ & -0.7538286671970170D+00, 0.5459519080612600D+00,-0.3656211900262870D+00, & 0.5459519080612580D+00,-0.3656211900262900D+00,-0.7538286671970170D+00, & 0.7538286671970160D+00,-0.5459519080612610D+00,-0.3656211900262880D+00, & -0.3656211900262890D+00,-0.7538286671970170D+00, 0.5459519080612590D+00, & -0.5459519080612580D+00,-0.3656211900262880D+00, 0.7538286671970180D+00, & -0.3656211900262890D+00, 0.7538286671970170D+00,-0.5459519080612590D+00, & -0.5459519080612580D+00, 0.3656211900262890D+00,-0.7538286671970170D+00, & 0.3656211900262870D+00, 0.7538286671970170D+00, 0.5459519080612600D+00, & 0.5459519080612590D+00, 0.3656211900262890D+00, 0.7538286671970170D+00, & 0.3656211900262870D+00,-0.7538286671970180D+00,-0.5459519080612590D+00, & -0.7538286671970170D+00,-0.5459519080612610D+00, 0.3656211900262880D+00, & 0.7538286671970160D+00, 0.5459519080612610D+00, 0.3656211900262870D+00, & 0.7001810193637300D+00,-0.7131510658477930D+00, 0.3408954976125600D-01, & -0.7131510658477940D+00, 0.3408954976125400D-01, 0.7001810193637290D+00, & -0.7001810193637300D+00, 0.7131510658477930D+00, 0.3408954976125600D-01, & 0.3408954976125500D-01, 0.7001810193637300D+00,-0.7131510658477930D+00, & 0.7131510658477930D+00, 0.3408954976125400D-01,-0.7001810193637300D+00, & 0.3408954976125700D-01,-0.7001810193637290D+00, 0.7131510658477940D+00, & 0.7131510658477940D+00,-0.3408954976125500D-01, 0.7001810193637280D+00, & -0.3408954976125600D-01,-0.7001810193637290D+00,-0.7131510658477940D+00, & -0.7131510658477940D+00,-0.3408954976125400D-01,-0.7001810193637290D+00, & -0.3408954976125700D-01, 0.7001810193637290D+00, 0.7131510658477940D+00, & 0.7001810193637300D+00, 0.7131510658477930D+00,-0.3408954976125700D-01, & -0.7001810193637290D+00,-0.7131510658477940D+00,-0.3408954976125700D-01, & 0.2762302182617920D+00, 0.7705072072573600D-01,-0.9579979399532590D+00, & 0.7705072072573500D-01,-0.9579979399532580D+00, 0.2762302182617930D+00, & -0.2762302182617920D+00,-0.7705072072573400D-01,-0.9579979399532590D+00, & -0.9579979399532590D+00, 0.2762302182617910D+00, 0.7705072072573800D-01, & -0.7705072072573500D-01,-0.9579979399532590D+00,-0.2762302182617920D+00, & -0.9579979399532580D+00,-0.2762302182617930D+00,-0.7705072072573600D-01, & -0.7705072072573600D-01, 0.9579979399532580D+00, 0.2762302182617940D+00, & 0.9579979399532590D+00,-0.2762302182617900D+00, 0.7705072072573700D-01, & 0.7705072072573400D-01, 0.9579979399532590D+00,-0.2762302182617920D+00, & 0.9579979399532580D+00, 0.2762302182617930D+00,-0.7705072072573800D-01, & 0.2762302182617930D+00,-0.7705072072573600D-01, 0.9579979399532580D+00, & -0.2762302182617910D+00, 0.7705072072573500D-01, 0.9579979399532590D+00, & 0.4518191025552430D+00,-0.7833559375218190D+00, 0.4268641162190700D+00, & -0.7833559375218180D+00, 0.4268641162190710D+00, 0.4518191025552430D+00, & -0.4518191025552430D+00, 0.7833559375218190D+00, 0.4268641162190700D+00, & 0.4268641162190710D+00, 0.4518191025552420D+00,-0.7833559375218190D+00, & 0.7833559375218180D+00, 0.4268641162190700D+00,-0.4518191025552440D+00, & 0.4268641162190720D+00,-0.4518191025552420D+00, 0.7833559375218180D+00, & 0.7833559375218190D+00,-0.4268641162190700D+00, 0.4518191025552420D+00, & -0.4268641162190720D+00,-0.4518191025552410D+00,-0.7833559375218190D+00, & -0.7833559375218180D+00,-0.4268641162190700D+00,-0.4518191025552430D+00, & -0.4268641162190720D+00, 0.4518191025552410D+00, 0.7833559375218190D+00, & 0.4518191025552430D+00, 0.7833559375218180D+00,-0.4268641162190710D+00, & -0.4518191025552420D+00,-0.7833559375218190D+00,-0.4268641162190710D+00, & -0.3385843599592600D+00,-0.9332100372395270D+00, 0.1203314488667840D+00, & -0.9332100372395260D+00, 0.1203314488667870D+00,-0.3385843599592600D+00, & 0.3385843599592610D+00, 0.9332100372395260D+00, 0.1203314488667860D+00, & 0.1203314488667850D+00,-0.3385843599592610D+00,-0.9332100372395260D+00, & 0.9332100372395260D+00, 0.1203314488667890D+00, 0.3385843599592600D+00, & 0.1203314488667850D+00, 0.3385843599592610D+00, 0.9332100372395260D+00, & 0.9332100372395260D+00,-0.1203314488667870D+00,-0.3385843599592620D+00, & -0.1203314488667840D+00, 0.3385843599592620D+00,-0.9332100372395260D+00, & -0.9332100372395260D+00,-0.1203314488667870D+00, 0.3385843599592610D+00, & -0.1203314488667840D+00,-0.3385843599592620D+00, 0.9332100372395260D+00, & -0.3385843599592620D+00, 0.9332100372395260D+00,-0.1203314488667840D+00, & 0.3385843599592610D+00,-0.9332100372395270D+00,-0.1203314488667830D+00 & /), (/ 3, 60 /) ) return end subroutine design_11_070_3d ( x ) !*****************************************************************************80 ! !! DESIGN_11_070_3D returns an 11-design of 70 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,70), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,70) x(1:3,1:70) = reshape ( (/ & -0.6749405204804370D+00, 0.7256290520645010D+00, 0.1338572844994640D+00, & 0.9672433446143000D-01, -0.9103273829899871D+00, -0.4024282034122290D+00, & 0.9069603159163580D+00, 0.1351270221350530D+00, 0.3989532218717040D+00, & -0.1327587047580260D+00, -0.3076585240607330D+00, 0.9421896618429551D+00, & -0.2260558011275870D+00, -0.9588311747087040D+00, -0.1718765637988270D+00, & 0.2757382640198530D+00, -0.1806927335075380D+00, -0.9440966824498920D+00, & 0.8308816505135890D+00, 0.3332786445281770D+00, -0.4456018715639280D+00, & -0.6164713286127870D+00, -0.2675443371664000D+00, 0.7405289519313720D+00, & 0.4302772932874360D+00, -0.8926444716153570D+00, -0.1343402329005700D+00, & -0.6909871985230760D+00, 0.1751093390532070D+00, 0.7013368740153190D+00, & 0.8105170415355070D+00, -0.3814493375472150D+00, 0.4444755654311270D+00, & -0.8673444385462600D-01, -0.7060085178359240D+00, -0.7028720431147840D+00, & 0.8713208520567370D+00, 0.4604578060039600D+00, 0.1696425113618090D+00, & -0.6007352667495490D+00, 0.3032661185525090D+00, -0.7396937208206140D+00, & -0.8991009470834190D+00, -0.4180812466398280D+00, 0.1297133692484600D+00, & 0.8969270870795710D+00, -0.1880663273448430D+00, -0.4001910256139910D+00, & 0.1504949609669910D+00, 0.9030721531392540D+00, 0.4022585647913240D+00, & 0.2486017164026210D+00, -0.2242836122819530D+00, 0.9422812997525900D+00, & 0.8425846747084230D+00, -0.5107563820855460D+00, -0.1708185707275000D+00, & 0.2600345004183370D+00, 0.2093564899576840D+00, -0.9426303192157490D+00, & -0.5880246157243400D-01, 0.8945952131887460D+00, -0.4429917324880950D+00, & 0.6161176918013200D-01, -0.6712901087901590D+00, 0.7386295280714080D+00, & 0.9823375360976140D+00, 0.1337840147101790D+00, -0.1308235551485130D+00, & -0.3822775825325760D+00, -0.6052438479001370D+00, 0.6982433203920290D+00, & 0.6118392782163570D+00, 0.6515716084972490D+00, 0.4484497035699710D+00, & 0.6468653485695820D+00, -0.2984641292976520D+00, -0.7017723165974470D+00, & -0.1692010168812820D+00, 0.9704309127468180D+00, 0.1721477838129720D+00, & -0.4717254508623250D+00, -0.4752957036627900D+00, -0.7426769776211120D+00, & 0.1193697559557230D+00, -0.9846926044113470D+00, 0.1270091972286680D+00, & 0.4572892122317290D+00, 0.7961559905587140D+00, -0.3962602870260380D+00, & -0.8136314363509790D+00, 0.4209422727934990D+00, 0.4010130780372200D+00, & 0.2871545553868710D+00, 0.1641733239706600D+00, 0.9437104328219510D+00, & 0.7466675770451550D+00, 0.6440359890663980D+00, -0.1664487133527440D+00, & -0.1157796447409060D+00, 0.3149524646461050D+00, -0.9420191181058980D+00, & -0.8675792121114660D+00, 0.2219163150406650D+00, -0.4450387172267380D+00, & 0.6551400224339120D+00, -0.1511626316805080D+00, 0.7402306463452570D+00, & 0.1767365123580470D+00, 0.9760026717210610D+00, -0.1272123814448300D+00, & 0.4552846077010780D+00, -0.5527863541042300D+00, 0.6979564260801880D+00, & -0.4320239302197420D+00, 0.7818380260588590D+00, 0.4495382349988430D+00, & 0.4859612670925570D+00, 0.5251632870762940D+00, -0.6986022965844150D+00, & -0.9757586399688970D+00, 0.1384318633541960D+00, 0.1695042646494000D+00, & 0.3086023784018720D+00, -0.5931881528186310D+00, -0.7435673388472140D+00, & 0.9729796935790060D+00, -0.1911673832241180D+00, 0.1294818422565370D+00, & -0.6146246897809310D+00, 0.6821777798642300D+00, -0.3960681347586600D+00, & -0.6530289643965320D+00, -0.6449755112599790D+00, 0.3969379819746680D+00, & -0.7037890092249301D-01, 0.3208780019654030D+00, 0.9445020477265430D+00, & -0.3812529252505450D+00, 0.9096621317590370D+00, -0.1648059860305650D+00, & -0.3323417963042340D+00, -0.9834857390798001D-02, -0.9431077382830541D+00, & -0.4777466211688960D+00, -0.7551386761927890D+00, -0.4489139624465980D+00, & 0.3438775584320710D+00, 0.5740395992766760D+00, 0.7431196157208280D+00, & -0.8732125444955480D+00, 0.4700939413920300D+00, -0.1284972311068120D+00, & 0.6642168929664370D+00, 0.2599873469743290D+00, 0.7008726692568790D+00, & -0.8784891093226410D+00, -0.1706738463406710D+00, 0.4462368462787390D+00, & -0.3470827166082120D+00, 0.6266486359259690D+00, -0.6977428429758250D+00, & -0.4337167959777130D+00, -0.8857449345235880D+00, 0.1653652075033670D+00, & 0.6618616833629820D+00, 0.1125121287996140D+00, -0.7411343555448630D+00, & 0.4820689451276740D+00, 0.8658695321747411D+00, 0.1337141929452020D+00, & -0.8374660393934000D+00, -0.3724869462279710D+00, -0.3998801167256170D+00, & 0.4103552192662560D+00, -0.8216190506679300D+00, 0.3956649208616700D+00, & -0.3298995680158790D+00, 0.2926988290883000D-01, 0.9435621595726690D+00, & -0.9824290346165530D+00, -0.8096425490319800D-01, -0.1681605820944880D+00, & -0.9037042148768300D-01, -0.3161604362075780D+00, -0.9443917436621160D+00, & 0.5719599204934041D+00, -0.6863129715022710D+00, -0.4492620109656520D+00, & -0.4420214769968210D+00, 0.5021117496198080D+00, 0.7433039787107850D+00, & -0.7165157243440931D+00, -0.6847935061717610D+00, -0.1329024855773800D+00, & -0.4421804362881600D-01, 0.7098516256115680D+00, 0.7029619009834420D+00, & -0.1105565563628060D+00, -0.8896249757307140D+00, 0.4431079444123340D+00, & -0.7010281311842810D+00, -0.1342573855036490D+00, -0.7003816914554510D+00, & 0.7078411100140820D+00, -0.6867219567092810D+00, 0.1654506486763020D+00, & 0.9986011140880301D-01, 0.6665513378697570D+00, -0.7387403279457930D+00 & /), (/ 3, 70 /) ) return end subroutine design_12_084_3d ( x ) !*****************************************************************************80 ! !! DESIGN_12_084_3D returns a 12-design of 84 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,84), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,84) x(1:3,1:84) = reshape ( (/ & -0.8938049777611360D+00, -0.4268621911244970D+00, 0.1374821134468340D+00, & -0.4268621912410920D+00, 0.1374821134452880D+00, -0.8938049777056910D+00, & 0.8938049777701570D+00, 0.4268621911281570D+00, 0.1374821133768230D+00, & 0.1374821132964000D+00, -0.8938049777394910D+00, -0.4268621912182710D+00, & 0.4268621912727310D+00, 0.1374821133773450D+00, 0.8938049777010320D+00, & 0.1374821135290330D+00, 0.8938049777077750D+00, 0.4268621912097560D+00, & 0.4268621911859830D+00, -0.1374821134749930D+00, -0.8938049777274411D+00, & -0.1374821133242910D+00, 0.8938049777252790D+00, -0.4268621912390470D+00, & -0.4268621912174140D+00, -0.1374821133472880D+00, 0.8938049777320730D+00, & -0.1374821135010710D+00, -0.8938049776936550D+00, 0.4268621912483280D+00, & -0.8938049776725480D+00, 0.4268621913280710D+00, -0.1374821133907030D+00, & 0.8938049776635530D+00, -0.4268621913332600D+00, -0.1374821134330650D+00, & 0.9830866003855740D+00, 0.2230038010752200D-01, -0.1817785168533230D+00, & 0.2230038023239400D-01, -0.1817785168087260D+00, 0.9830866003909881D+00, & -0.9830866003966130D+00, -0.2230038011332300D-01, -0.1817785167929150D+00, & -0.1817785167104710D+00, 0.9830866004096310D+00, 0.2230038021145500D-01, & -0.2230038027285400D-01, -0.1817785168366860D+00, -0.9830866003849000D+00, & -0.1817785169360100D+00, -0.9830866003681790D+00, -0.2230038020037600D-01, & -0.2230038017080000D-01, 0.1817785168418750D+00, 0.9830866003862559D+00, & 0.1817785167109790D+00, -0.9830866004090441D+00, 0.2230038023321200D-01, & 0.2230038021255800D-01, 0.1817785168040810D+00, -0.9830866003922970D+00, & 0.1817785169343840D+00, 0.9830866003675029D+00, -0.2230038024343100D-01, & 0.9830866003916290D+00, -0.2230038033237200D-01, 0.1817785167929960D+00, & -0.9830866003805701D+00, 0.2230038033786500D-01, 0.1817785168521280D+00, & -0.8979519869718749D+00, 0.3766956030353650D+00, 0.2275580184196640D+00, & 0.3766956029275280D+00, 0.2275580183392060D+00, -0.8979519870375030D+00, & 0.8979519869860531D+00, -0.3766956030285690D+00, 0.2275580183749660D+00, & 0.2275580183055540D+00, -0.8979519870419040D+00, 0.3766956029373660D+00, & -0.3766956028752610D+00, 0.2275580184552540D+00, 0.8979519870300200D+00, & 0.2275580184865670D+00, 0.8979519869904800D+00, -0.3766956029506000D+00, & -0.3766956029825110D+00, -0.2275580183689100D+00, -0.8979519870069100D+00, & -0.2275580182809390D+00, 0.8979519870547670D+00, 0.3766956029215730D+00, & 0.3766956029314370D+00, -0.2275580184255800D+00, 0.8979519870139741D+00, & -0.2275580185113490D+00, -0.8979519870023480D+00, -0.3766956029073390D+00, & -0.8979519870721940D+00, -0.3766956028306370D+00, -0.2275580183627070D+00, & 0.8979519870578190D+00, 0.3766956028230510D+00, -0.2275580184319890D+00, & -0.1713301512452210D+00, 0.4597861949530550D+00, -0.8713453013615680D+00, & 0.4597861948431170D+00, -0.8713453014146491D+00, -0.1713301512702920D+00, & 0.1713301511912190D+00, -0.4597861949823340D+00, -0.8713453013567360D+00, & -0.8713453013647540D+00, -0.1713301511629810D+00, 0.4597861949776620D+00, & -0.4597861950424320D+00, -0.8713453013037380D+00, 0.1713301512994720D+00, & -0.8713453013534070D+00, 0.1713301513627270D+00, -0.4597861949247340D+00, & -0.4597861948552020D+00, 0.8713453014084100D+00, -0.1713301512695920D+00, & 0.8713453013928349D+00, 0.1713301511781830D+00, 0.4597861949187800D+00, & 0.4597861950544120D+00, 0.8713453013090380D+00, 0.1713301512403680D+00, & 0.8713453013254860D+00, -0.1713301513773550D+00, -0.4597861949721960D+00, & -0.1713301512966100D+00, -0.4597861949130030D+00, 0.8713453013725970D+00, & 0.1713301513507360D+00, 0.4597861949429830D+00, 0.8713453013461350D+00, & -0.3971917022972230D+00, -0.5480955906492260D+00, -0.7360910100912190D+00, & -0.5480955907789020D+00, -0.7360910100565570D+00, -0.3971917021825150D+00, & 0.3971917022502210D+00, 0.5480955906252050D+00, -0.7360910101344670D+00, & -0.7360910101747640D+00, -0.3971917021370830D+00, -0.5480955906530750D+00, & 0.5480955906102120D+00, -0.7360910101691310D+00, 0.3971917022066690D+00, & -0.7360910100491940D+00, 0.3971917023058890D+00, 0.5480955906993850D+00, & 0.5480955907525290D+00, 0.7360910100441170D+00, -0.3971917022419620D+00, & 0.7360910101399250D+00, 0.3971917021196020D+00, -0.5480955907125310D+00, & -0.5480955905843859D+00, 0.7360910101826250D+00, 0.3971917022173000D+00, & 0.7360910100837820D+00, -0.3971917022879800D+00, 0.5480955906659120D+00, & -0.3971917021252600D+00, 0.5480955907404190D+00, 0.7360910101161060D+00, & 0.3971917021713860D+00, -0.5480955907162950D+00, 0.7360910101091790D+00, & 0.3794747255349560D+00, 0.6962772780944900D+00, 0.6092592918368150D+00, & 0.6962772782104411D+00, 0.6092592917871140D+00, 0.3794747254020010D+00, & -0.3794747254955760D+00, -0.6962772780741610D+00, 0.6092592918845759D+00, & 0.6092592919259530D+00, 0.3794747253762130D+00, 0.6962772781030080D+00, & -0.6962772780710560D+00, 0.6092592919338879D+00, -0.3794747254221020D+00, & 0.6092592917959100D+00, -0.3794747255155420D+00, -0.6962772781408640D+00, & -0.6962772781859060D+00, -0.6092592917748491D+00, 0.3794747254667130D+00, & -0.6092592918828780D+00, -0.3794747253530890D+00, 0.6962772781533030D+00, & 0.6962772780465480D+00, -0.6092592919465890D+00, -0.3794747254466760D+00, & -0.6092592918387369D+00, 0.3794747254930950D+00, -0.6962772781156230D+00, & 0.3794747253362900D+00, -0.6962772781815950D+00, -0.6092592918610080D+00, & -0.3794747253752370D+00, 0.6962772781612160D+00, -0.6092592918600390D+00, & -0.6787014464703280D+00, 0.7297642134790810D+00, 0.8251387328409700D-01, & 0.7297642133897720D+00, 0.8251387317923400D-01, -0.6787014465791040D+00, & 0.6787014464747720D+00, -0.7297642134772200D+00, 0.8251387326399499D-01, & 0.8251387321767099D-01, -0.6787014465525470D+00, 0.7297642134101250D+00, & -0.7297642133709740D+00, 0.8251387336840200D-01, 0.6787014465763180D+00, & 0.8251387332689200D-01, 0.6787014465346920D+00, -0.7297642134143810D+00, & -0.7297642134312839D+00, -0.8251387320173600D-01, -0.6787014465317330D+00, & -0.8251387317169399D-01, 0.6787014465773989D+00, 0.7297642133922100D+00, & 0.7297642134127970D+00, -0.8251387334668001D-01, 0.6787014465339900D+00, & -0.8251387337365500D-01, -0.6787014465583360D+00, -0.7297642133871040D+00, & -0.6787014466415410D+00, -0.7297642133248270D+00, -0.8251387324006100D-01, & 0.6787014466370160D+00, 0.7297642133213440D+00, -0.8251387330807500D-01 & /), (/ 3, 84 /) ) return end subroutine design_13_094_3d ( x ) !*****************************************************************************80 ! !! DESIGN_13_094_3D returns a 13-design of 94 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,94), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,94) x(1:3,1:94) = reshape ( (/ & 0.2767901292869220D+00, -0.2352564669166030D+00, 0.9316875115097590D+00, & 0.1988867806345010D+00, 0.3605486031395280D+00, 0.9112896099830060D+00, & -0.2588713390623730D+00, 0.2042300774414090D+00, 0.9440739935409350D+00, & -0.2002829139273100D+00, -0.2283461619503540D+00, 0.9527564141538640D+00, & -0.8835451666675250D+00, -0.4142776966390410D+00, -0.2184534928214830D+00, & 0.3977500579085590D+00, -0.9016195359986890D+00, -0.1699326447132700D+00, & 0.8765394870692820D+00, 0.4343921041923270D+00, -0.2073210732744830D+00, & -0.4117423575176250D+00, 0.8848959788397900D+00, -0.2177781845341660D+00, & 0.5011140938672040D+00, 0.3778689327520590D+00, 0.7785240745079570D+00, & -0.3942388477903860D+00, 0.4736871338809520D+00, 0.7875253837741090D+00, & -0.4953642920021360D+00, -0.4064298087406120D+00, 0.7677428142133880D+00, & 0.3701865838021720D+00, -0.5593062709682520D+00, 0.7417131443007230D+00, & 0.4117423575179610D+00, -0.8848959788392530D+00, 0.2177781845357150D+00, & 0.8835451666683970D+00, 0.4142776966391570D+00, 0.2184534928177370D+00, & -0.3977500579105900D+00, 0.9016195359972180D+00, 0.1699326447163240D+00, & -0.8765394870698781D+00, -0.4343921041912780D+00, 0.2073210732741600D+00, & -0.6910143013156500D+00, -0.7028152268879870D+00, -0.1689674294993920D+00, & 0.6844003444601270D+00, -0.7144410442516540D+00, -0.1454997003140040D+00, & 0.6607104894827650D+00, 0.7313577150350630D+00, -0.1690489329931910D+00, & -0.7736112879563090D+00, 0.6152223578577780D+00, -0.1517465832844280D+00, & 0.6836297846860220D+00, -0.2199673313208400D+00, -0.6958912922588780D+00, & 0.2565740995035260D+00, 0.6814727910714180D+00, -0.6853937309994060D+00, & -0.6444745096378920D+00, 0.3540622279855340D+00, -0.6777112549905880D+00, & -0.2205350804161410D+00, -0.7315477541408590D+00, -0.6451373200469120D+00, & 0.3942388477920410D+00, -0.4736871338825220D+00, -0.7875253837723361D+00, & 0.4953642920009680D+00, 0.4064298087412850D+00, -0.7677428142137850D+00, & -0.3701865838024390D+00, 0.5593062709700030D+00, -0.7417131442992700D+00, & -0.5011140938646400D+00, -0.3778689327522390D+00, -0.7785240745095200D+00, & -0.4885748739685340D+00, -0.6884557978444000D-02, -0.8724948110952140D+00, & 0.5554204872744400D-01, -0.5841317202499910D+00, -0.8097562684048490D+00, & 0.5268121074647910D+00, 0.4970781961767100D-01, -0.8485270391079840D+00, & 0.4245864108125000D-02, 0.4886546223943000D+00, -0.8724669808369020D+00, & -0.7103173615146130D+00, -0.4795309146254010D+00, 0.5152662882912530D+00, & 0.5214043844765620D+00, -0.7280391654517230D+00, 0.4450802640164760D+00, & 0.7380993553888520D+00, 0.4078032732059310D+00, 0.5374996111041300D+00, & -0.4960579912625540D+00, 0.6996701137033650D+00, 0.5141869322482640D+00, & -0.9732208093073270D+00, 0.1942607517895710D+00, -0.1228983996858520D+00, & -0.3762035726666050D+00, -0.9088659640035350D+00, -0.1800931186603390D+00, & 0.9144779003707620D+00, -0.3686579885340490D+00, -0.1667976535311930D+00, & 0.2874621878541300D+00, 0.9468179143405530D+00, -0.1445729146068610D+00, & -0.9890066992933400D-01, 0.9950970592800400D+00, 0.7071773083110000D-03, & -0.9860682014252020D+00, -0.1613285617797790D+00, 0.4052896855503000D-01, & 0.9890066992737100D-01, -0.9950970592802360D+00, -0.7071773075220000D-03, & 0.9860682014253800D+00, 0.1613285617778720D+00, -0.4052896855829700D-01, & 0.8152324408482650D+00, 0.1318323811749280D+00, -0.5639293312661871D+00, & -0.1136445670803390D+00, 0.7872516056155811D+00, -0.6060691559787640D+00, & -0.7605044417053100D+00, -0.1056987489052100D-01, -0.6492466957107240D+00, & 0.1798482279122410D+00, -0.8324827854052400D+00, -0.5240486875479800D+00, & -0.9276898918095100D+00, -0.4724128948218800D-01, -0.3703508136922610D+00, & -0.6227375977374500D-01, -0.9444346856864090D+00, -0.3227461902425130D+00, & 0.9398402607408960D+00, -0.4456984180221600D-01, -0.3386942773242700D+00, & -0.8242738781550000D-02, 0.9340870594601500D+00, -0.3569501682398660D+00, & -0.2874621878541230D+00, -0.9468179143395800D+00, 0.1445729146132430D+00, & 0.9732208093070031D+00, -0.1942607517949340D+00, 0.1228983996799320D+00, & 0.3762035726640970D+00, 0.9088659640033310D+00, 0.1800931186666110D+00, & -0.9144779003729060D+00, 0.3686579885319570D+00, 0.1667976535240600D+00, & -0.1988867806309870D+00, -0.3605486031400650D+00, -0.9112896099835600D+00, & 0.2588713390641120D+00, -0.2042300774442540D+00, -0.9440739935398430D+00, & 0.2002829139245270D+00, 0.2283461619513520D+00, -0.9527564141542090D+00, & -0.2767901292844010D+00, 0.2352564669207660D+00, -0.9316875115094559D+00, & 0.4960579912585950D+00, -0.6996701137077700D+00, -0.5141869322460890D+00, & 0.7103173615128360D+00, 0.4795309146287610D+00, -0.5152662882905760D+00, & -0.5214043844766950D+00, 0.7280391654524530D+00, -0.4450802640151260D+00, & -0.7380993553844990D+00, -0.4078032732091310D+00, -0.5374996111076790D+00, & -0.8152324408494460D+00, -0.1318323811707120D+00, 0.5639293312654660D+00, & 0.1136445670801830D+00, -0.7872516056137170D+00, 0.6060691559812150D+00, & 0.7605044417099350D+00, 0.1056987488986400D-01, 0.6492466957053170D+00, & -0.1798482279168390D+00, 0.8324827854024680D+00, 0.5240486875508060D+00, & 0.6444745096387340D+00, -0.3540622279858040D+00, 0.6777112549896470D+00, & 0.2205350804135180D+00, 0.7315477541423701D+00, 0.6451373200460950D+00, & -0.6836297846853430D+00, 0.2199673313253120D+00, 0.6958912922581320D+00, & -0.2565740995009430D+00, -0.6814727910693790D+00, 0.6853937310024000D+00, & 0.8242738783470000D-02, -0.9340870594594580D+00, 0.3569501682416340D+00, & 0.9276898918126020D+00, 0.4724128947913300D-01, 0.3703508136849070D+00, & 0.6227375976878800D-01, 0.9444346856860160D+00, 0.3227461902446170D+00, & -0.9398402607419311D+00, 0.4456984180196000D-01, 0.3386942773214330D+00, & -0.6844003444607160D+00, 0.7144410442513049D+00, 0.1454997003129530D+00, & -0.6607104894826710D+00, -0.7313577150342460D+00, 0.1690489329970960D+00, & 0.7736112879557430D+00, -0.6152223578588770D+00, 0.1517465832828550D+00, & 0.6910143013134310D+00, 0.7028152268893190D+00, 0.1689674295029260D+00, & 0.8235860235780980D+00, -0.3946345889044380D+00, 0.4073936707989480D+00, & 0.4940686203583030D+00, 0.7088396084166290D+00, 0.5034308372726120D+00, & -0.7588751305010500D+00, 0.4506058872740210D+00, 0.4701732347348220D+00, & -0.4314990726013570D+00, -0.7879350487114470D+00, 0.4392799897061760D+00, & 0.4885748739746180D+00, 0.6884557979146000D-02, 0.8724948110918010D+00, & -0.5554204872563400D-01, 0.5841317202474440D+00, 0.8097562684068100D+00, & -0.5268121074675800D+00, -0.4970781961590400D-01, 0.8485270391063560D+00, & -0.4245864106237000D-02, -0.4886546223892350D+00, 0.8724669808397481D+00, & -0.4940686203577400D+00, -0.7088396084202140D+00, -0.5034308372681180D+00, & 0.7588751304965180D+00, -0.4506058872758780D+00, -0.4701732347403580D+00, & 0.4314990726012260D+00, 0.7879350487142150D+00, -0.4392799897013400D+00, & -0.8235860235777540D+00, 0.3946345889034440D+00, -0.4073936708006050D+00, & -0.5223814787874000D-01, -0.5618483004750600D-01, -0.9970528776242170D+00, & 0.5223814788153800D-01, 0.5618483004769000D-01, 0.9970528776240600D+00 & /), (/ 3, 94 /) ) return end subroutine design_14_108_3d ( x ) !*****************************************************************************80 ! !! DESIGN_14_108_3D returns a 14-design of 108 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,108), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,108) x(1:3,1:108) = reshape ( (/ & -0.6255209881602540D+00, -0.7673610045544000D+00, 0.1409985179364700D+00, & -0.7672427413700500D+00, 0.1411116382934610D+00, -0.6256405368525180D+00, & 0.6254929286339920D+00, 0.7673366024979470D+00, 0.1412555651851610D+00, & 0.1412599782857530D+00, -0.6254973365383090D+00, -0.7673321969774169D+00, & 0.7672171775977220D+00, 0.1411424450656330D+00, 0.6256649363676060D+00, & 0.1409941047324360D+00, 0.6255228978850370D+00, 0.7673602587130800D+00, & 0.7673671213519559D+00, -0.1410034708464950D+00, -0.6255123678051890D+00, & -0.1411074217386270D+00, 0.6256300078253110D+00, -0.7672521025313510D+00, & -0.7673415575795750D+00, -0.1412506125124700D+00, 0.6254879682905210D+00, & -0.1411466612794940D+00, -0.6256555691710420D+00, 0.7672240407957190D+00, & -0.6256319160421340D+00, 0.7672476983004790D+00, -0.1411229077154560D+00, & 0.6256599755682080D+00, -0.7672232962446800D+00, -0.1411311754058560D+00, & 0.5571880480715090D+00, -0.4475345647833600D-01, -0.8291794782913420D+00, & -0.4485587811454200D-01, -0.8292832145924940D+00, 0.5570254035443200D+00, & -0.5570231765605090D+00, 0.4489683751235000D-01, -0.8292824939402920D+00, & -0.8292769813537900D+00, 0.5570303319042050D+00, -0.4490988238058500D-01, & 0.4500608513203000D-01, -0.8291787590336450D+00, -0.5571687696457081D+00, & -0.8291849902738060D+00, -0.5571805246676550D+00, 0.4474499788474600D-01, & 0.4474511262700700D-01, 0.8291868870343290D+00, 0.5571776927213760D+00, & 0.8292859037266960D+00, -0.5570225725340150D+00, -0.4484131541096600D-01, & -0.4489531964450300D-01, 0.8292750865916650D+00, -0.5570343266194679D+00, & 0.8291760679008749D+00, 0.5571727652985150D+00, 0.4500619990677900D-01, & 0.5570150477881700D+00, 0.4485436035495000D-01, 0.8292902525019150D+00, & -0.5571799192991600D+00, -0.4499774138897400D-01, 0.8291717197297990D+00, & -0.2560655654109130D+00, 0.8607703824921130D+00, 0.4398917762759060D+00, & 0.8608171937494520D+00, 0.4399428934535150D+00, -0.2558202678543430D+00, & 0.2559781138440990D+00, -0.8608464350244180D+00, 0.4397938386773610D+00, & 0.4397804109275130D+00, -0.2559784640674000D+00, 0.8608531907927870D+00, & -0.8608968669378900D+00, 0.4397427222396980D+00, 0.2558963124660960D+00, & 0.4399052037052120D+00, 0.2560581296938110D+00, -0.8607657323399810D+00, & -0.8607663058143230D+00, -0.4398986385971350D+00, -0.2560674804326980D+00, & -0.4399515654145460D+00, 0.2558296190221410D+00, 0.8608099825863290D+00, & 0.8608459790026740D+00, -0.4397869770942070D+00, 0.2559914358201620D+00, & -0.4397340492183260D+00, -0.2559092846502780D+00, -0.8608974410391970D+00, & -0.2558221826891520D+00, -0.8608167398021460D+00, -0.4399426682200340D+00, & 0.2559096342559870D+00, 0.8608927923347450D+00, -0.4397429467325800D+00, & -0.2148474707463120D+00, -0.3239846898907800D-01, 0.9761100878082740D+00, & -0.3236168906853200D-01, 0.9761498728347380D+00, -0.2146721846102970D+00, & 0.2146533912587150D+00, 0.3222968714374900D-01, 0.9761583728513260D+00, & 0.9761568906841750D+00, -0.2146621640237140D+00, -0.3221614667304400D-01, & 0.3218487197615700D-01, 0.9761185894661391D+00, 0.2148409488773350D+00, & 0.9761115692642290D+00, 0.2148389643385310D+00, 0.3241024144420400D-01, & 0.3240438731759700D-01, -0.9761127401675890D+00, -0.2148345274044470D+00, & -0.9761504697174500D+00, 0.2146677480389070D+00, -0.3237311264470600D-01, & -0.3222757022524800D-01, -0.9761557221333460D+00, 0.2146657631381940D+00, & -0.9761179902309550D+00, -0.2148445483512370D+00, 0.3217901787241900D-01, & -0.2146592411125010D+00, 0.3235957222646100D-01, -0.9761527894189130D+00, & 0.2148533206001200D+00, -0.3219079038364400D-01, -0.9761156712406470D+00, & -0.5316579534180750D+00, -0.8273339530941490D+00, -0.1812687248945930D+00, & -0.8272321878124060D+00, -0.1811732915871120D+00, -0.5318487998130590D+00, & 0.5316939693677690D+00, 0.8273652744794220D+00, -0.1810199589093320D+00, & -0.1810139371845850D+00, -0.5316957717831260D+00, -0.8273654336584780D+00, & 0.8272650003242400D+00, -0.1811153925205650D+00, 0.5318174816895900D+00, & -0.1812747464882220D+00, 0.5316629623847160D+00, 0.8273294148729020D+00, & 0.8273379129643990D+00, 0.1812652358037510D+00, -0.5316529808631980D+00, & 0.1811784250571890D+00, 0.5318388191844240D+00, -0.8272374802330430D+00, & -0.8273707254762440D+00, 0.1810234483059680D+00, 0.5316842989860900D+00, & 0.1811102586154330D+00, -0.5318060097873780D+00, 0.8272734990184800D+00, & -0.5318438268704770D+00, 0.8272376408031750D+00, 0.1811629914142640D+00, & 0.5318078109206920D+00, -0.8272689621881080D+00, 0.1811256923905430D+00, & -0.6600529784314531D+00, -0.6410703014238900D+00, -0.3916106922647170D+00, & -0.6409432781620240D+00, -0.3914906263600130D+00, -0.6602475321053179D+00, & 0.6601308161987530D+00, 0.6411379932874460D+00, -0.3913685974476170D+00, & -0.3913661271946650D+00, -0.6601299904344370D+00, -0.6411403514158810D+00, & 0.6410141926516800D+00, -0.3914886640099250D+00, 0.6601798472922650D+00, & -0.3916131622327060D+00, 0.6600590826742870D+00, 0.6410625075180110D+00, & 0.6410745322055420D+00, 0.3916047718584040D+00, -0.6600523818722061D+00, & 0.3914935825762260D+00, 0.6602408324135060D+00, -0.6409483739081900D+00, & -0.6411454466952770D+00, 0.3913745185133930D+00, 0.6601200666850870D+00, & 0.3914857068510000D+00, -0.6601699246530760D+00, 0.6410262178062029D+00, & -0.6602469350652650D+00, 0.6409507331154460D+00, 0.3914794278831230D+00, & 0.6601690972978630D+00, -0.6410184249773530D+00, 0.3914998618294490D+00, & -0.8878095447864510D+00, -0.2962340013095760D+00, 0.3521926016460220D+00, & -0.2960667920239880D+00, 0.3523564025129960D+00, -0.8878003268014290D+00, & 0.8877394975983099D+00, 0.2961730845485540D+00, 0.3524203291424820D+00, & 0.3524163117776850D+00, -0.8877438356676250D+00, -0.2961648619048960D+00, & 0.2960029751816160D+00, 0.3522565288678170D+00, 0.8878612372176340D+00, & 0.3521966187543760D+00, 0.8878076464546130D+00, 0.2962349146112020D+00, & 0.2962433087631600D+00, -0.3522028941708700D+00, -0.8878023560177780D+00, & -0.3523574352146160D+00, 0.8877950378573241D+00, -0.2960814222555850D+00, & -0.2961794919204630D+00, -0.3524100372105150D+00, 0.8877414456017140D+00, & -0.3522554953175360D+00, -0.8878588486443870D+00, 0.2960113695493150D+00, & -0.8877931373857750D+00, 0.2960732001098630D+00, -0.3523691322852040D+00, & 0.8878631845739851D+00, -0.2960122833498500D+00, -0.3522437985034650D+00, & -0.2622350441333200D+00, -0.9631968323160830D+00, -0.5903087196255600D-01, & -0.9631467531778790D+00, -0.5898607663367000D-01, -0.2624289896453470D+00, & 0.2622467591823790D+00, 0.9632070203456430D+00, -0.5881218628148700D-01, & -0.5880201827075200D-01, -0.2622502022937960D+00, -0.9632067036956030D+00, & 0.9631574262320920D+00, -0.5885698170873300D-01, 0.2624188026763920D+00, & -0.5904103993067700D-01, 0.2622395302518100D+00, 0.9631949877243000D+00, & 0.9631989104617060D+00, 0.5903143538266000D-01, -0.2622272840915770D+00, & 0.5899352177531600D-01, 0.2624167439246010D+00, -0.9631496337000580D+00, & -0.9632095835159220D+00, 0.5881162296107100D-01, 0.2622374710596640D+00, & 0.5884953642617700D-01, -0.2624060718844200D+00, 0.9631613496712860D+00, & -0.2624212294120090D+00, 0.9631493186694710D+00, 0.5897871056935700D-01, & 0.2624095143628140D+00, -0.9631595066989900D+00, 0.5886434767523200D-01, & -0.7155075639675860D+00, -0.5512037701387860D+00, 0.4292124528598390D+00, & -0.5510696074929951D+00, 0.4293431805847270D+00, -0.7155324737444880D+00, & 0.7154222023624230D+00, 0.5511295351454290D+00, 0.4294500062373780D+00, & 0.4294508199172310D+00, -0.7154282714111380D+00, -0.5511210227691270D+00, & 0.5509918383737999D+00, 0.4293192792340250D+00, 0.7156067010051250D+00, & 0.4292116388668970D+00, 0.7155060331503000D+00, 0.5512063910974860D+00, & 0.5512118694302640D+00, -0.4292194620397150D+00, -0.7154971197744480D+00, & -0.4293412473175590D+00, 0.7155235615710740D+00, -0.5510826854369940D+00, & -0.5511341003106980D+00, -0.4294429977797900D+00, 0.7154228925136690D+00, & -0.4293212114664220D+00, -0.7156013233111630D+00, 0.5509973171080930D+00, & -0.7155220290294890D+00, 0.5510741739856940D+00, -0.4293547260011320D+00, & 0.7156073906348100D+00, -0.5509999389950980D+00, -0.4293077330962450D+00 & /), (/ 3, 108 /) ) return end subroutine design_15_120_3d ( x ) !*****************************************************************************80 ! !! DESIGN_15_120_3D returns a 15-design of 120 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,120), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,120) x = reshape ( (/ & .8544032798674690D+00, & -.5053541340072060D+00, & .1208810762424740D+00, & -.5054349175556900D+00, & .1208162198059960D+00, & .8543646675438200D+00, & -.8543867766655620D+00, & .5053247652039460D+00, & .1211202606115420D+00, & .1208333586366210D+00, & .8543977898340150D+00, & -.5053748273977881D+00, & .5053979097545750D+00, & .1211845075248970D+00, & -.8543344005432850D+00, & .1211678917817770D+00, & -.8543595921698920D+00, & .5053593070959080D+00, & .5055502439906060D+00, & -.1210290998832230D+00, & .8542662979320300D+00, & -.1209012300582570D+00, & -.8542712308899000D+00, & -.5055725038451520D+00, & -.5055127430804750D+00, & -.1209718018932920D+00, & -.8542966052431350D+00, & -.1211002217099370D+00, & .8542333106705450D+00, & .5055889508708080D+00, & .8542280860180770D+00, & .5056228637404400D+00, & -.1209954409091880D+00, & -.8542442519150010D+00, & -.5055933391780419D+00, & -.1210046835828190D+00, & -.2649878987783750D+00, & .8838136985753620D+00, & -.3855577255244170D+00, & .8838495436614180D+00, & -.3855147723237870D+00, & -.2649308296322680D+00, & .2649355775867900D+00, & -.8837175312641120D+00, & -.3858140286008490D+00, & -.3855859698282140D+00, & -.2649918821321200D+00, & .8838001823238730D+00, & -.8837292209025740D+00, & -.3858572041282210D+00, & .2648336877084960D+00, & -.3857855708943260D+00, & .2648743412798280D+00, & -.8837483106752260D+00, & -.8838853342582800D+00, & .3855799306321350D+00, & -.2647165143636620D+00, & .3857055651888330D+00, & .2647713629504000D+00, & .8838140881111540D+00, & .8837647796773470D+00, & .3857919932898060D+00, & .2648100302392800D+00, & .3856671594356380D+00, & -.2646531289490930D+00, & -.8838662588142510D+00, & -.2646569340063700D+00, & -.8838970336929440D+00, & .3855940107304080D+00, & .2647089806257060D+00, & .8838007658213300D+00, & .3857789028831530D+00, & -.9733521648930310D+00, & .2280262392539510D+00, & .2428162494047400D-01, & .2281123458869260D+00, & .2434454914745700D-01, & -.9733304169606380D+00, & .9733554428092021D+00, & -.2280321804631890D+00, & .2409370513031500D-01, & .2435678809852000D-01, & -.9733428333593150D+00, & .2280580531829240D+00, & -.2281196464116130D+00, & .2403159282358800D-01, & .9733364831687970D+00, & .2401983627102700D-01, & .9733503190618910D+00, & -.2280618421560860D+00, & -.2282433008441410D+00, & -.2410134029444700D-01, & -.9733057695379100D+00, & -.2432552782893800D-01, & .9732858593312290D+00, & .2283044124008880D+00, & .2282511002315090D+00, & -.2427391182448100D-01, & .9732996519304030D+00, & -.2405063099923100D-01, & -.9732933319509690D+00, & -.2283016800821220D+00, & -.9732840684188760D+00, & -.2283305892567410D+00, & -.2415086275208100D-01, & .9732807668974420D+00, & .2283367922856280D+00, & -.2422515379181500D-01, & .1764941643605970D+00, & -.6439159617576869D+00, & .7444609084030720D+00, & -.6439554747347650D+00, & .7444298798961950D+00, & .1764808785020650D+00, & -.1763945457957580D+00, & .6437302463633759D+00, & .7446451061616230D+00, & .7444821769849420D+00, & .1765407684656810D+00, & -.6438785950948430D+00, & .6437221735449969D+00, & .7446755654083020D+00, & -.1762953929356370D+00, & .7446232602985000D+00, & -.1763087546453890D+00, & .6437790174103420D+00, & .6439796469919020D+00, & -.7444731281720970D+00, & .1762100328864330D+00, & -.7445688151995640D+00, & -.1763819450908200D+00, & -.6438219387985280D+00, & -.6437453004186960D+00, & -.7446323171490740D+00, & -.1763935952523330D+00, & -.7445366073086630D+00, & .1761491742553250D+00, & .6439229059339890D+00, & .1761957070235800D+00, & .6439711185624940D+00, & -.7444838959197370D+00, & -.1762957903461830D+00, & -.6437847990685870D+00, & -.7446213311438470D+00, & -.5382680916693570D+00, & -.7144430972079280D+00, & .4470330215348440D+00, & -.7144085713422340D+00, & .4470401217976180D+00, & -.5383080184206051D+00, & .5383282112194380D+00, & .7143316685506840D+00, & .4471386857686070D+00, & .4471201978095990D+00, & -.5382526202846060D+00, & -.7144001997952280D+00, & .7142693383899240D+00, & .4471311678997200D+00, & .5384171532637620D+00, & .4470507600303290D+00, & .5383906456527831D+00, & .7143396465476940D+00, & .7143353780369111D+00, & -.4469080442623900D+00, & -.5385147794243240D+00, & -.4472169730606700D+00, & .5384307414350790D+00, & -.7142053736035060D+00, & -.7141970542020020D+00, & -.4472644791667070D+00, & .5384023156013700D+00, & -.4469551475186520D+00, & -.5385686162875321D+00, & .7142653160112939D+00, & -.5385531358716100D+00, & .7142734853459000D+00, & -.4469607454516280D+00, & .5384924234925470D+00, & -.7141627498790130D+00, & -.4472107741617700D+00, & -.8542625591715191D+00, & -.1211967864813340D+00, & .5055163884033080D+00, & -.1211352259177740D+00, & .5055628408198980D+00, & -.8542438006939030D+00, & .8543309414267219D+00, & .1210719374548010D+00, & .5054307355927919D+00, & .5056009275709700D+00, & -.8542205006653180D+00, & -.1211405724044100D+00, & .1209804838657270D+00, & .5053853800887540D+00, & .8543707275627840D+00, & .5053474173253790D+00, & .8543800060355530D+00, & .1210735028371490D+00, & .1210201874725290D+00, & -.5053486406391460D+00, & -.8543868360574620D+00, & -.5056166845687340D+00, & .8542509316027460D+00, & -.1208598946779390D+00, & -.1208657310279020D+00, & -.5055984823689040D+00, & .8542608791752970D+00, & -.5053307603087600D+00, & -.8544103980209300D+00, & .1209284682756170D+00, & -.8543680863966890D+00, & .1209284250962590D+00, & -.5054023040614261D+00, & .8543001861498960D+00, & -.1208044489183890D+00, & -.5055467110621700D+00, & .7444632973046910D+00, & .6439458791655150D+00, & .1763748952430030D+00, & .6438740925860960D+00, & .1763526546392770D+00, & .7445306535651250D+00, & -.7444394383517560D+00, & -.6439896739556450D+00, & .1763156897868830D+00, & .1762725384812350D+00, & .7444994067041240D+00, & .6439321591554410D+00, & -.6439272279941870D+00, & .1763373486135990D+00, & -.7444883239726789D+00, & .1764171560921340D+00, & -.7444453593527010D+00, & -.6439550403363510D+00, & -.6437732764367420D+00, & -.1765375814946370D+00, & .7445740063055700D+00, & -.1761862890755310D+00, & -.7446599355387570D+00, & .6437701235264091D+00, & .6438276103095120D+00, & -.1761530561791470D+00, & -.7446180960746860D+00, & -.1765038577432470D+00, & .7446064263886381D+00, & -.6437450255982510D+00, & .7446425637975280D+00, & -.6437113057563880D+00, & -.1764743806405070D+00, & -.7446660895457380D+00, & .6437550517808310D+00, & -.1762153466282650D+00, & -.2283360715319860D+00, & .9732852424848990D+00, & -.2405151135480500D-01, & .9733064515717200D+00, & -.2400727784605000D-01, & -.2282503054527860D+00, & .2283331981063670D+00, & -.9732791056699990D+00, & -.2432556492099900D-01, & -.2408702223521400D-01, & -.2283063346170040D+00, & .9732913402129850D+00, & -.9732989361424010D+00, & -.2436839438764100D-01, & .2282440848281470D+00, & -.2428855197935900D-01, & .2282995532397890D+00, & -.9732879225764900D+00, & -.9733371742315819D+00, & .2406421970882400D-01, & -.2281132582483620D+00, & .2421829905214100D-01, & .2280639631571710D+00, & .9733449042862790D+00, & .9733297074180570D+00, & .2431156319331100D-01, & .2281188912663230D+00, & .2415744670603300D-01, & -.2280561009745210D+00, & -.9733482586294300D+00, & -.2280265869033880D+00, & -.9733570472797149D+00, & .2408182252460000D-01, & .2280302830467680D+00, & .9733509134612350D+00, & .2429381151218600D-01, & .7141882216414780D+00, & .5385777705783590D+00, & -.4470672981861090D+00, & .5385274786436980D+00, & -.4470918044321980D+00, & .7142108044234720D+00, & -.7142487687796120D+00, & -.5384672206579190D+00, & -.4471037335716890D+00, & -.4471613350315650D+00, & .7141657502813270D+00, & .5385294992643370D+00, & -.5383895843123190D+00, & -.4470804117879510D+00, & -.7143218888565050D+00, & -.4470109694572090D+00, & -.7143051130884170D+00, & -.5384694964440140D+00, & -.5384320154150550D+00, & .4469136813525910D+00, & .7143942372359660D+00, & .4472327482991630D+00, & -.7143004980317570D+00, & .5382914334822370D+00, & .5382930399261510D+00, & .4472578784152030D+00, & -.7142835524934020D+00, & .4469384509810380D+00, & .7144394343055710D+00, & -.5383514797451630D+00, & .7144166137892710D+00, & -.5383541595409150D+00, & .4469717002751600D+00, & -.7143571305190370D+00, & .5382432243381430D+00, & .4472003147703360D+00, & -.8838226247963700D+00, & .2646362039654440D+00, & .3857787545327180D+00, & .2647035776402780D+00, & .3858370666585500D+00, & -.8837769933511300D+00, & .8838743905237270D+00, & -.2647316811045790D+00, & .3855946042099840D+00, & .3858426620308480D+00, & -.8837798638376970D+00, & .2646858372334470D+00, & -.2648227787783800D+00, & .3855352957965560D+00, & .8838729725108470D+00, & .3855302957630180D+00, & .8839005568372480D+00, & -.2647379773883650D+00, & -.2648230436494790D+00, & -.3856162950490920D+00, & -.8838375577813140D+00, & -.3858066593719060D+00, & .8837138197040250D+00, & .2649586881919760D+00, & .2649414017790060D+00, & -.3857554525672310D+00, & .8837413560754220D+00, & -.3855652667157560D+00, & -.8838346669397210D+00, & -.2649069772919530D+00, & -.8837919506738630D+00, & -.2649137504830990D+00, & -.3856585182892600D+00, & .8837398529718770D+00, & .2650088880255400D+00, & -.3857125374378070D+00 /), (/ 3, 120 /) ) return end subroutine design_16_144_3d ( x ) !*****************************************************************************80 ! !! DESIGN_16_144_3D returns a 16-design of 144 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,144), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,144) x = reshape ( (/ & 0.9383118258138560D+00, -0.1750792557749200D+00, -0.2981915017822760D+00, & -0.1751096322456290D+00, -0.2982825311210240D+00, 0.9382772235980340D+00, & -0.9383116523013461D+00, 0.1751477614500080D+00, -0.2981518150449020D+00, & -0.2981827578157150D+00, 0.9383270575537280D+00, -0.1750125024219040D+00, & 0.1750977124101310D+00, -0.2980583478457380D+00, -0.9383506873169580D+00, & -0.2981854777577620D+00, -0.9383236127415390D+00, 0.1750263369497320D+00, & 0.1751212256614090D+00, 0.2980709997422250D+00, 0.9383422805328110D+00, & 0.2981590222823750D+00, -0.9382974848874340D+00, -0.1752113788700180D+00, & -0.1751366381351110D+00, 0.2982885002265250D+00, -0.9382702854803310D+00, & 0.2981750565054620D+00, 0.9382926280748330D+00, 0.1752101018160420D+00, & 0.9383097216767580D+00, 0.1750911370548140D+00, 0.2981911466354040D+00, & -0.9383070207140820D+00, -0.1751442959881740D+00, 0.2981684263322820D+00, & 0.3183193898656830D+00, -0.1895522954118680D+00, 0.9288394335619220D+00, & -0.1894661062614570D+00, 0.9288339463361680D+00, 0.3183867062421130D+00, & -0.3182933144730710D+00, 0.1893628596173800D+00, 0.9288870078536330D+00, & 0.9288529435535660D+00, 0.3183507003489590D+00, -0.1894334733863170D+00, & 0.1894416073975330D+00, 0.9288927988957520D+00, -0.3182295485121640D+00, & 0.9288662644063450D+00, -0.3183138373071290D+00, 0.1894301027466670D+00, & 0.1894518259149400D+00, -0.9288871565521020D+00, 0.3182399347191460D+00, & -0.9288657503320540D+00, -0.3182891226867960D+00, -0.1894741466251780D+00, & -0.1894810419822530D+00, -0.9288341329001750D+00, -0.3183772735119440D+00, & -0.9288638749080860D+00, 0.3182773954415380D+00, 0.1895030380803610D+00, & 0.3182754841245910D+00, 0.1895728238082200D+00, -0.9288502897015400D+00, & -0.3183459025831120D+00, -0.1893534180173150D+00, -0.9288709110493790D+00, & 0.4152709071162880D+00, 0.6265468605244530D+00, 0.6595370385882560D+00, & 0.6266126549472570D+00, 0.6594514158910070D+00, 0.4153076097777360D+00, & -0.4152418281129630D+00, -0.6266763943801670D+00, 0.6594322716641020D+00, & 0.6594942179223080D+00, 0.4151968471621200D+00, 0.6266410093775210D+00, & -0.6266189964270690D+00, 0.6595218123324770D+00, -0.4151862381804330D+00, & 0.6594787856877939D+00, -0.4151922150229020D+00, -0.6266603193215040D+00, & -0.6266022331854350D+00, -0.6595287758101400D+00, 0.4152004759696260D+00, & -0.6594726936833410D+00, -0.4153261780732930D+00, 0.6265779537240910D+00, & 0.6266060528732360D+00, -0.6594438365947900D+00, -0.4153296051087230D+00, & -0.6594986338231030D+00, 0.4153157815166040D+00, -0.6265575421369630D+00, & 0.4152509631584860D+00, -0.6265428543902710D+00, -0.6595534013318720D+00, & -0.4152672330732850D+00, 0.6266741585574390D+00, -0.6594183983875370D+00, & 0.8147686975402800D-01, 0.8847674930322230D+00, 0.4588551001880220D+00, & 0.8848021501705901D+00, 0.4587806295976860D+00, 0.8151986849505800D-01, & -0.8148097265168000D-01, -0.8848439651039500D+00, 0.4587068873636580D+00, & 0.4587780511560210D+00, 0.8139667888042000D-01, 0.8848148283368230D+00, & -0.8848095158928860D+00, 0.4587845170278200D+00, -0.8141798032957800D-01, & 0.4587323275728680D+00, -0.8138617295209800D-01, -0.8848395009784490D+00, & -0.8848064690255750D+00, -0.4587847998886890D+00, 0.8144949208920500D-01, & -0.4587707681467430D+00, -0.8156762447812400D-01, 0.8848028621851550D+00, & 0.8848211768135870D+00, -0.4587411019232240D+00, -0.8153579869288199D-01, & -0.4588108991167440D+00, 0.8157388735636099D-01, -0.8847814757064350D+00, & 0.8147060004176100D-01, -0.8847779034947540D+00, -0.4588361393964780D+00, & -0.8149754501781800D-01, 0.8848501095969900D+00, -0.4586920902983440D+00, & -0.7225816121467720D+00, 0.6911694469079300D+00, -0.1267317830534700D-01, & 0.6911462318877840D+00, -0.1272247709073500D-01, -0.7226029509516230D+00, & 0.7225897391740940D+00, -0.6911572322235680D+00, -0.1287436155202900D-01, & -0.1271999109003300D-01, -0.7226498291394290D+00, 0.6910972625263569D+00, & -0.6911640879369000D+00, -0.1283280970189800D-01, 0.7225839207604250D+00, & -0.1274089462228200D-01, 0.7226581266795230D+00, -0.6910882009904870D+00, & -0.6911848254516650D+00, 0.1280693240541800D-01, -0.7225645435168510D+00, & 0.1278690708865000D-01, 0.7225094351193579D+00, 0.6912428018942500D+00, & 0.6911667589030220D+00, 0.1267926954320300D-01, 0.7225840764307940D+00, & 0.1279840273451600D-01, -0.7225177746588930D+00, -0.6912338722815931D+00, & -0.7225871989731150D+00, -0.6911634956048890D+00, 0.1267920436853000D-01, & 0.7225783528006580D+00, 0.6911703359443890D+00, 0.1280979212978900D-01, & 0.5601175739954590D+00, 0.8068680228904130D+00, 0.1877026822886580D+00, & 0.8068834787163790D+00, 0.1875714426539700D+00, 0.5601392734626480D+00, & -0.5601340935408990D+00, -0.8068916312063850D+00, 0.1875518401461700D+00, & 0.1876521312373620D+00, 0.5600251494167630D+00, 0.8069439321680340D+00, & -0.8068854415129990D+00, 0.1876857418815800D+00, -0.5600981579765580D+00, & 0.1876302229010670D+00, -0.5600047208391949D+00, -0.8069632036790220D+00, & -0.8068746775941580D+00, -0.1876975169586680D+00, 0.5601097185238560D+00, & -0.1876148088020380D+00, -0.5602157603217920D+00, 0.8068202931293010D+00, & 0.8068923207022480D+00, -0.1875761333114300D+00, -0.5601249655243670D+00, & -0.1876364876816170D+00, 0.5602249742367100D+00, -0.8068088539003420D+00, & 0.5600918210887200D+00, -0.8068808651992270D+00, -0.1877243226778800D+00, & -0.5601293840974760D+00, 0.8068962450831860D+00, -0.1875460549871360D+00, & -0.9948563422103200D-01, -0.3588951295179950D+00, -0.9280608248341810D+00, & -0.3590507942888110D+00, -0.9279946080877720D+00, -0.9954162185034500D-01, & 0.9943438966061501D-01, 0.3591437619459990D+00, -0.9279701290494740D+00, & -0.9280190267200989D+00, -0.9942019096838001D-01, -0.3590213248169130D+00, & 0.3589908155319930D+00, -0.9280357484444770D+00, 0.9937426212442400D-01, & -0.9280072072034909D+00, 0.9942025966856401D-01, 0.3590518560679110D+00, & 0.3590029825622480D+00, 0.9280313484672880D+00, -0.9937139816565700D-01, & 0.9280179389220590D+00, 0.9951094937970199D-01, -0.3589989916314580D+00, & -0.3590428637423850D+00, 0.9279961920762100D+00, 0.9955545935668900D-01, & 0.9280136656320840D+00, -0.9948954910509600D-01, 0.3590159690305810D+00, & -0.9945187531254500D-01, 0.3589267513480540D+00, 0.9280522138670591D+00, & 0.9946550331739699D-01, -0.3591200632919870D+00, 0.9279759661709870D+00, & 0.7878331994376070D+00, 0.5574500823251660D+00, -0.2618554096816970D+00, & 0.5574053886878521D+00, -0.2619772920486170D+00, 0.7878243021845780D+00, & -0.7878614778767180D+00, -0.5573641116878389D+00, -0.2619533139327300D+00, & -0.2618610280706080D+00, 0.7878026576023160D+00, 0.5574906049903740D+00, & -0.5574272044780030D+00, -0.2618353048552930D+00, -0.7878560686059189D+00, & -0.2618500918686550D+00, -0.7878041469245110D+00, -0.5574936371627220D+00, & -0.5573980474810630D+00, 0.2618066241900950D+00, 0.7878862279507650D+00, & 0.2619281468060600D+00, -0.7878933745001880D+00, 0.5573308499710470D+00, & 0.5573998343635920D+00, 0.2619357785378840D+00, -0.7878420352920970D+00, & 0.2619095353283640D+00, 0.7879083713373940D+00, -0.5573183952468600D+00, & 0.7878589677335660D+00, -0.5574443214494930D+00, 0.2617901362647470D+00, & -0.7878560239272930D+00, 0.5573693294883240D+00, 0.2619586152567080D+00, & -0.5072827321686140D+00, -0.7170499460473529D+00, -0.4780205063771150D+00, & -0.7170643140017600D+00, -0.4779062710060660D+00, -0.5073700481091310D+00, & 0.5073317531927670D+00, 0.7171162628030801D+00, -0.4778689755839950D+00, & -0.4778916169164080D+00, -0.5072575026701600D+00, -0.7171536993321960D+00, & 0.7171087443614590D+00, -0.4779885176598600D+00, 0.5072297563685140D+00, & -0.4779136769269750D+00, 0.5072353404128420D+00, 0.7171546742805260D+00, & 0.7171036377589221D+00, 0.4779420929439370D+00, -0.5072807196270021D+00, & 0.4779497913306490D+00, 0.5073623117813870D+00, -0.7170407809538000D+00, & -0.7170736052366210D+00, 0.4778899243873540D+00, 0.5073723138307850D+00, & 0.4779668855044820D+00, -0.5073968954290570D+00, 0.7170049141185160D+00, & -0.5072894944901550D+00, 0.7170398743210130D+00, 0.4780284378712520D+00, & 0.5073429733358930D+00, -0.7171476166924809D+00, 0.4778100075123900D+00, & -0.4697053900856580D+00, -0.3362487640635100D+00, 0.8162803533040850D+00, & -0.3361804588591880D+00, 0.8163540175197370D+00, -0.4696262526314000D+00, & 0.4697292672795090D+00, 0.3360874275716510D+00, 0.8163330548797630D+00, & 0.8162993201022140D+00, -0.4696884809882010D+00, -0.3362263386881820D+00, & 0.3361661885920780D+00, 0.8162610446461910D+00, 0.4697980423975660D+00, & 0.8163081878411480D+00, 0.4696844879905110D+00, 0.3362103868184210D+00, & 0.3361611964247630D+00, -0.8162545204851160D+00, -0.4698129498065010D+00, & -0.8163147475476900D+00, 0.4697491969062010D+00, -0.3361040388661380D+00, & -0.3361667115393140D+00, -0.8163550823770680D+00, 0.4696342422885870D+00, & -0.8163020291364350D+00, -0.4697523383167870D+00, 0.3361305369549900D+00, & -0.4697259147648690D+00, 0.3362543092749910D+00, -0.8162662583326020D+00, & 0.4697157090205860D+00, -0.3360825711370180D+00, -0.8163428557151829D+00, & 0.2209757831175440D+00, 0.5619818996413199D+00, -0.7970859726222270D+00, & 0.5618985433861100D+00, -0.7971884426164270D+00, 0.2208180560990520D+00, & -0.2209098087123600D+00, -0.5618199356383180D+00, -0.7972184244822900D+00, & -0.7971433029261999D+00, 0.2209065606243460D+00, 0.5619277943588750D+00, & -0.5619110464580350D+00, -0.7971135607042630D+00, -0.2210564344456110D+00, & -0.7971666081668140D+00, -0.2209214541641100D+00, -0.5618888768376120D+00, & -0.5619031892145560D+00, 0.7971171958991410D+00, 0.2210632985196790D+00, & 0.7971490712061960D+00, -0.2210199177081820D+00, 0.5618750343727400D+00, & 0.5618715451873800D+00, 0.7971902229922720D+00, -0.2208803184402730D+00, & 0.7971513117794931D+00, 0.2209665013294830D+00, -0.5618928647157230D+00, & 0.2209896742277390D+00, -0.5619592284389200D+00, 0.7970981052900900D+00, & -0.2209345147362070D+00, 0.5618214796441769D+00, 0.7972104899013210D+00, & -0.2558632109166300D-01, 0.9914006599926770D+00, -0.1283357765359230D+00, & 0.9913910231541920D+00, -0.1284105096544480D+00, -0.2558476538037500D-01, & 0.2555531861480000D-01, -0.9913788670530650D+00, -0.1285101850091180D+00, & -0.1284273557345780D+00, -0.2568716764003100D-01, 0.9913861930235141D+00, & -0.9913884123582900D+00, -0.1284322897283740D+00, 0.2557660643704000D-01, & -0.1284710461501210D+00, 0.2569665752758400D-01, -0.9913802863144920D+00, & -0.9913887704923130D+00, 0.1284336117570290D+00, -0.2555607780520200D-01, & 0.1284346430688090D+00, 0.2546732329508000D-01, 0.9913909208299070D+00, & 0.9913860541495390D+00, 0.1284483459343360D+00, 0.2558738096289900D-01, & 0.1283929891583590D+00, -0.2544975216483000D-01, -0.9913967674194480D+00, & -0.2558970505166500D-01, -0.9913982207318930D+00, 0.1283539439402070D+00, & 0.2557174693595500D-01, 0.9913764768665120D+00, 0.1285253549864190D+00 & /), (/ 3, 144 /) ) return end subroutine design_17_156_3d ( x ) !*****************************************************************************80 ! !! DESIGN_17_156_3D returns a 17-design of 156 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,156), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,156) x = reshape ( (/ & -0.5389531643378300D-01, -0.1406035066764100D+00, -0.9885979712586910D+00, & -0.1406020108260560D+00, -0.9885983027651530D+00, -0.5389313798182900D-01, & 0.5389273741214000D-01, 0.1406029924863770D+00, -0.9885981849862470D+00, & -0.9885980986472160D+00, -0.5389529965961800D-01, -0.1406026174212630D+00, & 0.1406044785163560D+00, -0.9885978843519180D+00, 0.5389437518011600D-01, & -0.9885980658632700D+00, 0.5389281342061700D-01, 0.1406038009195060D+00, & 0.1406018124157300D+00, 0.9885982766199050D+00, -0.5389413520565200D-01, & 0.9885980025098750D+00, 0.5389565255163500D-01, -0.1406031581064890D+00, & -0.1406048466356920D+00, 0.9885978500353531D+00, 0.5389404427235700D-01, & 0.9885981673609280D+00, -0.5389245072780000D-01, 0.1406032262971460D+00, & -0.5389283591088400D-01, 0.1406027515842190D+00, 0.9885982138788380D+00, & 0.5389609745044300D-01, -0.1406037938522460D+00, 0.9885978878360840D+00, & -0.7121378206194821D+00, 0.4847259556271390D+00, 0.5078390221169400D+00, & 0.4847271423032010D+00, 0.5078385891579620D+00, -0.7121373216428299D+00, & 0.7121378784277490D+00, -0.4847268954121660D+00, 0.5078380440381630D+00, & 0.5078397607384350D+00, -0.7121379693767980D+00, 0.4847249632369050D+00, & -0.4847276428704660D+00, 0.5078383908533870D+00, 0.7121371223385879D+00, & 0.5078396078145549D+00, 0.7121361919763640D+00, -0.4847277347255500D+00, & -0.4847267680672810D+00, -0.5078392368723810D+00, -0.7121371144744010D+00, & -0.5078401126327480D+00, 0.7121362579125310D+00, 0.4847271089669900D+00, & 0.4847285491224600D+00, -0.5078400954114271D+00, 0.7121352899261120D+00, & -0.5078388610154030D+00, -0.7121374661713530D+00, -0.4847266451492260D+00, & -0.7121366718689039D+00, -0.4847276329363270D+00, -0.5078390320243490D+00, & 0.7121377658573640D+00, 0.4847264015559710D+00, -0.5078386732755610D+00, & -0.7030054480395250D+00, 0.2617901117095170D+00, 0.6612482721624801D+00, & 0.2617908536144600D+00, 0.6612471360366760D+00, -0.7030062404041000D+00, & 0.7030064339445450D+00, -0.2617905690855730D+00, 0.6612470429199860D+00, & 0.6612494875894130D+00, -0.7030044337858500D+00, 0.2617897653464990D+00, & -0.2617917110517330D+00, 0.6612472320733990D+00, 0.7030058307723160D+00, & 0.6612474233592150D+00, 0.7030059089592190D+00, -0.2617910179307560D+00, & -0.2617910421351510D+00, -0.6612484386100850D+00, -0.7030049449993340D+00, & -0.6612490446749040D+00, 0.7030044247634020D+00, 0.2617909083211350D+00, & 0.2617922462673400D+00, -0.6612489166629980D+00, 0.7030040469345190D+00, & -0.6612472545309420D+00, -0.7030062605254830D+00, -0.2617905002807940D+00, & -0.7030054694842800D+00, -0.2617913268030810D+00, -0.6612477683031300D+00, & 0.7030051584635270D+00, 0.2617910253582180D+00, -0.6612482183080450D+00, & 0.6280044724638099D-01, 0.7862188199982439D+00, -0.6147487868277770D+00, & 0.7862200431089770D+00, -0.6147474493883091D+00, 0.6279822676069299D-01, & -0.6279950225219801D-01, -0.7862192395650210D+00, -0.6147483467685591D+00, & -0.6147477070961400D+00, 0.6279957151438099D-01, 0.7862197341949950D+00, & -0.7862185345191240D+00, -0.6147492348350890D+00, -0.6279963573361200D-01, & -0.6147493306965400D+00, -0.6279995662861700D-01, -0.7862184339327080D+00, & -0.7862195385712860D+00, 0.6147479435281090D+00, 0.6279970618335100D-01, & 0.6147493951504540D+00, -0.6279977045814100D-01, 0.7862183984062920D+00, & 0.7862177984580020D+00, 0.6147501790519670D+00, -0.6279960782860100D-01, & 0.6147478002140190D+00, 0.6280005812960500D-01, -0.7862196225171070D+00, & 0.6280052636345899D-01, -0.7862189094158020D+00, 0.6147486643869180D+00, & -0.6280141239775699D-01, 0.7862187128109530D+00, 0.6147488253154580D+00, & 0.8295436077392320D+00, 0.3214653682205850D+00, 0.4566370767839410D+00, & 0.3214635955020470D+00, 0.4566373806324790D+00, 0.8295441274435040D+00, & -0.8295434450385299D+00, -0.3214644597436460D+00, 0.4566380119036660D+00, & 0.4566350005565370D+00, 0.8295451597747751D+00, 0.3214643124220390D+00, & -0.3214642002577900D+00, 0.4566374595738670D+00, -0.8295438496345710D+00, & 0.4566370685585350D+00, -0.8295440001148970D+00, -0.3214643673747430D+00, & -0.3214629541954330D+00, -0.4566364211023370D+00, 0.8295449041509410D+00, & -0.4566367130348990D+00, -0.8295441909717370D+00, 0.3214643798832610D+00, & 0.3214629551065890D+00, -0.4566366887995170D+00, -0.8295447564395499D+00, & -0.4566371123963230D+00, 0.8295440985782710D+00, -0.3214640510170780D+00, & 0.8295447019765780D+00, -0.3214622985067580D+00, -0.4566372499712900D+00, & -0.8295448614467950D+00, 0.3214635893905120D+00, -0.4566360515151940D+00, & -0.2495004234484620D+00, 0.9540250943623850D+00, -0.1660893073797370D+00, & 0.9540254708554060D+00, -0.1660875670107380D+00, -0.2495001423718530D+00, & 0.2495009434844220D+00, -0.9540250296126640D+00, -0.1660888981026120D+00, & -0.1660868774086620D+00, -0.2495001374496830D+00, 0.9540255921961580D+00, & -0.9540248553834940D+00, -0.1660899370031220D+00, 0.2495009181081350D+00, & -0.1660901519985670D+00, 0.2495001071183790D+00, -0.9540250300474360D+00, & -0.9540255936888939D+00, 0.1660878626585790D+00, -0.2494994758793280D+00, & 0.1660896928744990D+00, 0.2494996878225310D+00, 0.9540252196337971D+00, & 0.9540249314198170D+00, 0.1660906479136480D+00, 0.2495001541182640D+00, & 0.1660879561220760D+00, -0.2495000023520480D+00, -0.9540254397328820D+00, & -0.2494997599825380D+00, -0.9540252259304089D+00, 0.1660895483077950D+00, & 0.2494983747081790D+00, 0.9540257202571130D+00, 0.1660887898260250D+00, & 0.8607872157664440D+00, 0.4186303330445690D+00, -0.2894719562030950D+00, & 0.4186314255109590D+00, -0.2894739329391020D+00, 0.8607860197072390D+00, & -0.8607867717364260D+00, -0.4186306871370340D+00, -0.2894727645060190D+00, & -0.2894744465036730D+00, 0.8607866519642620D+00, 0.4186297703479170D+00, & -0.4186298892623020D+00, -0.2894728382265150D+00, -0.8607871349789790D+00, & -0.2894723996931710D+00, -0.8607875560309860D+00, -0.4186293267296020D+00, & -0.4186292573884460D+00, 0.2894725949335480D+00, 0.8607875240968800D+00, & 0.2894731851561890D+00, -0.8607878175448240D+00, 0.4186282458720980D+00, & 0.4186289424246520D+00, 0.2894727563162090D+00, -0.8607876230029770D+00, & 0.2894735427627720D+00, 0.8607866031630780D+00, -0.4186304956107950D+00, & 0.8607880122617581D+00, -0.4186286769276890D+00, 0.2894719827551960D+00, & -0.8607878393611090D+00, 0.4186284794744050D+00, 0.2894727824528270D+00, & -0.1691041295942500D+00, -0.8789173916920940D+00, 0.4459910446806490D+00, & -0.8789181750666980D+00, 0.4459895744177420D+00, -0.1691039356375460D+00, & 0.1691023334883080D+00, 0.8789182061434621D+00, 0.4459901206510830D+00, & 0.4459894690345720D+00, -0.1691039825983880D+00, -0.8789182195060160D+00, & 0.8789165401716110D+00, 0.4459924431702190D+00, 0.1691048670140030D+00, & 0.4459902241080140D+00, 0.1691041980057110D+00, 0.8789177949142500D+00, & 0.8789182853189470D+00, -0.4459889689345610D+00, -0.1691049594798730D+00, & -0.4459914048704790D+00, 0.1691046511475080D+00, -0.8789171085727800D+00, & -0.8789175018706720D+00, -0.4459905576894250D+00, 0.1691048413183190D+00, & -0.4459886381463000D+00, -0.1691042398787440D+00, 0.8789185916223650D+00, & -0.1691036827918700D+00, 0.8789180817859170D+00, -0.4459898541177720D+00, & 0.1691042183067080D+00, -0.8789179365417820D+00, -0.4459899373035370D+00, & 0.6991597494364490D+00, 0.6826055934699530D+00, 0.2126622875159000D+00, & 0.6826036001105980D+00, 0.2126624328405610D+00, 0.6991616513899950D+00, & -0.6991612748012420D+00, -0.6826040568208240D+00, 0.2126622049872900D+00, & 0.2126608434125310D+00, 0.6991621013472430D+00, 0.6826036344166200D+00, & -0.6826047628202950D+00, 0.2126619851953860D+00, -0.6991606523738350D+00, & 0.2126625942230910D+00, -0.6991616996786060D+00, -0.6826035003725280D+00, & -0.6826025624027640D+00, -0.2126623681590730D+00, 0.6991626841945699D+00, & -0.2126612305468040D+00, -0.6991609500605910D+00, 0.6826046930198260D+00, & 0.6826033952274170D+00, -0.2126617699757900D+00, -0.6991620530426170D+00, & -0.2126618767979380D+00, 0.6991622129811330D+00, -0.6826031981291210D+00, & 0.6991620776118520D+00, -0.6826026485321290D+00, -0.2126640859346050D+00, & -0.6991619828086700D+00, 0.6826032705883111D+00, -0.2126624009485230D+00, & -0.8932543729812280D+00, -0.1723424150411760D+00, -0.4152044281166660D+00, & -0.1723416913831600D+00, -0.4152051051823470D+00, -0.8932541978864180D+00, & 0.8932544795129390D+00, 0.1723418657252130D+00, -0.4152044269374090D+00, & -0.4152037601443590D+00, -0.8932549611410140D+00, -0.1723409758558650D+00, & 0.1723436211169660D+00, -0.4152058956702590D+00, 0.8932534581298580D+00, & -0.4152063344497700D+00, 0.8932537022259500D+00, 0.1723412988590360D+00, & 0.1723405995636110D+00, 0.4152042778538470D+00, -0.8932547930987670D+00, & 0.4152040134618810D+00, 0.8932546529877980D+00, -0.1723419627391880D+00, & -0.1723411946292100D+00, 0.4152043287490480D+00, 0.8932546546320540D+00, & 0.4152061427713250D+00, -0.8932540159148660D+00, 0.1723401347827120D+00, & -0.8932544146433700D+00, 0.1723402748586850D+00, 0.4152052268237520D+00, & 0.8932547053896590D+00, -0.1723405507866280D+00, 0.4152044867939110D+00, & -0.3011910729024200D-01, 0.5380310043275850D+00, -0.8423867744440729D+00, & 0.5380327159133010D+00, -0.8423856691864600D+00, -0.3011944481952300D-01, & 0.3011808764135300D-01, -0.5380315902624120D+00, -0.8423864366646280D+00, & -0.8423861832095870D+00, -0.3011934729278300D-01, 0.5380319165776710D+00, & -0.5380303041055450D+00, -0.8423872336366450D+00, 0.3011877271892400D-01, & -0.8423873127238229D+00, 0.3011790102264100D-01, -0.5380302290763280D+00, & -0.5380317233086820D+00, 0.8423863249999340D+00, -0.3011883408426400D-01, & 0.8423871030981440D+00, 0.3011978965830300D-01, 0.5380304515590700D+00, & 0.5380291730323310D+00, 0.8423879687460450D+00, 0.3011841758887700D-01, & 0.8423863305324070D+00, -0.3011744117944100D-01, -0.5380317926191250D+00, & -0.3011705911664400D-01, -0.5380307391371790D+00, 0.8423870170495660D+00, & 0.3011834681252400D-01, 0.5380307101818240D+00, 0.8423869895045401D+00, & 0.9519058810513840D+00, -0.3017741210977390D+00, 0.5298654032370100D-01, & -0.3017744053434990D+00, 0.5298679853019400D-01, 0.9519057765667240D+00, & -0.9519056988554310D+00, 0.3017745999000000D+00, 0.5298708655653000D-01, & 0.5298761295842300D-01, 0.9519052380669770D+00, -0.3017759609604790D+00, & 0.3017745623980470D+00, 0.5298683421290300D-01, -0.9519057247908330D+00, & 0.5298676682920600D-01, -0.9519052521736440D+00, 0.3017760650303790D+00, & 0.3017772936453360D+00, -0.5298799472785900D-01, 0.9519047943228440D+00, & -0.5298657470118700D-01, -0.9519055918473009D+00, -0.3017750273155070D+00, & -0.3017769418414010D+00, -0.5298652631673400D-01, -0.9519049875915860D+00, & -0.5298879489611200D-01, 0.9519052401764430D+00, 0.3017757467724780D+00, & 0.9519048618122400D+00, 0.3017774620595560D+00, -0.5298582311443300D-01, & -0.9519050185943480D+00, -0.3017768243243040D+00, -0.5298663865195100D-01, & 0.5536061463002190D+00, 0.4544011566904800D+00, 0.6978823852032480D+00, & 0.4543996192985590D+00, 0.6978827381162330D+00, 0.5536069633058400D+00, & -0.5536058140188820D+00, -0.4544011281970970D+00, 0.6978826673429410D+00, & 0.6978809697722890D+00, 0.5536078912532600D+00, 0.4544012046328750D+00, & -0.4544007964543470D+00, 0.6978822688342310D+00, -0.5536065886786770D+00, & 0.6978825206532540D+00, -0.5536051035411030D+00, -0.4544022190745830D+00, & -0.4544001809008960D+00, -0.6978821410325180D+00, 0.5536072550329360D+00, & -0.6978816353908840D+00, -0.5536070488477030D+00, 0.4544012086804820D+00, & 0.4544012718187750D+00, -0.6978817257806300D+00, -0.5536068830776310D+00, & -0.6978823219163300D+00, 0.5536063431414440D+00, -0.4544010140726250D+00, & 0.5536060294802920D+00, -0.4544003068459940D+00, -0.6978830312175039D+00, & -0.5536069608106720D+00, 0.4544008427162030D+00, -0.6978819435124930D+00 & /), (/ 3, 156 /) ) return end subroutine design_18_180_3d ( x ) !*****************************************************************************80 ! !! DESIGN_18_180_3D returns an 18-design of 180 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,180), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,180) x = reshape ( (/ & -0.8663763436416970D+00, 0.2236968045802250D+00, 0.4464882650178410D+00, & 0.2236968062120170D+00, 0.4464882653478410D+00, -0.8663763430503050D+00, & 0.8663763431155790D+00, -0.2236968062252930D+00, 0.4464882652145300D+00, & 0.4464882636797900D+00, -0.8663763440671450D+00, 0.2236968056031530D+00, & -0.2236968042860020D+00, 0.4464882650235440D+00, 0.8663763437147251D+00, & 0.4464882628495670D+00, 0.8663763449479410D+00, -0.2236968038489200D+00, & -0.2236968055886000D+00, -0.4464882635825370D+00, -0.8663763441210220D+00, & -0.4464882648104650D+00, 0.8663763438297410D+00, 0.2236968042658440D+00, & 0.2236968038013990D+00, -0.4464882628087740D+00, 0.8663763449812339D+00, & -0.4464882650149240D+00, -0.8663763432190640D+00, -0.2236968062229010D+00, & -0.8663763449212000D+00, -0.2236968040744080D+00, -0.4464882627884830D+00, & 0.8663763441725580D+00, 0.2236968054822140D+00, -0.4464882635358360D+00, & -0.8068447839335680D+00, -0.4617580792431280D+00, -0.3684846956019890D+00, & -0.4617580817749450D+00, -0.3684846983908350D+00, -0.8068447812109451D+00, & 0.8068447813350600D+00, 0.4617580816135860D+00, -0.3684846983212730D+00, & -0.3684846979683570D+00, -0.8068447817064940D+00, -0.4617580812461950D+00, & 0.4617580789457650D+00, -0.3684846957167930D+00, 0.8068447840513190D+00, & -0.3684846977025540D+00, 0.8068447831055050D+00, 0.4617580790137720D+00, & 0.4617580812172950D+00, 0.3684846984438830D+00, -0.8068447815058620D+00, & 0.3684846954813280D+00, 0.8068447841829690D+00, -0.4617580789036290D+00, & -0.4617580789678360D+00, 0.3684846979360600D+00, 0.8068447830251509D+00, & 0.3684846984606100D+00, -0.8068447812133080D+00, 0.4617580817151360D+00, & -0.8068447827741030D+00, 0.4617580793147060D+00, 0.3684846980510910D+00, & 0.8068447817099870D+00, -0.4617580810987120D+00, 0.3684846981455260D+00, & -0.1348424188581120D+00, -0.4002166950757200D-01, 0.9900584770842180D+00, & -0.4002166997561800D-01, 0.9900584770162760D+00, -0.1348424192180460D+00, & 0.1348424189813570D+00, 0.4002166978894200D-01, 0.9900584770560580D+00, & 0.9900584769241170D+00, -0.1348424204361430D+00, -0.4002166815140200D-01, & 0.4002166967746100D-01, 0.9900584771169210D+00, 0.1348424185675700D+00, & 0.9900584772860210D+00, 0.1348424178553970D+00, 0.4002166789374300D-01, & 0.4002166803783600D-01, -0.9900584769275680D+00, -0.1348424204445070D+00, & -0.9900584771152390D+00, 0.1348424186351910D+00, -0.4002166949123500D-01, & -0.4002166783779800D-01, -0.9900584772708920D+00, 0.1348424179830820D+00, & -0.9900584770420310D+00, -0.1348424190874290D+00, 0.4002166977857500D-01, & -0.1348424181227450D+00, 0.4002166767089100D-01, -0.9900584772586170D+00, & 0.1348424203781130D+00, -0.4002166786721200D-01, -0.9900584769435080D+00, & 0.4979407731320700D-01, -0.2797381565618790D+00, -0.9587841851156540D+00, & -0.2797381571299750D+00, -0.9587841850685120D+00, 0.4979407502941500D-01, & -0.4979407508500500D-01, 0.2797381571118340D+00, -0.9587841850709180D+00, & -0.9587841844609550D+00, 0.4979407742976100D-01, -0.2797381587850680D+00, & 0.2797381566842330D+00, -0.9587841850838510D+00, -0.4979407723819100D-01, & -0.9587841843069630D+00, -0.4979407685685800D-01, 0.2797381594148460D+00, & 0.2797381590129380D+00, 0.9587841843903790D+00, 0.4979407750856700D-01, & 0.9587841850340860D+00, -0.4979407733711300D-01, -0.2797381568371920D+00, & -0.2797381595759920D+00, 0.9587841842431500D+00, -0.4979407718026100D-01, & 0.9587841850167220D+00, 0.4979407490928900D-01, 0.2797381573288650D+00, & 0.4979407703190500D-01, 0.2797381595171780D+00, 0.9587841842680150D+00, & -0.4979407778562100D-01, -0.2797381588886910D+00, 0.9587841844122410D+00, & 0.2054707686707770D+00, -0.1929017430722870D+00, 0.9594637464446030D+00, & -0.1929017443858980D+00, 0.9594637463317141D+00, 0.2054707679646680D+00, & -0.2054707680866780D+00, 0.1929017446304500D+00, 0.9594637462564179D+00, & 0.9594637457387351D+00, 0.2054707703405020D+00, -0.1929017448046460D+00, & 0.1929017429028800D+00, 0.9594637464476850D+00, -0.2054707688154330D+00, & 0.9594637453663400D+00, -0.2054707716940410D+00, 0.1929017452151490D+00, & 0.1929017448926260D+00, -0.9594637456856751D+00, 0.2054707705056730D+00, & -0.9594637463725330D+00, -0.2054707690642030D+00, -0.1929017430116920D+00, & -0.1929017451220650D+00, -0.9594637453488920D+00, -0.2054707718629080D+00, & -0.9594637462205630D+00, 0.2054707682606520D+00, 0.1929017446234780D+00, & 0.2054707717264440D+00, 0.1929017454605980D+00, -0.9594637453100530D+00, & -0.2054707706527430D+00, -0.1929017449496980D+00, -0.9594637456427050D+00, & -0.2789053920740190D+00, 0.7720048541378570D+00, -0.5711569726963190D+00, & 0.7720048542681720D+00, -0.5711569724663990D+00, -0.2789053921841520D+00, & 0.2789053921602380D+00, -0.7720048542493390D+00, -0.5711569725035320D+00, & -0.5711569716753651D+00, -0.2789053920788350D+00, 0.7720048548914560D+00, & -0.7720048542660130D+00, -0.5711569725335670D+00, 0.2789053920525770D+00, & -0.5711569705671390D+00, 0.2789053915822340D+00, -0.7720048558907699D+00, & -0.7720048550783650D+00, 0.5711569714219210D+00, -0.2789053920804900D+00, & 0.5711569725564700D+00, 0.2789053919520950D+00, 0.7720048542853700D+00, & 0.7720048559951140D+00, 0.5711569703761280D+00, 0.2789053916845750D+00, & 0.5711569724069060D+00, -0.2789053920341090D+00, -0.7720048543663940D+00, & -0.2789053917242620D+00, -0.7720048559398009D+00, 0.5711569704315110D+00, & 0.2789053921850090D+00, 0.7720048549702490D+00, 0.5711569715170171D+00, & 0.9123638599455530D+00, -0.3931981490415770D+00, -0.1139622861104940D+00, & -0.3931981469939110D+00, -0.1139622871432540D+00, 0.9123638606990270D+00, & -0.9123638607831750D+00, 0.3931981468241020D+00, -0.1139622870554620D+00, & -0.1139622883699460D+00, 0.9123638601627020D+00, -0.3931981478828430D+00, & 0.3931981490357560D+00, -0.1139622856815620D+00, -0.9123638600016380D+00, & -0.1139622852598250D+00, -0.9123638610299220D+00, 0.3931981467719950D+00, & 0.3931981476774950D+00, 0.1139622880664040D+00, 0.9123638602891160D+00, & 0.1139622860313020D+00, -0.9123638599292250D+00, -0.3931981491024160D+00, & -0.3931981469123200D+00, 0.1139622850908440D+00, -0.9123638609905540D+00, & 0.1139622873708760D+00, 0.9123638607624560D+00, 0.3931981467807590D+00, & 0.9123638610767000D+00, 0.3931981466410030D+00, 0.1139622853372880D+00, & -0.9123638601990820D+00, -0.3931981478492600D+00, 0.1139622881945650D+00, & 0.8486623369817879D+00, -0.1290998447282500D-01, -0.5287774296332260D+00, & -0.1290998542466000D-01, -0.5287774320354930D+00, 0.8486623354705241D+00, & -0.8486623355528130D+00, 0.1290998521418300D-01, -0.5287774319085620D+00, & -0.5287774335254600D+00, 0.8486623345656930D+00, -0.1290998387823900D-01, & 0.1290998463968200D-01, -0.5287774292818080D+00, -0.8486623371982090D+00, & -0.5287774303861490D+00, -0.8486623365210670D+00, 0.1290998392038200D-01, & 0.1290998394444800D-01, 0.5287774334464800D+00, 0.8486623346138960D+00, & 0.5287774294968270D+00, -0.8486623370671220D+00, -0.1290998444996100D-01, & -0.1290998387164700D-01, 0.5287774304196710D+00, -0.8486623365009220D+00, & 0.5287774322403560D+00, 0.8486623353443260D+00, 0.1290998532959400D-01, & 0.8486623363435590D+00, 0.1290998374355700D-01, 0.5287774306753590D+00, & -0.8486623346681990D+00, -0.1290998365530300D-01, 0.5287774333663861D+00, & -0.6958511320861700D+00, 0.2111647821010340D+00, 0.6864405558929479D+00, & 0.2111647810997110D+00, 0.6864405555444100D+00, -0.6958511327338580D+00, & 0.6958511327414010D+00, -0.2111647813350840D+00, 0.6864405554643570D+00, & 0.6864405538891910D+00, -0.6958511343847570D+00, 0.2111647810401820D+00, & -0.2111647819305030D+00, 0.6864405559602180D+00, 0.6958511320715590D+00, & 0.6864405535989390D+00, 0.6958511345259980D+00, -0.2111647815182800D+00, & -0.2111647808703600D+00, -0.6864405537379060D+00, -0.6958511345855301D+00, & -0.6864405557764750D+00, 0.6958511322246790D+00, 0.2111647820232230D+00, & 0.2111647814985050D+00, -0.6864405534995720D+00, 0.6958511346300230D+00, & -0.6864405552923321D+00, -0.6958511329810830D+00, -0.2111647811044670D+00, & -0.6958511347448820D+00, -0.2111647815311530D+00, -0.6864405533730940D+00, & 0.6958511344958130D+00, 0.2111647811014860D+00, -0.6864405537577530D+00, & -0.2617181692630290D+00, -0.5816300983962440D+00, 0.7702012909085409D+00, & -0.5816300982908330D+00, 0.7702012915065020D+00, -0.2617181677375720D+00, & 0.2617181678578640D+00, 0.5816300981264260D+00, 0.7702012915897810D+00, & 0.7702012927267941D+00, -0.2617181683217910D+00, -0.5816300964120250D+00, & 0.5816300984506260D+00, 0.7702012908883760D+00, 0.2617181692015180D+00, & 0.7702012932631270D+00, 0.2617181680777750D+00, 0.5816300958116080D+00, & 0.5816300962138030D+00, -0.7702012928812780D+00, -0.2617181683076860D+00, & -0.7702012910515680D+00, 0.2617181691302900D+00, -0.5816300982665770D+00, & -0.5816300957166070D+00, -0.7702012933042760D+00, 0.2617181681678060D+00, & -0.7702012917059650D+00, -0.2617181676410450D+00, 0.5816300980701370D+00, & -0.2617181680763480D+00, 0.5816300956377460D+00, -0.7702012933949069D+00, & 0.2617181684945420D+00, -0.5816300961299260D+00, -0.7702012928811240D+00, & 0.5061364370868440D+00, 0.7009928815969670D+00, 0.5024289870254460D+00, & 0.7009928835685090D+00, 0.5024289853021360D+00, 0.5061364360669680D+00, & -0.5061364361231960D+00, -0.7009928835031120D+00, 0.5024289853367360D+00, & 0.5024289862814260D+00, 0.5061364357644880D+00, 0.7009928830850130D+00, & -0.7009928816351710D+00, 0.5024289869389250D+00, -0.5061364371198200D+00, & 0.5024289861990811D+00, -0.5061364363423220D+00, -0.7009928827268210D+00, & -0.7009928831784340D+00, -0.5024289861247950D+00, 0.5061364357905840D+00, & -0.5024289870991430D+00, -0.5061364370344130D+00, 0.7009928815820030D+00, & 0.7009928826710060D+00, -0.5024289861979140D+00, -0.5061364364207820D+00, & -0.5024289852778980D+00, 0.5061364359559350D+00, -0.7009928836660509D+00, & 0.5061364363001890D+00, -0.7009928827898670D+00, -0.5024289861535630D+00, & -0.5061364358522460D+00, 0.7009928829915320D+00, -0.5024289863234450D+00, & -0.4407481491825780D+00, 0.6022420241579790D+00, 0.6656167165345470D+00, & 0.6022420222600990D+00, 0.6656167183423400D+00, -0.4407481490457330D+00, & 0.4407481491000160D+00, -0.6022420223379980D+00, 0.6656167182359140D+00, & 0.6656167156340270D+00, -0.4407481493907860D+00, 0.6022420250008870D+00, & -0.6022420238049980D+00, 0.6656167168141670D+00, 0.4407481492426140D+00, & 0.6656167165860120D+00, 0.4407481497832090D+00, -0.6022420236615290D+00, & -0.6022420249402080D+00, -0.6656167157609320D+00, -0.4407481492820460D+00, & -0.6656167164623710D+00, 0.4407481494240080D+00, 0.6022420240610620D+00, & 0.6022420238520260D+00, -0.6656167164607441D+00, 0.4407481497120920D+00, & -0.6656167182662930D+00, -0.4407481491798800D+00, -0.6022420222459740D+00, & -0.4407481496557820D+00, -0.6022420238831940D+00, -0.6656167164698300D+00, & 0.4407481492825400D+00, 0.6022420253069330D+00, -0.6656167154287970D+00, & -0.8902578367755301D+00, -0.2935185477582290D+00, 0.3482640466394050D+00, & -0.2935185468996730D+00, 0.3482640436499220D+00, -0.8902578382280660D+00, & 0.8902578381784459D+00, 0.2935185467624440D+00, 0.3482640438924220D+00, & 0.3482640444627600D+00, -0.8902578373223530D+00, -0.2935185486823070D+00, & 0.2935185479087850D+00, 0.3482640466866250D+00, 0.8902578367074200D+00, & 0.3482640471787870D+00, 0.8902578362705020D+00, 0.2935185486500240D+00, & 0.2935185489325450D+00, -0.3482640443361840D+00, -0.8902578372893650D+00, & -0.3482640469012240D+00, 0.8902578366266271D+00, -0.2935185478992100D+00, & -0.2935185489462000D+00, -0.3482640470802280D+00, 0.8902578362114080D+00, & -0.3482640437865890D+00, -0.8902578381927660D+00, 0.2935185468445850D+00, & -0.8902578363572190D+00, 0.2935185486920580D+00, -0.3482640469216880D+00, & 0.8902578371864430D+00, -0.2935185488110970D+00, -0.3482640447016380D+00, & 0.6619719465221540D+00, 0.3138965556450800D-01, 0.7488710379906620D+00, & 0.3138965429721000D-01, 0.7488710401727520D+00, 0.6619719441137080D+00, & -0.6619719441960080D+00, -0.3138965411214200D-01, 0.7488710401077590D+00, & 0.7488710391643290D+00, 0.6619719452186930D+00, 0.3138965505254900D-01, & -0.3138965576897200D-01, 0.7488710377831830D+00, -0.6619719467471750D+00, & 0.7488710374229330D+00, -0.6619719471714430D+00, -0.3138965541621500D-01, & -0.3138965502676800D-01, -0.7488710390441610D+00, 0.6619719453558580D+00, & -0.7488710377677350D+00, -0.6619719467611250D+00, 0.3138965584333200D-01, & 0.3138965553856000D-01, -0.7488710372221780D+00, -0.6619719473927510D+00, & -0.7488710402389310D+00, 0.6619719440450870D+00, -0.3138965416549700D-01, & 0.6619719470714800D+00, -0.3138965522335800D-01, -0.7488710375193790D+00, & -0.6619719455513510D+00, 0.3138965496147900D-01, -0.7488710388740900D+00, & -0.1257325468629560D+00, -0.8776970906645390D+00, -0.4624274469561240D+00, & -0.8776970918317050D+00, -0.4624274453820790D+00, -0.1257325445044810D+00, & 0.1257325444036380D+00, 0.8776970919764240D+00, -0.4624274451348200D+00, & -0.4624274461677560D+00, -0.1257325478951010D+00, -0.8776970909320440D+00, & 0.8776970907904780D+00, -0.4624274466873070D+00, 0.1257325469724930D+00, & -0.4624274432329320D+00, 0.1257325471315280D+00, 0.8776970925876830D+00, & 0.8776970911119200D+00, 0.4624274458629050D+00, -0.1257325477606500D+00, & 0.4624274466781810D+00, 0.1257325473667960D+00, -0.8776970907388010D+00, & -0.8776970925085100D+00, 0.4624274433572250D+00, 0.1257325472270750D+00, & 0.4624274449492740D+00, -0.1257325447342650D+00, 0.8776970920268180D+00, & -0.1257325468959420D+00, 0.8776970926167950D+00, 0.4624274432417320D+00, & 0.1257325478895730D+00, -0.8776970910219350D+00, 0.4624274459986440D+00 & /), (/ 3,180 /) ) return end subroutine design_19_204_3d ( x ) !*****************************************************************************80 ! !! DESIGN_19_204_3D returns a 19-design of 204 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,204), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,204) x = reshape ( (/ & 0.5530359455875240D+00, -0.4720502222559440D+00, 0.6865273705805380D+00, & -0.4720502274596730D+00, 0.6865273657666380D+00, 0.5530359471216960D+00, & -0.5530359455874700D+00, 0.4720502225054740D+00, 0.6865273704090060D+00, & 0.6865273723664030D+00, 0.5530359415017250D+00, -0.4720502244454320D+00, & 0.4720502285674120D+00, 0.6865273648053050D+00, -0.5530359473695520D+00, & 0.6865273720316900D+00, -0.5530359415181640D+00, 0.4720502249129640D+00, & 0.4720502283409270D+00, -0.6865273652682360D+00, 0.5530359469881980D+00, & -0.6865273717321450D+00, -0.5530359429652730D+00, -0.4720502236532300D+00, & -0.4720502275804660D+00, -0.6865273656085270D+00, -0.5530359472148680D+00, & -0.6865273710216550D+00, 0.5530359429830480D+00, 0.4720502246657080D+00, & 0.5530359466448860D+00, 0.4720502216916090D+00, -0.6865273701168060D+00, & -0.5530359472128320D+00, -0.4720502224652870D+00, -0.6865273691273200D+00, & 0.5341516544244360D+00, 0.7920823931523260D+00, 0.2954445676158600D+00, & 0.7920823974890390D+00, 0.2954445683760440D+00, 0.5341516475731480D+00, & -0.5341516546059200D+00, -0.7920823927601730D+00, 0.2954445683390990D+00, & 0.2954445679493510D+00, 0.5341516453418870D+00, 0.7920823991528760D+00, & -0.7920823976007660D+00, 0.2954445692975700D+00, -0.5341516468977650D+00, & 0.2954445678295920D+00, -0.5341516453644880D+00, -0.7920823991823051D+00, & -0.7920823978659109D+00, -0.2954445679335430D+00, 0.5341516472590450D+00, & -0.2954445679620420D+00, -0.5341516456003500D+00, 0.7920823989738450D+00, & 0.7920823971287770D+00, -0.2954445690843200D+00, -0.5341516477156180D+00, & -0.2954445674894440D+00, 0.5341516447626100D+00, -0.7920823997150640D+00, & 0.5341516544647930D+00, -0.7920823931259270D+00, -0.2954445676136700D+00, & -0.5341516544608670D+00, 0.7920823927136630D+00, -0.2954445687260430D+00, & -0.9877839019893629D+00, -0.8366313346394001D-02, -0.1556051662756040D+00, & -0.8366316491905000D-02, -0.1556051662541940D+00, -0.9877839019660940D+00, & 0.9877839020420180D+00, 0.8366312354304999D-02, -0.1556051659946880D+00, & -0.1556051675072520D+00, -0.9877839018153201D+00, -0.8366310987655000D-02, & 0.8366315777747000D-02, -0.1556051667664770D+00, 0.9877839018914431D+00, & -0.1556051684243930D+00, 0.9877839016674920D+00, 0.8366311383277999D-02, & 0.8366317026602000D-02, 0.1556051667060530D+00, -0.9877839018903840D+00, & 0.1556051668358580D+00, 0.9877839019198360D+00, -0.8366311135093001D-02, & -0.8366315838957000D-02, 0.1556051656859480D+00, 0.9877839020611400D+00, & 0.1556051677615080D+00, -0.9877839017739820D+00, 0.8366311139443000D-02, & -0.9877839021131400D+00, 0.8366313179595001D-02, 0.1556051654988360D+00, & 0.9877839021656430D+00, -0.8366312162939000D-02, 0.1556051652202080D+00, & 0.9507649813879450D+00, 0.2027274941124910D+00, -0.2344088592550430D+00, & 0.2027274967899960D+00, -0.2344088607327570D+00, 0.9507649804527050D+00, & -0.9507649809862370D+00, -0.2027274948474850D+00, -0.2344088602487210D+00, & -0.2344088502156700D+00, 0.9507649832330110D+00, 0.2027274959113820D+00, & -0.2027274972989600D+00, -0.2344088615925410D+00, -0.9507649801322030D+00, & -0.2344088501157700D+00, -0.9507649833228990D+00, -0.2027274956053300D+00, & -0.2027274967590510D+00, 0.2344088604914850D+00, 0.9507649805187890D+00, & 0.2344088505693270D+00, -0.9507649832537470D+00, 0.2027274954052100D+00, & 0.2027274975652030D+00, 0.2344088616793760D+00, -0.9507649800540250D+00, & 0.2344088503414610D+00, 0.9507649833052240D+00, -0.2027274954272670D+00, & 0.9507649813805390D+00, -0.2027274936956060D+00, 0.2344088596456210D+00, & -0.9507649809706660D+00, 0.2027274944264320D+00, 0.2344088606760230D+00, & 0.5120729891159830D+00, -0.1240516071700760D+00, -0.8499367344551850D+00, & -0.1240516096533600D+00, -0.8499367347162670D+00, 0.5120729880810549D+00, & -0.5120729889353700D+00, 0.1240516066744210D+00, -0.8499367346363440D+00, & -0.8499367347256450D+00, 0.5120729893519020D+00, -0.1240516043431770D+00, & 0.1240516097069470D+00, -0.8499367344192840D+00, -0.5120729885610040D+00, & -0.8499367341856190D+00, -0.5120729901339510D+00, 0.1240516048149250D+00, & 0.1240516099052720D+00, 0.8499367341592090D+00, 0.5120729889446310D+00, & 0.8499367344868650D+00, -0.5120729897186670D+00, -0.1240516044651950D+00, & -0.1240516097769130D+00, 0.8499367349119090D+00, -0.5120729877263990D+00, & 0.8499367338659730D+00, 0.5120729907276490D+00, 0.1240516045542460D+00, & 0.5120729896570440D+00, 0.1240516068374590D+00, 0.8499367341777510D+00, & -0.5120729893965740D+00, -0.1240516069700320D+00, 0.8499367343153300D+00, & 0.3918836979149760D+00, 0.8504231947935850D+00, -0.3510093404249470D+00, & 0.8504231953303970D+00, -0.3510093359232440D+00, 0.3918837007822100D+00, & -0.3918836974663060D+00, -0.8504231950076681D+00, -0.3510093404071850D+00, & -0.3510093358722430D+00, 0.3918837053260610D+00, 0.8504231932575950D+00, & -0.8504231945931280D+00, -0.3510093374446540D+00, -0.3918837010194270D+00, & -0.3510093379994500D+00, -0.3918837052641880D+00, -0.8504231924081080D+00, & -0.8504231944686730D+00, 0.3510093377603370D+00, 0.3918837010067490D+00, & 0.3510093355274980D+00, -0.3918837058668520D+00, 0.8504231931506850D+00, & 0.8504231953614160D+00, 0.3510093361703770D+00, -0.3918837004935390D+00, & 0.3510093368734830D+00, 0.3918837061640700D+00, -0.8504231924581730D+00, & 0.3918836983231810D+00, -0.8504231947867780D+00, 0.3510093399857000D+00, & -0.3918836981670360D+00, 0.8504231948119020D+00, 0.3510093400991560D+00, & -0.6371433781201160D+00, -0.6284993741332821D+00, 0.4461354642165980D+00, & -0.6284993752049540D+00, 0.4461354680865760D+00, -0.6371433743531780D+00, & 0.6371433775770700D+00, 0.6284993748268231D+00, 0.4461354640151090D+00, & 0.4461354669911080D+00, -0.6371433738119600D+00, -0.6284993765312260D+00, & 0.6284993759118970D+00, 0.4461354682922670D+00, 0.6371433735117989D+00, & 0.4461354673116640D+00, 0.6371433734809540D+00, 0.6284993766392390D+00, & 0.6284993752708899D+00, -0.4461354686663920D+00, -0.6371433738821430D+00, & -0.4461354670064370D+00, 0.6371433732864240D+00, -0.6284993770531080D+00, & -0.6284993761951470D+00, -0.4461354678872510D+00, 0.6371433735159890D+00, & -0.4461354676330940D+00, -0.6371433733821350D+00, 0.6284993765112530D+00, & -0.6371433778168560D+00, 0.6284993739350579D+00, -0.4461354649289460D+00, & 0.6371433775429560D+00, -0.6284993747841900D+00, -0.4461354641238870D+00, & -0.4203787081845960D+00, 0.9035659576472320D+00, -0.8276655052671900D-01, & 0.9035659605471290D+00, -0.8276654807481700D-01, -0.4203787024342720D+00, & 0.4203787076875200D+00, -0.9035659579043220D+00, -0.8276655024474300D-01, & -0.8276654593922000D-01, -0.4203787012835850D+00, 0.9035659612781000D+00, & -0.9035659607605540D+00, -0.8276654714625301D-01, 0.4203787021583580D+00, & -0.8276654503910599D-01, 0.4203787012540780D+00, -0.9035659613742790D+00, & -0.9035659605096850D+00, 0.8276654772283600D-01, -0.4203787025840560D+00, & 0.8276654605224100D-01, 0.4203787004398820D+00, 0.9035659616602750D+00, & 0.9035659609093500D+00, 0.8276654786268300D-01, 0.4203787016974780D+00, & 0.8276654567939599D-01, -0.4203787012705280D+00, -0.9035659613079750D+00, & -0.4203787075055050D+00, -0.9035659579454950D+00, 0.8276655071972200D-01, & 0.4203787069560330D+00, 0.9035659582334380D+00, 0.8276655036706200D-01, & 0.4918482984737960D+00, 0.3553670079722870D+00, 0.7948581891968250D+00, & 0.3553670126265960D+00, 0.7948581877088960D+00, 0.4918482975155830D+00, & -0.4918482983446310D+00, -0.3553670081567440D+00, 0.7948581891942840D+00, & 0.7948581920629110D+00, 0.4918482946262250D+00, 0.3553670068869010D+00, & -0.3553670125488890D+00, 0.7948581876341591D+00, -0.4918482976925080D+00, & 0.7948581920911830D+00, -0.4918482946181820D+00, -0.3553670068347960D+00, & -0.3553670126054030D+00, -0.7948581876190900D+00, 0.4918482976760280D+00, & -0.7948581926084100D+00, -0.4918482939269670D+00, 0.3553670066345830D+00, & 0.3553670125079900D+00, -0.7948581879865350D+00, -0.4918482971525960D+00, & -0.7948581923580540D+00, 0.4918482945785170D+00, -0.3553670062927780D+00, & 0.4918482979798090D+00, -0.3553670078685580D+00, -0.7948581895488740D+00, & -0.4918482975718080D+00, 0.3553670074172150D+00, -0.7948581900031270D+00, & 0.6066725580591500D-01, 0.9779826388870601D+00, 0.1996733385018680D+00, & 0.9779826388105760D+00, 0.1996733414823130D+00, 0.6066724722937100D-01, & -0.6066725576913000D-01, -0.9779826389361820D+00, 0.1996733382724510D+00, & 0.1996733379093700D+00, 0.6066725097601800D-01, 0.9779826393076430D+00, & -0.9779826390721680D+00, 0.1996733402420810D+00, -0.6066724709436200D-01, & 0.1996733373731380D+00, -0.6066725108681100D-01, -0.9779826394102520D+00, & -0.9779826389789210D+00, -0.1996733407028710D+00, 0.6066724708094300D-01, & -0.1996733379900360D+00, -0.6066725130688600D-01, 0.9779826392706490D+00, & 0.9779826388978650D+00, -0.1996733410525940D+00, -0.6066724723656200D-01, & -0.1996733370842010D+00, 0.6066725078957500D-01, -0.9779826394876820D+00, & 0.6066725570001000D-01, -0.9779826388984560D+00, -0.1996733384782320D+00, & -0.6066725607420900D-01, 0.9779826389619390D+00, -0.1996733380536040D+00, & -0.7083129618733460D+00, 0.7024145919905340D+00, 0.7004633467198600D-01, & 0.7024145841583940D+00, 0.7004632814692500D-01, -0.7083129702855400D+00, & 0.7083129618062400D+00, -0.7024145919500020D+00, 0.7004633575700700D-01, & 0.7004632573079300D-01, -0.7083129711293000D+00, 0.7024145835484910D+00, & -0.7024145842416020D+00, 0.7004632881992701D-01, 0.7083129701364700D+00, & 0.7004632513807500D-01, 0.7083129713932310D+00, -0.7024145833414500D+00, & -0.7024145843408820D+00, -0.7004632732975700D-01, -0.7083129701853820D+00, & -0.7004632609498600D-01, 0.7083129705424070D+00, 0.7024145841039930D+00, & 0.7024145841262820D+00, -0.7004632899964500D-01, 0.7083129702330579D+00, & -0.7004632593766000D-01, -0.7083129702923990D+00, -0.7024145843717901D+00, & -0.7083129612904699D+00, -0.7024145924889560D+00, -0.7004633556796400D-01, & 0.7083129610595130D+00, 0.7024145926403830D+00, -0.7004633638491400D-01, & -0.6087782464978910D+00, -0.7295294625447331D+00, -0.3117303480095350D+00, & -0.7295294611628020D+00, -0.3117303415315250D+00, -0.6087782514710520D+00, & 0.6087782466796730D+00, 0.7295294620234890D+00, -0.3117303488743800D+00, & -0.3117303430694020D+00, -0.6087782534161340D+00, -0.7295294588825280D+00, & 0.7295294609550670D+00, -0.3117303419927740D+00, 0.6087782514838040D+00, & -0.3117303424530460D+00, 0.6087782538377420D+00, 0.7295294587940750D+00, & 0.7295294612856030D+00, 0.3117303412869020D+00, -0.6087782514491540D+00, & 0.3117303426760670D+00, 0.6087782545845650D+00, -0.7295294580755680D+00, & -0.7295294607371670D+00, 0.3117303426257060D+00, 0.6087782514208260D+00, & 0.3117303425000450D+00, -0.6087782544496140D+00, 0.7295294582633970D+00, & -0.6087782472925320D+00, 0.7295294620208300D+00, 0.3117303476837500D+00, & 0.6087782473304520D+00, -0.7295294616178460D+00, 0.3117303485527810D+00, & 0.2301027741906510D+00, -0.8077565541706230D+00, 0.5427541455430510D+00, & -0.8077565520843450D+00, 0.5427541494247280D+00, 0.2301027723584630D+00, & -0.2301027736836010D+00, 0.8077565541973331D+00, 0.5427541457182660D+00, & 0.5427541442060190D+00, 0.2301027725135640D+00, -0.8077565555467580D+00, & 0.8077565521327510D+00, 0.5427541491807930D+00, -0.2301027727639210D+00, & 0.5427541438768900D+00, -0.2301027734329550D+00, 0.8077565555060050D+00, & 0.8077565522293090D+00, -0.5427541488826160D+00, 0.2301027731282830D+00, & -0.5427541450840050D+00, -0.2301027725000650D+00, -0.8077565549606600D+00, & -0.8077565522377380D+00, -0.5427541493469090D+00, -0.2301027720035430D+00, & -0.5427541442887860D+00, 0.2301027732279550D+00, 0.8077565552876390D+00, & 0.2301027740976750D+00, 0.8077565540258960D+00, -0.5427541457978590D+00, & -0.2301027735623570D+00, -0.8077565537617610D+00, -0.5427541464179090D+00, & -0.4963838094741050D+00, -0.8625182307751310D+00, -0.9831284388376600D-01, & -0.8625182242873330D+00, -0.9831283897578499D-01, -0.4963838217193900D+00, & 0.4963838095962310D+00, 0.8625182306862210D+00, -0.9831284404717300D-01, & -0.9831283935004100D-01, -0.4963838230195620D+00, -0.8625182234964180D+00, & 0.8625182243002610D+00, -0.9831283833314700D-01, 0.4963838218242060D+00, & -0.9831283829978200D-01, 0.4963838230786360D+00, 0.8625182235821330D+00, & 0.8625182244705150D+00, 0.9831283852450599D-01, -0.4963838214904720D+00, & 0.9831283917121000D-01, 0.4963838243140410D+00, -0.8625182227718220D+00, & -0.8625182240785880D+00, 0.9831283937838700D-01, 0.4963838220023670D+00, & 0.9831283847005599D-01, -0.4963838238101500D+00, 0.8625182231417350D+00, & -0.4963838100694220D+00, 0.8625182304143790D+00, 0.9831284404294301D-01, & 0.4963838104038140D+00, -0.8625182302154630D+00, 0.9831284409972600D-01, & 0.2786925513279580D+00, 0.9193131884651310D+00, 0.2778375844776740D+00, & 0.9193131917449719D+00, 0.2778375815265590D+00, 0.2786925434509230D+00, & -0.2786925515665470D+00, -0.9193131884136300D+00, 0.2778375844087580D+00, & 0.2778375830510050D+00, 0.2786925449083510D+00, 0.9193131908424260D+00, & -0.9193131921803260D+00, 0.2778375803459510D+00, -0.2786925431918220D+00, & 0.2778375820085320D+00, -0.2786925450460710D+00, -0.9193131911157350D+00, & -0.9193131921965040D+00, -0.2778375802556450D+00, 0.2786925432284890D+00, & -0.2778375828255750D+00, -0.2786925452655750D+00, 0.9193131908022630D+00, & 0.9193131918140520D+00, -0.2778375810866550D+00, -0.2786925436616070D+00, & -0.2778375815286020D+00, 0.2786925445358110D+00, -0.9193131914154680D+00, & 0.2786925512993890D+00, -0.9193131886048900D+00, -0.2778375840438940D+00, & -0.2786925517195550D+00, 0.9193131885016330D+00, -0.2778375839640920D+00, & 0.7117238189820730D+00, -0.1473551783591070D+00, -0.6868301514234280D+00, & -0.1473551800456200D+00, -0.6868301516966510D+00, 0.7117238183692320D+00, & -0.7117238189949870D+00, 0.1473551790836350D+00, -0.6868301512546030D+00, & -0.6868301560317550D+00, 0.7117238162213120D+00, -0.1473551702138960D+00, & 0.1473551798781810D+00, -0.6868301515087299D+00, -0.7117238185852460D+00, & -0.6868301558996560D+00, -0.7117238164804049D+00, 0.1473551695782020D+00, & 0.1473551798320490D+00, 0.6868301511512620D+00, 0.7117238189397620D+00, & 0.6868301563077071D+00, -0.7117238161174280D+00, -0.1473551694294310D+00, & -0.1473551804107280D+00, 0.6868301515960830D+00, -0.7117238183906890D+00, & 0.6868301558133360D+00, 0.7117238166677690D+00, 0.1473551690755790D+00, & 0.7117238191679540D+00, 0.1473551778536360D+00, 0.6868301513392560D+00, & -0.7117238189582320D+00, -0.1473551779327430D+00, 0.6868301515396070D+00, & 0.9108668157709010D+00, -0.4075474740818870D+00, 0.6501307789093599D-01, & -0.4075474700550140D+00, 0.6501307469253000D-01, 0.9108668178009230D+00, & -0.9108668160296600D+00, 0.4075474735167600D+00, 0.6501307780819900D-01, & 0.6501307141712299D-01, 0.9108668172437731D+00, -0.4075474718227450D+00, & 0.4075474695472240D+00, 0.6501307442432699D-01, -0.9108668180472660D+00, & 0.6501307150394400D-01, -0.9108668171938550D+00, 0.4075474719204620D+00, & 0.4075474699947020D+00, -0.6501307473049800D-01, 0.9108668178251990D+00, & -0.6501307123716700D-01, -0.9108668170028290D+00, -0.4075474723899620D+00, & -0.4075474694929540D+00, -0.6501307476090901D-01, -0.9108668180475250D+00, & -0.6501307089406901D-01, 0.9108668170468960D+00, 0.4075474723462040D+00, & 0.9108668155710270D+00, 0.4075474746073930D+00, -0.6501307739703199D-01, & -0.9108668158269980D+00, -0.4075474740697620D+00, -0.6501307718099700D-01 & /), (/ 3, 204 /) ) return end subroutine design_20_216_3d ( x ) !*****************************************************************************80 ! !! DESIGN_20_216_3D returns a 20-design of 216 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,216), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,216) x = reshape ( (/ & -0.2515812993559380D+00, 0.9657024628131560D+00, -0.6423085809004400D-01, & 0.9657024628129730D+00, -0.6423085809016300D-01, -0.2515812993566090D+00, & 0.2515812993562100D+00, -0.9657024628130760D+00, -0.6423085809018400D-01, & -0.6423085809003699D-01, -0.2515812993564690D+00, 0.9657024628130180D+00, & -0.9657024628129880D+00, -0.6423085809021201D-01, 0.2515812993565400D+00, & -0.6423085809028301D-01, 0.2515812993562130D+00, -0.9657024628130680D+00, & -0.9657024628131290D+00, 0.6423085809034999D-01, -0.2515812993559620D+00, & 0.6423085809020899D-01, 0.2515812993563220D+00, 0.9657024628130449D+00, & 0.9657024628130900D+00, 0.6423085808991100D-01, 0.2515812993562260D+00, & 0.6423085809020000D-01, -0.2515812993563000D+00, -0.9657024628130511D+00, & -0.2515812993566000D+00, -0.9657024628129920D+00, 0.6423085808991900D-01, & 0.2515812993565160D+00, 0.9657024628129810D+00, 0.6423085809040199D-01, & -0.7742655338457720D+00, 0.3815151823433970D+00, -0.5049346975005830D+00, & 0.3815151823431970D+00, -0.5049346975006570D+00, -0.7742655338458230D+00, & 0.7742655338455831D+00, -0.3815151823433860D+00, -0.5049346975008830D+00, & -0.5049346975007970D+00, -0.7742655338456810D+00, 0.3815151823433000D+00, & -0.3815151823431530D+00, -0.5049346975008050D+00, 0.7742655338457480D+00, & -0.5049346975006220D+00, 0.7742655338458870D+00, -0.3815151823431140D+00, & -0.3815151823432720D+00, 0.5049346975008830D+00, -0.7742655338456390D+00, & 0.5049346975008080D+00, 0.7742655338456150D+00, 0.3815151823434190D+00, & 0.3815151823434900D+00, 0.5049346975006210D+00, 0.7742655338457030D+00, & 0.5049346975005800D+00, -0.7742655338458060D+00, -0.3815151823433330D+00, & -0.7742655338457189D+00, -0.3815151823433210D+00, 0.5049346975007230D+00, & 0.7742655338458940D+00, 0.3815151823429800D+00, 0.5049346975007110D+00, & 0.6218920898658570D+00, 0.4517167996942610D+00, -0.6396891131137470D+00, & 0.4517167996941910D+00, -0.6396891131139180D+00, 0.6218920898657310D+00, & -0.6218920898656480D+00, -0.4517167996942250D+00, -0.6396891131139760D+00, & -0.6396891131139010D+00, 0.6218920898654990D+00, 0.4517167996945350D+00, & -0.4517167996940080D+00, -0.6396891131137999D+00, -0.6218920898659860D+00, & -0.6396891131138790D+00, -0.6218920898656550D+00, -0.4517167996943510D+00, & -0.4517167996943470D+00, 0.6396891131136750D+00, 0.6218920898658690D+00, & 0.6396891131137880D+00, -0.6218920898659950D+00, 0.4517167996940130D+00, & 0.4517167996945870D+00, 0.6396891131139550D+00, -0.6218920898654060D+00, & 0.6396891131140610D+00, 0.6218920898659001D+00, -0.4517167996937570D+00, & 0.6218920898658890D+00, -0.4517167996942810D+00, 0.6396891131137010D+00, & -0.6218920898658979D+00, 0.4517167996937130D+00, 0.6396891131140940D+00, & 0.2818110426750910D+00, 0.8580478476961970D+00, -0.4293441827838140D+00, & 0.8580478476964080D+00, -0.4293441827836590D+00, 0.2818110426746880D+00, & -0.2818110426751140D+00, -0.8580478476963060D+00, -0.4293441827835810D+00, & -0.4293441827833150D+00, 0.2818110426749470D+00, 0.8580478476964950D+00, & -0.8580478476963860D+00, -0.4293441827833290D+00, -0.2818110426752570D+00, & -0.4293441827839790D+00, -0.2818110426747930D+00, -0.8580478476962130D+00, & -0.8580478476961360D+00, 0.4293441827839480D+00, 0.2818110426750750D+00, & 0.4293441827835740D+00, -0.2818110426750020D+00, 0.8580478476963470D+00, & 0.8580478476964301D+00, 0.4293441827834320D+00, -0.2818110426749640D+00, & 0.4293441827834070D+00, 0.2818110426754000D+00, -0.8580478476963000D+00, & 0.2818110426747800D+00, -0.8580478476965150D+00, 0.4293441827833830D+00, & -0.2818110426751930D+00, 0.8580478476962270D+00, 0.4293441827836880D+00, & -0.6496120041073690D+00, -0.6153110840694710D+00, 0.4465383678261700D+00, & -0.6153110840695750D+00, 0.4465383678265440D+00, -0.6496120041070140D+00, & 0.6496120041073380D+00, 0.6153110840692740D+00, 0.4465383678264870D+00, & 0.4465383678262900D+00, -0.6496120041072340D+00, -0.6153110840695259D+00, & 0.6153110840696310D+00, 0.4465383678261890D+00, 0.6496120041072050D+00, & 0.4465383678263000D+00, 0.6496120041072230D+00, 0.6153110840695309D+00, & 0.6153110840693370D+00, -0.4465383678262700D+00, -0.6496120041074280D+00, & -0.4465383678262480D+00, 0.6496120041073460D+00, -0.6153110840694390D+00, & -0.6153110840693730D+00, -0.4465383678265360D+00, 0.6496120041072110D+00, & -0.4465383678262860D+00, -0.6496120041073030D+00, 0.6153110840694570D+00, & -0.6496120041071211D+00, 0.6153110840697230D+00, -0.4465383678261830D+00, & 0.6496120041071251D+00, -0.6153110840695510D+00, -0.4465383678264150D+00, & 0.9933631163195030D+00, -0.1134687281482460D+00, -0.1882994605477500D-01, & -0.1134687281480350D+00, -0.1882994605463900D-01, 0.9933631163195290D+00, & -0.9933631163195230D+00, 0.1134687281482040D+00, -0.1882994605396400D-01, & -0.1882994605390300D-01, 0.9933631163195540D+00, -0.1134687281479430D+00, & 0.1134687281480660D+00, -0.1882994605432300D-01, -0.9933631163195320D+00, & -0.1882994605474300D-01, -0.9933631163195330D+00, 0.1134687281479860D+00, & 0.1134687281482190D+00, 0.1882994605448500D-01, 0.9933631163195110D+00, & 0.1882994605434400D-01, -0.9933631163195100D+00, -0.1134687281482540D+00, & -0.1134687281481780D+00, 0.1882994605424600D-01, -0.9933631163195210D+00, & 0.1882994605448500D-01, 0.9933631163195030D+00, 0.1134687281482870D+00, & 0.9933631163195400D+00, 0.1134687281479850D+00, 0.1882994605438200D-01, & -0.9933631163195310D+00, -0.1134687281480370D+00, 0.1882994605454200D-01, & 0.2463988858915690D+00, -0.7208015696498040D+00, 0.6478677999575010D+00, & -0.7208015696493920D+00, 0.6478677999578860D+00, 0.2463988858917620D+00, & -0.2463988858916820D+00, 0.7208015696496320D+00, 0.6478677999576490D+00, & 0.6478677999574370D+00, 0.2463988858916630D+00, -0.7208015696498290D+00, & 0.7208015696498640D+00, 0.6478677999575770D+00, -0.2463988858911920D+00, & 0.6478677999576580D+00, -0.2463988858916790D+00, 0.7208015696496251D+00, & 0.7208015696496560D+00, -0.6478677999577340D+00, 0.2463988858913890D+00, & -0.6478677999579040D+00, -0.2463988858914330D+00, -0.7208015696494890D+00, & -0.7208015696498650D+00, -0.6478677999573730D+00, -0.2463988858917270D+00, & -0.6478677999574740D+00, 0.2463988858911660D+00, 0.7208015696499660D+00, & 0.2463988858917940D+00, 0.7208015696495070D+00, -0.6478677999577450D+00, & -0.2463988858914560D+00, -0.7208015696496660D+00, -0.6478677999576971D+00, & -0.7935442048021790D+00, -0.3876287734012690D+00, -0.4690751848651830D+00, & -0.3876287734013530D+00, -0.4690751848647940D+00, -0.7935442048023680D+00, & 0.7935442048021710D+00, 0.3876287734015360D+00, -0.4690751848649750D+00, & -0.4690751848650970D+00, -0.7935442048020340D+00, -0.3876287734016680D+00, & 0.3876287734016800D+00, -0.4690751848649880D+00, 0.7935442048020930D+00, & -0.4690751848651100D+00, 0.7935442048021040D+00, 0.3876287734015120D+00, & 0.3876287734014250D+00, 0.4690751848652980D+00, -0.7935442048020350D+00, & 0.4690751848650680D+00, 0.7935442048023370D+00, -0.3876287734010840D+00, & -0.3876287734014910D+00, 0.4690751848649310D+00, 0.7935442048022190D+00, & 0.4690751848647840D+00, -0.7935442048022960D+00, 0.3876287734015120D+00, & -0.7935442048022650D+00, 0.3876287734012240D+00, 0.4690751848650750D+00, & 0.7935442048021850D+00, -0.3876287734018230D+00, 0.4690751848647150D+00, & 0.1649450576530030D+00, -0.9583769097171539D+00, 0.2330382519605870D+00, & -0.9583769097169350D+00, 0.2330382519611260D+00, 0.1649450576535120D+00, & -0.1649450576532380D+00, 0.9583769097170010D+00, 0.2330382519610480D+00, & 0.2330382519606680D+00, 0.1649450576535040D+00, -0.9583769097170480D+00, & 0.9583769097171020D+00, 0.2330382519605140D+00, -0.1649450576534090D+00, & 0.2330382519607420D+00, -0.1649450576532880D+00, 0.9583769097170670D+00, & 0.9583769097170990D+00, -0.2330382519608270D+00, 0.1649450576529820D+00, & -0.2330382519611220D+00, -0.1649450576532260D+00, -0.9583769097169860D+00, & -0.9583769097170930D+00, -0.2330382519606320D+00, -0.1649450576532930D+00, & -0.2330382519604340D+00, 0.1649450576532610D+00, 0.9583769097171470D+00, & 0.1649450576534940D+00, 0.9583769097169650D+00, -0.2330382519610150D+00, & -0.1649450576534580D+00, -0.9583769097170310D+00, -0.2330382519607690D+00, & 0.5604842504669760D+00, 0.8132526494836950D+00, -0.1564529740408340D+00, & 0.8132526494836900D+00, -0.1564529740414460D+00, 0.5604842504668130D+00, & -0.5604842504672400D+00, -0.8132526494834310D+00, -0.1564529740412630D+00, & -0.1564529740410300D+00, 0.5604842504670470D+00, 0.8132526494836090D+00, & -0.8132526494838200D+00, -0.1564529740407260D+00, -0.5604842504668260D+00, & -0.1564529740410970D+00, -0.5604842504667780D+00, -0.8132526494837810D+00, & -0.8132526494836300D+00, 0.1564529740409670D+00, 0.5604842504670350D+00, & 0.1564529740412850D+00, -0.5604842504670530D+00, 0.8132526494835550D+00, & 0.8132526494834810D+00, 0.1564529740411510D+00, -0.5604842504671990D+00, & 0.1564529740408810D+00, 0.5604842504669960D+00, -0.8132526494836720D+00, & 0.5604842504668360D+00, -0.8132526494837370D+00, 0.1564529740411220D+00, & -0.5604842504667400D+00, 0.8132526494838230D+00, 0.1564529740410180D+00, & 0.3666300586513120D+00, 0.9220188325509330D+00, -0.1243530157042820D+00, & 0.9220188325508800D+00, -0.1243530157047620D+00, 0.3666300586512840D+00, & -0.3666300586517610D+00, -0.9220188325507080D+00, -0.1243530157046290D+00, & -0.1243530157043770D+00, 0.3666300586515770D+00, 0.9220188325508150D+00, & -0.9220188325509330D+00, -0.1243530157042030D+00, -0.3666300586513410D+00, & -0.1243530157045340D+00, -0.3666300586511110D+00, -0.9220188325509791D+00, & -0.9220188325508830D+00, 0.1243530157044780D+00, 0.3666300586513720D+00, & 0.1243530157046300D+00, -0.3666300586515370D+00, 0.9220188325507970D+00, & 0.9220188325507450D+00, 0.1243530157044630D+00, -0.3666300586517230D+00, & 0.1243530157042990D+00, 0.3666300586515630D+00, -0.9220188325508309D+00, & 0.3666300586512860D+00, -0.9220188325509230D+00, 0.1243530157044380D+00, & -0.3666300586512290D+00, 0.9220188325509380D+00, 0.1243530157044920D+00, & -0.8046719536517351D+00, -0.7083625075572700D-01, 0.5894788143650050D+00, & -0.7083625075605800D-01, 0.5894788143650030D+00, -0.8046719536517070D+00, & 0.8046719536519210D+00, 0.7083625075538300D-01, 0.5894788143647920D+00, & 0.5894788143647260D+00, -0.8046719536519410D+00, -0.7083625075571399D-01, & 0.7083625075593900D-01, 0.5894788143647760D+00, 0.8046719536518840D+00, & 0.5894788143650180D+00, 0.8046719536517150D+00, 0.7083625075584600D-01, & 0.7083625075560100D-01, -0.5894788143648110D+00, -0.8046719536518880D+00, & -0.5894788143647840D+00, 0.8046719536518840D+00, -0.7083625075587500D-01, & -0.7083625075555100D-01, -0.5894788143649440D+00, 0.8046719536517950D+00, & -0.5894788143649780D+00, -0.8046719536517590D+00, 0.7083625075567000D-01, & -0.8046719536518360D+00, 0.7083625075619300D-01, -0.5894788143648110D+00, & 0.8046719536517311D+00, -0.7083625075576400D-01, -0.5894788143650060D+00, & -0.8305971377714630D+00, -0.4813562216367220D+00, 0.2800081831259090D+00, & -0.4813562216367630D+00, 0.2800081831263240D+00, -0.8305971377712990D+00, & 0.8305971377714670D+00, 0.4813562216366280D+00, 0.2800081831260560D+00, & 0.2800081831258640D+00, -0.8305971377713430D+00, -0.4813562216369560D+00, & 0.4813562216370750D+00, 0.2800081831258990D+00, 0.8305971377712620D+00, & 0.2800081831260040D+00, 0.8305971377713510D+00, 0.4813562216368590D+00, & 0.4813562216366530D+00, -0.2800081831258590D+00, -0.8305971377715200D+00, & -0.2800081831260120D+00, 0.8305971377715200D+00, -0.4813562216365640D+00, & -0.4813562216367410D+00, -0.2800081831261120D+00, 0.8305971377713840D+00, & -0.2800081831260530D+00, -0.8305971377713141D+00, 0.4813562216368940D+00, & -0.8305971377713660D+00, 0.4813562216368400D+00, -0.2800081831259940D+00, & 0.8305971377711940D+00, -0.4813562216370290D+00, -0.2800081831261780D+00, & 0.6225761054046420D+00, 0.2744190843023600D-01, -0.7820779594393990D+00, & 0.2744190843027600D-01, -0.7820779594394310D+00, 0.6225761054046010D+00, & -0.6225761054049630D+00, -0.2744190843004500D-01, -0.7820779594391500D+00, & -0.7820779594391180D+00, 0.6225761054049880D+00, 0.2744190843039700D-01, & -0.2744190843020100D-01, -0.7820779594392960D+00, -0.6225761054047740D+00, & -0.7820779594394081D+00, -0.6225761054046280D+00, -0.2744190843028900D-01, & -0.2744190843023800D-01, 0.7820779594392210D+00, 0.6225761054048660D+00, & 0.7820779594392629D+00, -0.6225761054048200D+00, 0.2744190843008300D-01, & 0.2744190843041900D-01, 0.7820779594392691D+00, -0.6225761054047980D+00, & 0.7820779594394510D+00, 0.6225761054045910D+00, -0.2744190842992800D-01, & 0.6225761054047880D+00, -0.2744190843038000D-01, 0.7820779594392780D+00, & -0.6225761054045720D+00, 0.2744190842986800D-01, 0.7820779594394680D+00, & -0.9318695934738700D+00, 0.3187128632820320D+00, -0.1733238919982290D+00, & 0.3187128632819440D+00, -0.1733238919982580D+00, -0.9318695934738940D+00, & 0.9318695934737440D+00, -0.3187128632820510D+00, -0.1733238919988710D+00, & -0.1733238919988410D+00, -0.9318695934738360D+00, 0.3187128632817990D+00, & -0.3187128632819240D+00, -0.1733238919986170D+00, 0.9318695934738340D+00, & -0.1733238919982450D+00, 0.9318695934739750D+00, -0.3187128632817140D+00, & -0.3187128632819970D+00, 0.1733238919985150D+00, -0.9318695934738280D+00, & 0.1733238919985010D+00, 0.9318695934738011D+00, 0.3187128632820840D+00, & 0.3187128632820890D+00, 0.1733238919985390D+00, 0.9318695934737929D+00, & 0.1733238919984430D+00, -0.9318695934738240D+00, -0.3187128632820500D+00, & -0.9318695934738650D+00, -0.3187128632819280D+00, 0.1733238919984480D+00, & 0.9318695934738970D+00, 0.3187128632818020D+00, 0.1733238919985030D+00, & 0.8838481768527030D+00, 0.2014238044752130D+00, 0.4221858018276850D+00, & 0.2014238044757030D+00, 0.4221858018276610D+00, 0.8838481768526020D+00, & -0.8838481768525340D+00, -0.2014238044755540D+00, 0.4221858018278750D+00, & 0.4221858018279100D+00, 0.8838481768524840D+00, 0.2014238044757010D+00, & -0.2014238044754720D+00, 0.4221858018277440D+00, -0.8838481768526150D+00, & 0.4221858018276230D+00, -0.8838481768526470D+00, -0.2014238044755860D+00, & -0.2014238044753970D+00, -0.4221858018278330D+00, 0.8838481768525900D+00, & -0.4221858018279300D+00, -0.8838481768525230D+00, 0.2014238044754890D+00, & 0.2014238044754790D+00, -0.4221858018276820D+00, -0.8838481768526431D+00, & -0.4221858018275140D+00, 0.8838481768527690D+00, -0.2014238044752800D+00, & 0.8838481768524760D+00, -0.2014238044756140D+00, -0.4221858018279670D+00, & -0.8838481768527100D+00, 0.2014238044755630D+00, -0.4221858018275020D+00, & 0.2042750399564050D+00, 0.7187705698842261D+00, 0.6645604381236631D+00, & 0.7187705698843340D+00, 0.6645604381234740D+00, 0.2042750399566370D+00, & -0.2042750399562600D+00, -0.7187705698842650D+00, 0.6645604381236641D+00, & 0.6645604381238100D+00, 0.2042750399561560D+00, 0.7187705698841600D+00, & -0.7187705698843250D+00, 0.6645604381235790D+00, -0.2042750399563280D+00, & 0.6645604381234920D+00, -0.2042750399563730D+00, -0.7187705698843930D+00, & -0.7187705698843611D+00, -0.6645604381235540D+00, 0.2042750399562800D+00, & -0.6645604381235540D+00, -0.2042750399566620D+00, 0.7187705698842540D+00, & 0.7187705698840900D+00, -0.6645604381238021D+00, -0.2042750399564320D+00, & -0.6645604381235050D+00, 0.2042750399566820D+00, -0.7187705698842930D+00, & 0.2042750399563760D+00, -0.7187705698841650D+00, -0.6645604381237380D+00, & -0.2042750399563670D+00, 0.7187705698845380D+00, -0.6645604381233370D+00, & -0.8988479274720690D+00, 0.4377008233682800D+00, 0.2214480756061700D-01, & 0.4377008233679230D+00, 0.2214480756096300D-01, -0.8988479274722340D+00, & 0.8988479274721820D+00, -0.4377008233680650D+00, 0.2214480756027000D-01, & 0.2214480756031500D-01, -0.8988479274722930D+00, 0.4377008233678340D+00, & -0.4377008233677660D+00, 0.2214480756062300D-01, 0.8988479274723190D+00, & 0.2214480756055900D-01, 0.8988479274721600D+00, -0.4377008233680940D+00, & -0.4377008233682550D+00, -0.2214480756032700D-01, -0.8988479274720880D+00, & -0.2214480756066100D-01, 0.8988479274721030D+00, 0.4377008233682070D+00, & 0.4377008233680300D+00, -0.2214480756060700D-01, 0.8988479274721910D+00, & -0.2214480756073300D-01, -0.8988479274721950D+00, -0.4377008233680150D+00, & -0.8988479274722450D+00, -0.4377008233679080D+00, -0.2214480756079600D-01, & 0.8988479274723130D+00, 0.4377008233677780D+00, -0.2214480756063400D-01 & /), (/ 3, 216 /) ) return end subroutine design_21_240_3d ( x ) !*****************************************************************************80 ! !! DESIGN_21_240_3D returns a 21-design of 240 points in 3D. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 March 2020 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Ronald Hardin, Neil Sloane, ! McLaren's Improved Snub Cube and Other New Spherical Designs ! in Three Dimensions, ! Discrete and Computational Geometry, ! Volume 15, 1996, pages 429-441. ! ! Output: ! ! real ( kind = rk8 ) X(3,240), the coordinates of the points. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) x(3,240) x = reshape ( (/ & 0.8926535357627230D+00, 0.4125340536573610D+00, -0.1816186104542530D+00, & 0.4125340534250320D+00, -0.1816186106417820D+00, 0.8926535358319380D+00, & -0.8926535358064071D+00, -0.4125340536278530D+00, -0.1816186103065750D+00, & -0.1816186106138490D+00, 0.8926535357404750D+00, 0.4125340536352400D+00, & -0.4125340534774350D+00, -0.1816186104226540D+00, -0.8926535358523040D+00, & -0.1816186104513840D+00, -0.8926535357628120D+00, -0.4125340536584320D+00, & -0.4125340533170900D+00, 0.1816186106118270D+00, 0.8926535358879180D+00, & 0.1816186104001360D+00, -0.8926535358123000D+00, 0.4125340535739110D+00, & 0.4125340533279960D+00, 0.1816186104204000D+00, -0.8926535359218249D+00, & 0.1816186105807890D+00, 0.8926535358109040D+00, -0.4125340534973990D+00, & 0.8926535358676440D+00, -0.4125340534725580D+00, 0.1816186103583390D+00, & -0.8926535358550640D+00, 0.4125340535351600D+00, 0.1816186102779710D+00, & -0.2920937425934330D+00, -0.2957670279931700D+00, 0.9095070701703470D+00, & -0.2957670280268870D+00, 0.9095070700892600D+00, -0.2920937428117760D+00, & 0.2920937424478640D+00, 0.2957670280397130D+00, 0.9095070702019620D+00, & 0.9095070701476120D+00, -0.2920937429267210D+00, -0.2957670277339340D+00, & 0.2957670281453960D+00, 0.9095070700844410D+00, 0.2920937427067830D+00, & 0.9095070701888540D+00, 0.2920937426892070D+00, 0.2957670278416750D+00, & 0.2957670279073110D+00, -0.9095070701484190D+00, -0.2920937427486510D+00, & -0.9095070701012210D+00, 0.2920937431592720D+00, -0.2957670276469270D+00, & -0.2957670278353330D+00, -0.9095070700472930D+00, 0.2920937431364140D+00, & -0.9095070702185910D+00, -0.2920937427217760D+00, 0.2957670277180690D+00, & -0.2920937425408960D+00, 0.2957670277931470D+00, -0.9095070702522660D+00, & 0.2920937428619380D+00, -0.2957670277476140D+00, -0.9095070701639690D+00, & -0.5752257180381920D+00, 0.2412057282507800D-01, 0.8176390225974030D+00, & 0.2412057278614400D-01, 0.8176390225112380D+00, -0.5752257181623009D+00, & 0.5752257181164780D+00, -0.2412057297921300D-01, 0.8176390225377810D+00, & 0.8176390225560030D+00, -0.5752257181034800D+00, 0.2412057267146900D-01, & -0.2412057304150300D-01, 0.8176390224407570D+00, 0.5752257182517770D+00, & 0.8176390224583791D+00, 0.5752257182291181D+00, -0.2412057298452600D-01, & -0.2412057281823900D-01, -0.8176390225812600D+00, -0.5752257180614240D+00, & -0.8176390225435780D+00, 0.5752257181238820D+00, 0.2412057260611100D-01, & 0.2412057271295000D-01, -0.8176390225272960D+00, 0.5752257181425460D+00, & -0.8176390226004950D+00, -0.5752257180351740D+00, -0.2412057279222800D-01, & -0.5752257179254689D+00, -0.2412057271105200D-01, -0.8176390226800700D+00, & 0.5752257179082300D+00, 0.2412057259415500D-01, -0.8176390226956460D+00, & -0.1288331617248000D+00, 0.5224764072024000D-01, 0.9902889479738530D+00, & 0.5224764069440900D-01, 0.9902889479588950D+00, -0.1288331618502510D+00, & 0.1288331618403250D+00, -0.5224764032003800D-01, 0.9902889479799381D+00, & 0.9902889479497170D+00, -0.1288331619247960D+00, 0.5224764068455800D-01, & -0.5224764038851000D-01, 0.9902889479675810D+00, 0.1288331619075380D+00, & 0.9902889479777300D+00, 0.1288331618780010D+00, -0.5224764026899200D-01, & -0.5224764039040900D-01, -0.9902889479621900D+00, -0.1288331619482090D+00, & -0.9902889479606261D+00, 0.1288331618966490D+00, 0.5224764054718700D-01, & 0.5224764052780800D-01, -0.9902889479532510D+00, 0.1288331619612000D+00, & -0.9902889479708680D+00, -0.1288331619362050D+00, -0.5224764025552600D-01, & -0.1288331617904780D+00, -0.5224764033764300D-01, -0.9902889479854941D+00, & 0.1288331618574160D+00, 0.5224764055154500D-01, -0.9902889479655000D+00, & 0.7180063860347500D+00, 0.6574468762559930D+00, -0.2285397875962860D+00, & 0.6574468762867370D+00, -0.2285397878319220D+00, 0.7180063859315960D+00, & -0.7180063861094420D+00, -0.6574468761714340D+00, -0.2285397876048770D+00, & -0.2285397877372190D+00, 0.7180063859474221D+00, 0.6574468763023740D+00, & -0.6574468762410210D+00, -0.2285397877138000D+00, -0.7180063860110540D+00, & -0.2285397876789970D+00, -0.7180063860313590D+00, -0.6574468762309450D+00, & -0.6574468763611850D+00, 0.2285397878605490D+00, 0.7180063858543150D+00, & 0.2285397877030650D+00, -0.7180063858573850D+00, 0.6574468764125770D+00, & 0.6574468763044540D+00, 0.2285397878740170D+00, -0.7180063859019750D+00, & 0.2285397877849670D+00, 0.7180063858138530D+00, -0.6574468764316480D+00, & 0.7180063858807600D+00, -0.6574468763634850D+00, 0.2285397877708510D+00, & -0.7180063858910180D+00, 0.6574468763715580D+00, 0.2285397877154010D+00, & 0.8631764731178030D+00, 0.4681818166531380D+00, 0.1890295289400010D+00, & 0.4681818164384860D+00, 0.1890295291974920D+00, 0.8631764731778400D+00, & -0.8631764731944460D+00, -0.4681818165764200D+00, 0.1890295287800330D+00, & 0.1890295291255270D+00, 0.8631764730643890D+00, 0.4681818166767080D+00, & -0.4681818163926710D+00, 0.1890295288974430D+00, -0.8631764732683980D+00, & 0.1890295287921740D+00, -0.8631764731436880D+00, -0.4681818166651000D+00, & -0.4681818164112130D+00, -0.1890295291281380D+00, 0.8631764732078210D+00, & -0.1890295288978520D+00, -0.8631764730897200D+00, 0.4681818167219310D+00, & 0.4681818165088670D+00, -0.1890295289305550D+00, -0.8631764731981230D+00, & -0.1890295290018230D+00, 0.8631764731066590D+00, -0.4681818166487220D+00, & 0.8631764731352290D+00, -0.4681818166486420D+00, -0.1890295288715610D+00, & -0.8631764731233340D+00, 0.4681818166987620D+00, -0.1890295288017440D+00, & 0.7726328568471330D+00, -0.5170594506955900D+00, 0.3683585114621520D+00, & -0.5170594505671320D+00, 0.3683585115855150D+00, 0.7726328568742860D+00, & -0.7726328568060810D+00, 0.5170594506473910D+00, 0.3683585116159150D+00, & 0.3683585116480010D+00, 0.7726328568060540D+00, -0.5170594506245730D+00, & 0.5170594504940070D+00, 0.3683585118165880D+00, -0.7726328568130560D+00, & 0.3683585117204960D+00, -0.7726328568024760D+00, 0.5170594505782730D+00, & 0.5170594505834450D+00, -0.3683585114871170D+00, 0.7726328569102801D+00, & -0.3683585115673300D+00, -0.7726328568594670D+00, -0.5170594506022290D+00, & -0.5170594505023690D+00, -0.3683585116659560D+00, -0.7726328568792750D+00, & -0.3683585114698030D+00, 0.7726328568556510D+00, 0.5170594506774120D+00, & 0.7726328569347490D+00, 0.5170594506919191D+00, -0.3683585112835310D+00, & -0.7726328569274850D+00, -0.5170594506337780D+00, -0.3683585113803780D+00, & -0.8478192319146480D+00, -0.6632577590016700D-01, -0.5261211281130020D+00, & -0.6632577591363099D-01, -0.5261211282576860D+00, -0.8478192318238090D+00, & 0.8478192318830180D+00, 0.6632577581985200D-01, -0.5261211281740970D+00, & -0.5261211283487620D+00, -0.8478192317669569D+00, -0.6632577591791000D-01, & 0.6632577584612000D-01, -0.5261211284070980D+00, 0.8478192317363720D+00, & -0.5261211284592400D+00, 0.8478192317090800D+00, 0.6632577578136600D-01, & 0.6632577594578500D-01, 0.5261211283443800D+00, -0.8478192317674960D+00, & 0.5261211284495320D+00, 0.8478192317006920D+00, -0.6632577596561300D-01, & -0.6632577587721100D-01, 0.5261211283063880D+00, 0.8478192317964360D+00, & 0.5261211285046690D+00, -0.8478192316652130D+00, 0.6632577598176000D-01, & -0.8478192318217250D+00, 0.6632577594100500D-01, 0.5261211282575941D+00, & 0.8478192318502640D+00, -0.6632577599665500D-01, 0.5261211282045900D+00, & 0.9805743229230000D-02, 0.9429838158425931D+00, 0.3326941094438920D+00, & 0.9429838158089230D+00, 0.3326941095397480D+00, 0.9805743214949999D-02, & -0.9805743379690000D-02, -0.9429838157872910D+00, 0.3326941095962070D+00, & 0.3326941092265540D+00, 0.9805743204272001D-02, 0.9429838159195320D+00, & -0.9429838157740400D+00, 0.3326941096356470D+00, -0.9805743315804000D-02, & 0.3326941093979960D+00, -0.9805743298910000D-02, -0.9429838158580620D+00, & -0.9429838157761140D+00, -0.3326941096300980D+00, 0.9805743304667000D-02, & -0.3326941093190270D+00, -0.9805743188507000D-02, 0.9429838158870700D+00, & 0.9429838157750819D+00, -0.3326941096351990D+00, -0.9805743230762999D-02, & -0.3326941094557650D+00, 0.9805743389761999D-02, -0.9429838158367350D+00, & 0.9805743301140001D-02, -0.9429838157525240D+00, -0.3326941096970650D+00, & -0.9805743287713000D-02, 0.9429838157913790D+00, -0.3326941095873310D+00, & 0.7855992483711520D+00, -0.4051569453122690D+00, -0.4676341204658960D+00, & -0.4051569449321250D+00, -0.4676341206498590D+00, 0.7855992484576980D+00, & -0.7855992482017899D+00, 0.4051569454340510D+00, -0.4676341206449040D+00, & -0.4676341206112420D+00, 0.7855992483346230D+00, -0.4051569452153390D+00, & 0.4051569451364230D+00, -0.4676341208682010D+00, -0.7855992482223660D+00, & -0.4676341208118040D+00, -0.7855992481456090D+00, 0.4051569453503470D+00, & 0.4051569448419850D+00, 0.4676341208613320D+00, 0.7855992483783050D+00, & 0.4676341207867260D+00, -0.7855992482498571D+00, -0.4051569451771560D+00, & -0.4051569449996430D+00, 0.4676341208710980D+00, -0.7855992482911820D+00, & 0.4676341208937130D+00, 0.7855992482342400D+00, 0.4051569450839530D+00, & 0.7855992483133410D+00, 0.4051569451171040D+00, 0.4676341207321060D+00, & -0.7855992482810999D+00, -0.4051569451973700D+00, 0.4676341207167270D+00, & -0.7373319991314921D+00, 0.6208515010137640D+00, -0.2662422519918900D+00, & 0.6208515009491860D+00, -0.2662422521548950D+00, -0.7373319991270100D+00, & 0.7373319990600610D+00, -0.6208515010887370D+00, -0.2662422520148830D+00, & -0.2662422519486310D+00, -0.7373319991032550D+00, 0.6208515010658500D+00, & -0.6208515010792210D+00, -0.2662422522338000D+00, 0.7373319989890250D+00, & -0.2662422520116240D+00, 0.7373319989962220D+00, -0.6208515011659510D+00, & -0.6208515010721239D+00, 0.2662422522225600D+00, -0.7373319989990600D+00, & 0.2662422521138640D+00, 0.7373319988329740D+00, 0.6208515013159830D+00, & 0.6208515011873870D+00, 0.2662422523283740D+00, 0.7373319988637970D+00, & 0.2662422519322500D+00, -0.7373319989389900D+00, -0.6208515012679590D+00, & -0.7373319989479430D+00, -0.6208515011832970D+00, 0.2662422521048790D+00, & 0.7373319988350070D+00, 0.6208515013057860D+00, 0.2662422521320100D+00, & 0.7268714691656590D+00, -0.2748828235042800D-01, -0.6862231864680610D+00, & -0.2748828218275500D-01, -0.6862231864483250D+00, 0.7268714691906329D+00, & -0.7268714691729310D+00, 0.2748828237188500D-01, -0.6862231864594990D+00, & -0.6862231864497120D+00, 0.7268714691854060D+00, -0.2748828228634100D-01, & 0.2748828235160700D-01, -0.6862231864911200D+00, -0.7268714691438450D+00, & -0.6862231865456220D+00, -0.7268714690897941D+00, 0.2748828242028100D-01, & 0.2748828226683600D-01, 0.6862231864703350D+00, 0.7268714691666740D+00, & 0.6862231866611830D+00, -0.7268714689834220D+00, -0.2748828234818500D-01, & -0.2748828225102900D-01, 0.6862231865230920D+00, -0.7268714691174650D+00, & 0.6862231866091120D+00, 0.7268714690334980D+00, 0.2748828232394800D-01, & 0.7268714690701070D+00, 0.2748828233555000D-01, 0.6862231865698690D+00, & -0.7268714690801830D+00, -0.2748828230971600D-01, 0.6862231865602320D+00, & 0.6653633857205150D+00, 0.5808602677392710D+00, 0.4689274083527160D+00, & 0.5808602675770870D+00, 0.4689274084886380D+00, 0.6653633857663080D+00, & -0.6653633856773800D+00, -0.5808602677195750D+00, 0.4689274084383180D+00, & 0.4689274083407830D+00, 0.6653633858218631D+00, 0.5808602676328130D+00, & -0.5808602675284530D+00, 0.4689274086788320D+00, -0.6653633856747230D+00, & 0.4689274083726140D+00, -0.6653633856988030D+00, -0.5808602677480780D+00, & -0.5808602676408769D+00, -0.4689274085527620D+00, 0.6653633856654270D+00, & -0.4689274084683360D+00, -0.6653633858479470D+00, 0.5808602674999610D+00, & 0.5808602673867520D+00, -0.4689274086545190D+00, -0.6653633858155630D+00, & -0.4689274083756990D+00, 0.6653633856513560D+00, -0.5808602677999380D+00, & 0.6653633856518190D+00, -0.5808602677912120D+00, -0.4689274083858500D+00, & -0.6653633857517340D+00, 0.5808602675480170D+00, -0.4689274085453260D+00, & -0.5801253673053040D+00, -0.7790995979244340D+00, 0.2376097109187070D+00, & -0.7790995980535180D+00, 0.2376097109099340D+00, -0.5801253671355390D+00, & 0.5801253671868080D+00, 0.7790995979777320D+00, 0.2376097110332580D+00, & 0.2376097106959320D+00, -0.5801253672761100D+00, -0.7790995980141140D+00, & 0.7790995980647319D+00, 0.2376097111473200D+00, 0.5801253670232500D+00, & 0.2376097108192850D+00, 0.5801253670474260D+00, 0.7790995981467740D+00, & 0.7790995981702240D+00, -0.2376097108496420D+00, -0.5801253670034990D+00, & -0.2376097108118020D+00, 0.5801253671572560D+00, -0.7790995980672760D+00, & -0.7790995980749610D+00, -0.2376097110451280D+00, 0.5801253670513690D+00, & -0.2376097106092530D+00, -0.5801253670223590D+00, 0.7790995982294950D+00, & -0.5801253670900940D+00, 0.7790995981519661D+00, -0.2376097106980860D+00, & 0.5801253672184110D+00, -0.7790995979667160D+00, -0.2376097109922150D+00, & 0.9586680253602000D+00, 0.1011136059005390D+00, -0.2659542363899560D+00, & 0.1011136058898930D+00, -0.2659542364771990D+00, 0.9586680253371200D+00, & -0.9586680253264100D+00, -0.1011136060954320D+00, -0.2659542364376600D+00, & -0.2659542366341790D+00, 0.9586680252945550D+00, 0.1011136058805580D+00, & -0.1011136060031710D+00, -0.2659542366563170D+00, -0.9586680252754820D+00, & -0.2659542367154550D+00, -0.9586680252461620D+00, -0.1011136061256020D+00, & -0.1011136058254380D+00, 0.2659542364146640D+00, 0.9586680253612671D+00, & 0.2659542362867390D+00, -0.9586680253935830D+00, 0.1011136058555220D+00, & 0.1011136058024440D+00, 0.2659542362606640D+00, -0.9586680254064150D+00, & 0.2659542365158540D+00, 0.9586680253225770D+00, -0.1011136059261060D+00, & 0.9586680254495000D+00, -0.1011136059091010D+00, 0.2659542360648080D+00, & -0.9586680254786000D+00, 0.1011136057894970D+00, 0.2659542360053860D+00, & -0.7844318144170850D+00, 0.2843190250072290D+00, 0.5512072392025160D+00, & 0.2843190248228480D+00, 0.5512072393207090D+00, -0.7844318144008620D+00, & 0.7844318144434220D+00, -0.2843190248881310D+00, 0.5512072392264670D+00, & 0.5512072394346770D+00, -0.7844318142918880D+00, 0.2843190249025560D+00, & -0.2843190246401610D+00, 0.5512072393475040D+00, 0.7844318144482491D+00, & 0.5512072394083570D+00, 0.7844318144009980D+00, -0.2843190246525460D+00, & -0.2843190247149400D+00, -0.5512072391601370D+00, -0.7844318145528040D+00, & -0.5512072394176490D+00, 0.7844318144267430D+00, 0.2843190245635030D+00, & 0.2843190244771060D+00, -0.5512072393940670D+00, 0.7844318144746290D+00, & -0.5512072392271640D+00, -0.7844318145108320D+00, -0.2843190247007970D+00, & -0.7844318146549000D+00, -0.2843190247577290D+00, -0.5512072389927720D+00, & 0.7844318145421389D+00, 0.2843190246898840D+00, -0.5512072391882400D+00, & 0.1666638785351180D+00, 0.9794687788666500D+00, 0.1134198519532850D+00, & 0.9794687788923619D+00, 0.1134198520112480D+00, 0.1666638783445640D+00, & -0.1666638783223350D+00, -0.9794687788772219D+00, 0.1134198521746590D+00, & 0.1134198518526030D+00, 0.1666638784650920D+00, 0.9794687788902240D+00, & -0.9794687789080510D+00, 0.1134198522332290D+00, -0.1666638781012970D+00, & 0.1134198520235320D+00, -0.1666638782131650D+00, -0.9794687789132980D+00, & -0.9794687788914180D+00, -0.1134198520887550D+00, 0.1666638782973680D+00, & -0.1134198519422990D+00, -0.1666638783837850D+00, 0.9794687788936730D+00, & 0.9794687788877920D+00, -0.1134198522526510D+00, -0.1666638782071420D+00, & -0.1134198518873330D+00, 0.1666638784200610D+00, -0.9794687788938650D+00, & 0.1666638785133120D+00, -0.9794687788588400D+00, -0.1134198520527750D+00, & -0.1666638785259920D+00, 0.9794687788524030D+00, -0.1134198520897270D+00, & 0.9035426353908700D+00, 0.9900269067959901D-01, 0.4169042735078650D+00, & 0.9900269051118001D-01, 0.4169042737536920D+00, 0.9035426352958970D+00, & -0.9035426353835330D+00, -0.9900269064792300D-01, 0.4169042735312880D+00, & 0.4169042739582500D+00, 0.9035426351937680D+00, 0.9900269058185000D-01, & -0.9900269041493300D-01, 0.4169042736997320D+00, -0.9035426353313410D+00, & 0.4169042738439640D+00, -0.9035426352375170D+00, -0.9900269066384500D-01, & -0.9900269046419199D-01, -0.4169042739372540D+00, 0.9035426352163480D+00, & -0.4169042742060360D+00, -0.9035426351101470D+00, 0.9900269030157501D-01, & 0.9900269012804400D-01, -0.4169042740643800D+00, -0.9035426351945230D+00, & -0.4169042741137440D+00, 0.9035426351313860D+00, -0.9900269049639200D-01, & 0.9035426352792750D+00, -0.9900269046710200D-01, -0.4169042738001830D+00, & -0.9035426352343990D+00, 0.9900269024582900D-01, -0.4169042739499880D+00, & 0.2787624045360920D+00, 0.3493121855370630D+00, -0.8945795206981750D+00, & 0.3493121855860560D+00, -0.8945795206085150D+00, 0.2787624047624310D+00, & -0.2787624045405250D+00, -0.3493121855034730D+00, -0.8945795207099100D+00, & -0.8945795207341440D+00, 0.2787624047279170D+00, 0.3493121852918660D+00, & -0.3493121854667010D+00, -0.8945795206777230D+00, -0.2787624046898960D+00, & -0.8945795207888640D+00, -0.2787624046586770D+00, -0.3493121852069840D+00, & -0.3493121855510410D+00, 0.8945795206827980D+00, 0.2787624045679230D+00, & 0.8945795207852190D+00, -0.2787624046804690D+00, 0.3493121851989290D+00, & 0.3493121855496230D+00, 0.8945795206792300D+00, -0.2787624045811490D+00, & 0.8945795207818050D+00, 0.2787624045559080D+00, -0.3493121853070750D+00, & 0.2787624044379500D+00, -0.3493121855065000D+00, 0.8945795207406920D+00, & -0.2787624044432590D+00, 0.3493121854287870D+00, 0.8945795207693820D+00, & 0.5558962301794150D+00, -0.6768332117366710D+00, 0.4825724658147600D+00, & -0.6768332116815670D+00, 0.4825724660401160D+00, 0.5558962300508760D+00, & -0.5558962303148920D+00, 0.6768332115229870D+00, 0.4825724659584010D+00, & 0.4825724659102830D+00, 0.5558962301646720D+00, -0.6768332116806730D+00, & 0.6768332114576920D+00, 0.4825724660928950D+00, -0.5558962302776390D+00, & 0.4825724659029810D+00, -0.5558962303679090D+00, 0.6768332115189570D+00, & 0.6768332116355920D+00, -0.4825724660719810D+00, 0.5558962300791910D+00, & -0.4825724661505860D+00, -0.5558962302300841D+00, -0.6768332114556160D+00, & -0.6768332114382860D+00, -0.4825724663277370D+00, -0.5558962300974000D+00, & -0.4825724659723730D+00, 0.5558962302677700D+00, 0.6768332115517270D+00, & 0.5558962301926910D+00, 0.6768332115894530D+00, -0.4825724660059490D+00, & -0.5558962301943380D+00, -0.6768332114555370D+00, -0.4825724661918750D+00 & /), (/ 3, 240 /) ) return end subroutine design_order_max ( order_max ) !*****************************************************************************80 ! !! DESIGN_ORDER_MAX returns the maximum design order. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 09 June 2002 ! ! Author: ! ! John Burkardt ! ! Output: ! ! integer ORDER_MAX, the maximum design order. ! implicit none integer order_max order_max = 21 return end subroutine design_points ( order, n, x ) !*****************************************************************************80 ! !! DESIGN_POINTS returns the points for a given 3D spherical design. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 09 June 2002 ! ! Author: ! ! John Burkardt ! ! Input: ! ! integer ORDER, the desired order of the design. ! ! integer N, the number of points associated with the ! design, as returned by DESIGN_ORDER. ! ! Output: ! ! real ( kind = rk8 ) X(3,N), the points associated with the design. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) integer n integer order real ( kind = rk8 ) x(n) if ( order == 0 ) then call design_00_001_3d ( x ) else if ( order == 1 ) then call design_01_002_3d ( x ) else if ( order == 2 ) then call design_02_004_3d ( x ) else if ( order == 3 ) then call design_03_006_3d ( x ) else if ( order == 4 ) then call design_04_014_3d ( x ) else if ( order == 5 ) then call design_05_012_3d ( x ) else if ( order == 6 ) then call design_06_026_3d ( x ) else if ( order == 7 ) then call design_07_024_3d ( x ) else if ( order == 8 ) then call design_08_036_3d ( x ) else if ( order == 9 ) then call design_09_048_3d ( x ) else if ( order == 10 ) then call design_10_060_3d ( x ) else if ( order == 11 ) then call design_11_070_3d ( x ) else if ( order == 12 ) then call design_12_084_3d ( x ) else if ( order == 13 ) then call design_13_094_3d ( x ) else if ( order == 14 ) then call design_14_108_3d ( x ) else if ( order == 15 ) then call design_15_120_3d ( x ) else if ( order == 16 ) then call design_16_144_3d ( x ) else if ( order == 17 ) then call design_17_156_3d ( x ) else if ( order == 18 ) then call design_18_180_3d ( x ) else if ( order == 19 ) then call design_19_204_3d ( x ) else if ( order == 20 ) then call design_20_216_3d ( x ) else if ( order == 21 ) then call design_21_240_3d ( x ) else write ( *, '(a)' ) ' ' write ( *, '(a)' ) 'DESIGN_POINTS - Fatal error!' write ( *, '(a,i6)' ) ' There is no information for ORDER = ', order stop end if return end subroutine design_quad ( xc, r, func_3d, order, quad ) !*****************************************************************************80 ! !! DESIGN_QUAD approximates an integral with a design of a given order. ! ! Discussion: ! ! The design is defined for the unit sphere in 3D, and provides a ! guaranteed accuracy for integrals of polynomials of a given ! degree over the surface of the unit sphere. ! ! The same guarantee applies if the sphere is scaled to have ! a non-unit radius R. ! ! The same guarantee applies if the center of the sphere is ! translated by XC. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 08 June 2002 ! ! Author: ! ! John Burkardt ! ! Input: ! ! real ( kind = rk8 ) XC(3), the coordinates of the center of the sphere. ! ! real ( kind = rk8 ) R, the radius of the sphere. ! ! real ( kind = rk8 ) external FUNC_3D, evaluates the function. ! ! integer ORDER, the desired order of the design. ! ! Output: ! ! real ( kind = rk8 ) QUAD, the approximate integral of the function. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ), external :: func_3d integer i integer j integer n integer order real ( kind = rk8 ) quad real ( kind = rk8 ) r real ( kind = rk8 ), allocatable, dimension ( :, : ) :: x real ( kind = rk8 ) xc(3) call design_size ( order, n ) if ( n == 0 ) then quad = 0.0D+00 return end if allocate ( x(1:3,1:n) ) call design_points ( order, n, x ) x(1:3,1:n) = r * x(1:3,1:n) do i = 1, 3 x(i,1:n) = x(i,1:n) + xc(i) end do if ( .false. ) then write ( *, '(a)' ) ' ' write ( *, '(a)' ) 'The design points:' write ( *, '(a)' ) ' ' do j = 1, n write ( *, '(i6,3g14.6)' ) j, x(1:3,j) end do end if quad = 0.0D+00 do j = 1, n quad = quad + func_3d ( x(1:3,j) ) end do quad = quad * r**2 / real ( n, kind = rk8 ) deallocate ( x ) return end subroutine design_size ( order, n ) !*****************************************************************************80 ! !! DESIGN_SIZE returns the size of a 3D spherical design of a given order. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 09 June 2002 ! ! Author: ! ! John Burkardt ! ! Input: ! ! integer ORDER, the desired order of the design. ! ! Output: ! ! integer N, the number of points of the design. ! implicit none integer n integer order if ( order == 0 ) then n = 1 else if ( order == 1 ) then n = 2 else if ( order == 2 ) then n = 4 else if ( order == 3 ) then n = 6 else if ( order == 4 ) then n = 14 else if ( order == 5 ) then n = 12 else if ( order == 6 ) then n = 26 else if ( order == 7 ) then n = 24 else if ( order == 8 ) then n = 36 else if ( order == 9 ) then n = 48 else if ( order == 10 ) then n = 60 else if ( order == 11 ) then n = 70 else if ( order == 12 ) then n = 84 else if ( order == 13 ) then n = 94 else if ( order == 14 ) then n = 108 else if ( order == 15 ) then n = 120 else if ( order == 16 ) then n = 144 else if ( order == 17 ) then n = 156 else if ( order == 18 ) then n = 180 else if ( order == 19 ) then n = 204 else if ( order == 20 ) then n = 216 else if ( order == 21 ) then n = 240 else n = 0 end if return end subroutine get_unit ( iunit ) !*****************************************************************************80 ! !! GET_UNIT returns a free FORTRAN unit number. ! ! Discussion: ! ! A "free" FORTRAN unit number is an integer between 1 and 99 which ! is not currently associated with an I/O device. A free FORTRAN unit ! number is needed in order to open a file with the OPEN command. ! ! If IUNIT = 0, then no free FORTRAN unit could be found, although ! all 99 units were checked (except for units 5, 6 and 9, which ! are commonly reserved for console I/O). ! ! Otherwise, IUNIT is an integer between 1 and 99, representing a ! free FORTRAN unit. Note that GET_UNIT assumes that units 5 and 6 ! are special, and will never return those values. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 26 October 2008 ! ! Author: ! ! John Burkardt ! ! Output: ! ! integer IUNIT, the free unit number. ! implicit none integer i integer ios integer iunit logical lopen iunit = 0 do i = 1, 99 if ( i /= 5 .and. i /= 6 .and. i /= 9 ) then inquire ( unit = i, opened = lopen, iostat = ios ) if ( ios == 0 ) then if ( .not. lopen ) then iunit = i return end if end if end if end do return end subroutine r8mat_write ( output_filename, m, n, table ) !*****************************************************************************80 ! !! R8MAT_WRITE writes an R8MAT file. ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 31 May 2009 ! ! Author: ! ! John Burkardt ! ! Input: ! ! character ( len = * ) OUTPUT_FILENAME, the output file name. ! ! integer M, the spatial dimension. ! ! integer N, the number of points. ! ! real ( kind = rk8 ) TABLE(M,N), the table data. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) integer m integer n integer j character ( len = * ) output_filename integer output_status integer output_unit character ( len = 30 ) string real ( kind = rk8 ) table(m,n) ! ! Open the file. ! call get_unit ( output_unit ) open ( unit = output_unit, file = output_filename, & status = 'replace', iostat = output_status ) if ( output_status /= 0 ) then write ( *, '(a)' ) ' ' write ( *, '(a)' ) 'R8MAT_WRITE - Fatal error!' write ( *, '(a,i8)' ) ' Could not open the output file "' // & trim ( output_filename ) // '" on unit ', output_unit output_unit = -1 stop end if ! ! Create a format string. ! ! For less precision in the output file, try: ! ! '(', m, 'g', 14, '.', 6, ')' ! if ( 0 < m .and. 0 < n ) then write ( string, '(a1,i8,a1,i8,a1,i8,a1)' ) '(', m, 'g', 24, '.', 16, ')' ! ! Write the data. ! do j = 1, n write ( output_unit, string ) table(1:m,j) end do end if ! ! Close the file. ! close ( unit = output_unit ) return end subroutine sphere_area_3d ( r, area ) !*****************************************************************************80 ! !! SPHERE_AREA_3D computes the surface area of a sphere in 3D. ! ! Discussion: ! ! A sphere in 3D satisfies the equation: ! ! ( X - XC )^2 + ( Y - YC )^2 + ( Z - ZC )^2 = R^2 ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 26 January 1999 ! ! Author: ! ! John Burkardt ! ! Input: ! ! real ( kind = rk8 ) R, the radius of the sphere. ! ! Output: ! ! real ( kind = rk8 ) AREA, the area of the sphere. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) real ( kind = rk8 ) area real ( kind = rk8 ), parameter :: pi = 3.141592653589793D+00 real ( kind = rk8 ) r area = 4.0D+00 * pi * r * r return end subroutine sphere_monomial_int_3d ( r, e, integral ) !*****************************************************************************80 ! !! SPHERE_MONOMIAL_INT_3D integrates a monomial on a sphere in 3D. ! ! Discussion: ! ! The sphere may have nonunit radius, but it must be centered at 0. ! ! The integration region is ! ! X^2 + Y^2 + Z^2 = R^2. ! ! The monomial is F(X,Y,Z) = X^E(1) * Y^E(2) * Z^E(3). ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 24 June 2002 ! ! Author: ! ! John Burkardt ! ! Reference: ! ! Philip Davis, Philip Rabinowitz, ! Methods of Numerical Integration, ! Second Edition, ! Academic Press, 1984, page 263. ! ! Input: ! ! real ( kind = rk8 ) R, the radius of the sphere. ! ! integer E(3), the exponents of X, Y and Z in the ! monomial. Each exponent must be nonnegative. ! ! Output: ! ! real ( kind = rk8 ) INTEGRAL, the integral. ! implicit none integer, parameter :: rk8 = kind ( 1.0D+00 ) integer, parameter :: n = 3 integer e(n) integer i real ( kind = rk8 ) integral real ( kind = rk8 ), parameter :: pi = 3.141592653589793D+00 real ( kind = rk8 ) r if ( any ( e(1:n) < 0 ) ) then integral = - huge ( integral ) write ( *, '(a)' ) ' ' write ( *, '(a)' ) 'SPHERE_MONOMIAL_INT_3D - Fatal error!' write ( *, '(a)' ) ' All exponents must be nonnegative.' stop end if if ( all ( e(1:n) == 0 ) ) then integral = 2.0D+00 * sqrt ( pi**n ) & / gamma ( 0.5D+00 * real ( n, kind = rk8 ) ) else if ( any ( mod ( e(1:n), 2 ) == 1 ) ) then integral = 0.0D+00 else integral = 2.0D+00 do i = 1, n integral = integral * gamma ( 0.5D+00 * real ( e(i) + 1, kind = rk8 ) ) end do integral = integral & / gamma ( 0.5D+00 * ( real ( sum ( e(1:n) + 1 ), kind = rk8 ) ) ) end if integral = integral * r**( sum ( e(1:n) ) + 2 ) return end subroutine timestamp ( ) !*****************************************************************************80 ! !! TIMESTAMP prints the current YMDHMS date as a time stamp. ! ! Example: ! ! 31 May 2001 9:45:54.872 AM ! ! Licensing: ! ! This code is distributed under the MIT license. ! ! Modified: ! ! 06 August 2005 ! ! Author: ! ! John Burkardt ! implicit none character ( len = 8 ) ampm integer d integer h integer m integer mm character ( len = 9 ), parameter, dimension(12) :: month = (/ & 'January ', 'February ', 'March ', 'April ', & 'May ', 'June ', 'July ', 'August ', & 'September', 'October ', 'November ', 'December ' /) integer n integer s integer values(8) integer y call date_and_time ( values = values ) y = values(1) m = values(2) d = values(3) h = values(5) n = values(6) s = values(7) mm = values(8) if ( h < 12 ) then ampm = 'AM' else if ( h == 12 ) then if ( n == 0 .and. s == 0 ) then ampm = 'Noon' else ampm = 'PM' end if else h = h - 12 if ( h < 12 ) then ampm = 'PM' else if ( h == 12 ) then if ( n == 0 .and. s == 0 ) then ampm = 'Midnight' else ampm = 'AM' end if end if end if write ( *, '(i2,1x,a,1x,i4,2x,i2,a1,i2.2,a1,i2.2,a1,i3.3,1x,a)' ) & d, trim ( month(m) ), y, h, ':', n, ':', s, '.', mm, trim ( ampm ) return end