subroutine design_00_001_3d ( x ) !*****************************************************************************80 ! !! DESIGN_00_001_3D returns a 0-design of 1 point in 3D. ! ! 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 = 8 ) X(3,1), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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 = 8 ) X(3,2), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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 = 8 ) X(3,4), the coordinates of the points. ! implicit none real ( kind = 8 ), parameter :: a = 0.577350269189625763D+00 real ( kind = 8 ) 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. ! ! 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 = 8 ) X(3,6), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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 = 8 ) X(3,14), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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 = 8 ) X(3,12), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,26), the coordinates of the points. ! implicit none real ( kind = 8 ) x(3,26) x(1:3,1:26) = reshape ( (/ & -0.6156289423339210D+00, 0.1008195633550640D+00, 0.7815602478412730D+00, & 0.3951267761987170D+00, 0.4827405208135130D+00, 0.7815602473869220D+00, & 0.2205021660506020D+00, -0.5835600880219350D+00, 0.7815602461965580D+00, & -0.4164742585347590D+00, 0.6658869344429950D+00, 0.6189860923445970D+00, & -0.3684378761516490D+00, -0.6936207556288800D+00, 0.6189860893248600D+00, & 0.7849121331270050D+00, 0.2773382046045800D-01, 0.6189860890784880D+00, & -0.9753444939429849D+00, 0.2176679667586930D+00, 0.3638371039680900D-01, & 0.6761782314481880D+00, 0.7358391298492950D+00, 0.3638370924238100D-01, & 0.2991662611656320D+00, -0.9535070916939770D+00, 0.3638370884709700D-01, & 0.7506752044088399D+00, -0.5662857931564620D+00, 0.3403044783056840D+00, & 0.1150802794034220D+00, 0.9332466938773990D+00, 0.3403044778714140D+00, & -0.8657554869972310D+00, -0.3669608990989080D+00, 0.3403044743559930D+00, & 0.0D+00, -0.0D+00, 1.0000000000000000D+00, & -0.0D+00, -0.0D+00, -1.0000000000000000D+00, & 0.8657554869972310D+00, -0.3669608990989070D+00, -0.3403044743559920D+00, & -0.1150802794034210D+00, 0.9332466938773980D+00, -0.3403044778714150D+00, & -0.7506752044088399D+00, -0.5662857931564620D+00, -0.3403044783056840D+00, & -0.2991662611656320D+00, -0.9535070916939770D+00, -0.3638370884709700D-01, & -0.6761782314481870D+00, 0.7358391298492950D+00, -0.3638370924238000D-01, & 0.9753444939429849D+00, 0.2176679667586940D+00, -0.3638371039681000D-01, & -0.7849121331270040D+00, 0.2773382046045900D-01, -0.6189860890784890D+00, & 0.3684378761516490D+00, -0.6936207556288810D+00, -0.6189860893248590D+00, & 0.4164742585347590D+00, 0.6658869344429950D+00, -0.6189860923445970D+00, & -0.2205021660506030D+00, -0.5835600880219350D+00, -0.7815602461965580D+00, & -0.3951267761987180D+00, 0.4827405208135130D+00, -0.7815602473869220D+00, & 0.6156289423339210D+00, 0.1008195633550640D+00, -0.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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,24), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,36), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,48), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,60), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,70), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,84), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,94), the coordinates of the points. ! implicit none real ( kind = 8 ) 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. ! ! Discussion: ! ! The contortions necessary to get this routine to work arise because ! of the idiotic limit of 99 continuation statements, AND the ! inability of the RESHAPE command to deal with situations involving ! an array section. As far as I can tell, you must store the results ! of a RESHAPE command in an array that is exactly the dimensions ! given in the command. ! ! 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. ! ! Parameters: ! ! Output, real ( kind = 8 ) X(3,108), the coordinates of the points. ! implicit none real ( kind = 8 ) x(3,108) real ( kind = 8 ) temp(3,50) temp(1:3,1:50) = 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 & /), (/ 3, 50 /) ) x(1:3,1:50) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & 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 & /), (/ 3, 50 /) ) x(1:3,51:100) = temp(1:3,1:50) temp(1:3,1:8) = reshape ( (/ & 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, 8 /) ) x(1:3,101:108) = temp(1:3,1:8) return end subroutine design_15_120_3d ( x ) !*****************************************************************************80 ! !! DESIGN_15_120_3D returns a 15-design of 120 points in 3D. ! ! Discussion: ! ! The contortions necessary to get this routine to work arise because ! of the idiotic limit of 99 continuation statements, AND the ! inability of the RESHAPE command to deal with situations involving ! an array section. As far as I can tell, you must store the results ! of a RESHAPE command in an array that is exactly the dimensions ! given in the command. ! ! Modified: ! ! 09 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 = 8 ) X(3,120), the coordinates of the points. ! implicit none real ( kind = 8 ) x(3,120) real ( kind = 8 ) temp(3,50) temp(1:3,1:50) = reshape ( (/ & 0.8544032798674690D+00, -0.5053541340072060D+00, 0.1208810762424740D+00, & -0.5054349175556900D+00, 0.1208162198059960D+00, 0.8543646675438200D+00, & -0.8543867766655620D+00, 0.5053247652039460D+00, 0.1211202606115420D+00, & 0.1208333586366210D+00, 0.8543977898340150D+00, -0.5053748273977881D+00, & 0.5053979097545750D+00, 0.1211845075248970D+00, -0.8543344005432850D+00, & 0.1211678917817770D+00, -0.8543595921698920D+00, 0.5053593070959080D+00, & 0.5055502439906060D+00, -0.1210290998832230D+00, 0.8542662979320300D+00, & -0.1209012300582570D+00, -0.8542712308899000D+00, -0.5055725038451520D+00, & -0.5055127430804750D+00, -0.1209718018932920D+00, -0.8542966052431350D+00, & -0.1211002217099370D+00, 0.8542333106705450D+00, 0.5055889508708080D+00, & 0.8542280860180770D+00, 0.5056228637404400D+00, -0.1209954409091880D+00, & -0.8542442519150010D+00, -0.5055933391780419D+00, -0.1210046835828190D+00, & -0.2649878987783750D+00, 0.8838136985753620D+00, -0.3855577255244170D+00, & 0.8838495436614180D+00, -0.3855147723237870D+00, -0.2649308296322680D+00, & 0.2649355775867900D+00, -0.8837175312641120D+00, -0.3858140286008490D+00, & -0.3855859698282140D+00, -0.2649918821321200D+00, 0.8838001823238730D+00, & -0.8837292209025740D+00, -0.3858572041282210D+00, 0.2648336877084960D+00, & -0.3857855708943260D+00, 0.2648743412798280D+00, -0.8837483106752260D+00, & -0.8838853342582800D+00, 0.3855799306321350D+00, -0.2647165143636620D+00, & 0.3857055651888330D+00, 0.2647713629504000D+00, 0.8838140881111540D+00, & 0.8837647796773470D+00, 0.3857919932898060D+00, 0.2648100302392800D+00, & 0.3856671594356380D+00, -0.2646531289490930D+00, -0.8838662588142510D+00, & -0.2646569340063700D+00, -0.8838970336929440D+00, 0.3855940107304080D+00, & 0.2647089806257060D+00, 0.8838007658213300D+00, 0.3857789028831530D+00, & -0.9733521648930310D+00, 0.2280262392539510D+00, 0.2428162494047400D-01, & 0.2281123458869260D+00, 0.2434454914745700D-01, -0.9733304169606380D+00, & 0.9733554428092021D+00, -0.2280321804631890D+00, 0.2409370513031500D-01, & 0.2435678809852000D-01, -0.9733428333593150D+00, 0.2280580531829240D+00, & -0.2281196464116130D+00, 0.2403159282358800D-01, 0.9733364831687970D+00, & 0.2401983627102700D-01, 0.9733503190618910D+00, -0.2280618421560860D+00, & -0.2282433008441410D+00, -0.2410134029444700D-01, -0.9733057695379100D+00, & -0.2432552782893800D-01, 0.9732858593312290D+00, 0.2283044124008880D+00, & 0.2282511002315090D+00, -0.2427391182448100D-01, 0.9732996519304030D+00, & -0.2405063099923100D-01, -0.9732933319509690D+00, -0.2283016800821220D+00, & -0.9732840684188760D+00, -0.2283305892567410D+00, -0.2415086275208100D-01, & 0.9732807668974420D+00, 0.2283367922856280D+00, -0.2422515379181500D-01, & 0.1764941643605970D+00, -0.6439159617576869D+00, 0.7444609084030720D+00, & -0.6439554747347650D+00, 0.7444298798961950D+00, 0.1764808785020650D+00, & -0.1763945457957580D+00, 0.6437302463633759D+00, 0.7446451061616230D+00, & 0.7444821769849420D+00, 0.1765407684656810D+00, -0.6438785950948430D+00, & 0.6437221735449969D+00, 0.7446755654083020D+00, -0.1762953929356370D+00, & 0.7446232602985000D+00, -0.1763087546453890D+00, 0.6437790174103420D+00, & 0.6439796469919020D+00, -0.7444731281720970D+00, 0.1762100328864330D+00, & -0.7445688151995640D+00, -0.1763819450908200D+00, -0.6438219387985280D+00, & -0.6437453004186960D+00, -0.7446323171490740D+00, -0.1763935952523330D+00, & -0.7445366073086630D+00, 0.1761491742553250D+00, 0.6439229059339890D+00, & 0.1761957070235800D+00, 0.6439711185624940D+00, -0.7444838959197370D+00, & -0.1762957903461830D+00, -0.6437847990685870D+00, -0.7446213311438470D+00, & -0.5382680916693570D+00, -0.7144430972079280D+00, 0.4470330215348440D+00, & -0.7144085713422340D+00, 0.4470401217976180D+00, -0.5383080184206051D+00 & /), (/ 3, 50 /) ) x(1:3,1:50) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & 0.5383282112194380D+00, 0.7143316685506840D+00, 0.4471386857686070D+00, & 0.4471201978095990D+00, -0.5382526202846060D+00, -0.7144001997952280D+00, & 0.7142693383899240D+00, 0.4471311678997200D+00, 0.5384171532637620D+00, & 0.4470507600303290D+00, 0.5383906456527831D+00, 0.7143396465476940D+00, & 0.7143353780369111D+00, -0.4469080442623900D+00, -0.5385147794243240D+00, & -0.4472169730606700D+00, 0.5384307414350790D+00, -0.7142053736035060D+00, & -0.7141970542020020D+00, -0.4472644791667070D+00, 0.5384023156013700D+00, & -0.4469551475186520D+00, -0.5385686162875321D+00, 0.7142653160112939D+00, & -0.5385531358716100D+00, 0.7142734853459000D+00, -0.4469607454516280D+00, & 0.5384924234925470D+00, -0.7141627498790130D+00, -0.4472107741617700D+00, & -0.8542625591715191D+00, -0.1211967864813340D+00, 0.5055163884033080D+00, & -0.1211352259177740D+00, 0.5055628408198980D+00, -0.8542438006939030D+00, & 0.8543309414267219D+00, 0.1210719374548010D+00, 0.5054307355927919D+00, & 0.5056009275709700D+00, -0.8542205006653180D+00, -0.1211405724044100D+00, & 0.1209804838657270D+00, 0.5053853800887540D+00, 0.8543707275627840D+00, & 0.5053474173253790D+00, 0.8543800060355530D+00, 0.1210735028371490D+00, & 0.1210201874725290D+00, -0.5053486406391460D+00, -0.8543868360574620D+00, & -0.5056166845687340D+00, 0.8542509316027460D+00, -0.1208598946779390D+00, & -0.1208657310279020D+00, -0.5055984823689040D+00, 0.8542608791752970D+00, & -0.5053307603087600D+00, -0.8544103980209300D+00, 0.1209284682756170D+00, & -0.8543680863966890D+00, 0.1209284250962590D+00, -0.5054023040614261D+00, & 0.8543001861498960D+00, -0.1208044489183890D+00, -0.5055467110621700D+00, & 0.7444632973046910D+00, 0.6439458791655150D+00, 0.1763748952430030D+00, & 0.6438740925860960D+00, 0.1763526546392770D+00, 0.7445306535651250D+00, & -0.7444394383517560D+00, -0.6439896739556450D+00, 0.1763156897868830D+00, & 0.1762725384812350D+00, 0.7444994067041240D+00, 0.6439321591554410D+00, & -0.6439272279941870D+00, 0.1763373486135990D+00, -0.7444883239726789D+00, & 0.1764171560921340D+00, -0.7444453593527010D+00, -0.6439550403363510D+00, & -0.6437732764367420D+00, -0.1765375814946370D+00, 0.7445740063055700D+00, & -0.1761862890755310D+00, -0.7446599355387570D+00, 0.6437701235264091D+00, & 0.6438276103095120D+00, -0.1761530561791470D+00, -0.7446180960746860D+00, & -0.1765038577432470D+00, 0.7446064263886381D+00, -0.6437450255982510D+00, & 0.7446425637975280D+00, -0.6437113057563880D+00, -0.1764743806405070D+00, & -0.7446660895457380D+00, 0.6437550517808310D+00, -0.1762153466282650D+00, & -0.2283360715319860D+00, 0.9732852424848990D+00, -0.2405151135480500D-01, & 0.9733064515717200D+00, -0.2400727784605000D-01, -0.2282503054527860D+00, & 0.2283331981063670D+00, -0.9732791056699990D+00, -0.2432556492099900D-01, & -0.2408702223521400D-01, -0.2283063346170040D+00, 0.9732913402129850D+00, & -0.9732989361424010D+00, -0.2436839438764100D-01, 0.2282440848281470D+00, & -0.2428855197935900D-01, 0.2282995532397890D+00, -0.9732879225764900D+00, & -0.9733371742315819D+00, 0.2406421970882400D-01, -0.2281132582483620D+00, & 0.2421829905214100D-01, 0.2280639631571710D+00, 0.9733449042862790D+00, & 0.9733297074180570D+00, 0.2431156319331100D-01, 0.2281188912663230D+00, & 0.2415744670603300D-01, -0.2280561009745210D+00, -0.9733482586294300D+00, & -0.2280265869033880D+00, -0.9733570472797149D+00, 0.2408182252460000D-01, & 0.2280302830467680D+00, 0.9733509134612350D+00, 0.2429381151218600D-01, & 0.7141882216414780D+00, 0.5385777705783590D+00, -0.4470672981861090D+00, & 0.5385274786436980D+00, -0.4470918044321980D+00, 0.7142108044234720D+00, & -0.7142487687796120D+00, -0.5384672206579190D+00, -0.4471037335716890D+00, & -0.4471613350315650D+00, 0.7141657502813270D+00, 0.5385294992643370D+00 & /), (/ 3, 50 /) ) x(1:3,51:100) = temp(1:3,1:50) temp(1:3,1:20) = reshape ( (/ & -0.5383895843123190D+00, -0.4470804117879510D+00, -0.7143218888565050D+00, & -0.4470109694572090D+00, -0.7143051130884170D+00, -0.5384694964440140D+00, & -0.5384320154150550D+00, 0.4469136813525910D+00, 0.7143942372359660D+00, & 0.4472327482991630D+00, -0.7143004980317570D+00, 0.5382914334822370D+00, & 0.5382930399261510D+00, 0.4472578784152030D+00, -0.7142835524934020D+00, & 0.4469384509810380D+00, 0.7144394343055710D+00, -0.5383514797451630D+00, & 0.7144166137892710D+00, -0.5383541595409150D+00, 0.4469717002751600D+00, & -0.7143571305190370D+00, 0.5382432243381430D+00, 0.4472003147703360D+00, & -0.8838226247963700D+00, 0.2646362039654440D+00, 0.3857787545327180D+00, & 0.2647035776402780D+00, 0.3858370666585500D+00, -0.8837769933511300D+00, & 0.8838743905237270D+00, -0.2647316811045790D+00, 0.3855946042099840D+00, & 0.3858426620308480D+00, -0.8837798638376970D+00, 0.2646858372334470D+00, & -0.2648227787783800D+00, 0.3855352957965560D+00, 0.8838729725108470D+00, & 0.3855302957630180D+00, 0.8839005568372480D+00, -0.2647379773883650D+00, & -0.2648230436494790D+00, -0.3856162950490920D+00, -0.8838375577813140D+00, & -0.3858066593719060D+00, 0.8837138197040250D+00, 0.2649586881919760D+00, & 0.2649414017790060D+00, -0.3857554525672310D+00, 0.8837413560754220D+00, & -0.3855652667157560D+00, -0.8838346669397210D+00, -0.2649069772919530D+00, & -0.8837919506738630D+00, -0.2649137504830990D+00, -0.3856585182892600D+00, & 0.8837398529718770D+00, 0.2650088880255400D+00, -0.3857125374378070D+00 & /), (/ 3, 20 /) ) x(1:3,101:120) = temp(1:3,1:20) return end subroutine design_16_144_3d ( x ) !*****************************************************************************80 ! !! DESIGN_16_144_3D returns a 16-design of 144 points in 3D. ! ! Discussion: ! ! The contortions necessary to get this routine to work arise because ! of the idiotic limit of 99 continuation statements, AND the ! inability of the RESHAPE command to deal with situations involving ! an array section. As far as I can tell, you must store the results ! of a RESHAPE command in an array that is exactly the dimensions ! given in the command. ! ! Modified: ! ! 09 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 = 8 ) X(3,144), the coordinates of the points. ! implicit none real ( kind = 8 ) x(3,144) real ( kind = 8 ) temp(3,50) temp(1:3,1:50) = 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 & /), (/ 3, 50 /) ) x(1:3,1:50) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & 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 & /), (/ 3, 50 /) ) x(1:3,51:100) = temp(1:3,1:50) temp(1:3,1:44) = reshape ( (/ & 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, 44 /) ) x(1:3,101:144) = temp(1:3,1:44) return end subroutine design_17_156_3d ( x ) !*****************************************************************************80 ! !! DESIGN_17_156_3D returns a 17-design of 156 points in 3D. ! ! Discussion: ! ! The contortions necessary to get this routine to work arise because ! of the idiotic limit of 99 continuation statements, AND the ! inability of the RESHAPE command to deal with situations involving ! an array section. As far as I can tell, you must store the results ! of a RESHAPE command in an array that is exactly the dimensions ! given in the command. ! ! Modified: ! ! 09 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 = 8 ) X(3,156), the coordinates of the points. ! implicit none real ( kind = 8 ) x(3,156) real ( kind = 8 ) temp(3,50) temp(1:3,1:50) = 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 & /), (/ 3, 50 /) ) x(1:3,1:50) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & -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 & /), (/ 3, 50 /) ) x(1:3,51:100) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & -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 & /), (/ 3, 50 /) ) x(1:3,101:150) = temp(1:3,1:50) temp(1:3,1:6) = reshape ( (/ & -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, 6 /) ) x(1:3,151:156) = temp(1:3,1:6) return end subroutine design_18_180_3d ( x ) !*****************************************************************************80 ! !! DESIGN_18_180_3D returns an 18-design of 180 points in 3D. ! ! Discussion: ! ! The contortions necessary to get this routine to work arise because ! of the idiotic limit of 99 continuation statements, AND the ! inability of the RESHAPE command to deal with situations involving ! an array section. As far as I can tell, you must store the results ! of a RESHAPE command in an array that is exactly the dimensions ! given in the command. ! ! Modified: ! ! 09 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 = 8 ) X(3,180), the coordinates of the points. ! implicit none real ( kind = 8 ) x(3,180) real ( kind = 8 ) temp(3,50) temp(1:3,1:50) = 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 & /), (/ 3, 50 /) ) x(1:3,1:50) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & -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 & /), (/ 3, 50 /) ) x(1:3,51:100) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & -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 & /), (/ 3, 50 /) ) x(1:3,101:150) = temp(1:3,1:50) temp(1:3,1:30) = reshape ( (/ & 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, 30 /) ) x(1:3,151:180) = temp(1:3,1:30) return end subroutine design_19_204_3d ( x ) !*****************************************************************************80 ! !! DESIGN_19_204_3D returns a 19-design of 204 points in 3D. ! ! Discussion: ! ! The contortions necessary to get this routine to work arise because ! of the idiotic limit of 99 continuation statements, AND the ! inability of the RESHAPE command to deal with situations involving ! an array section. As far as I can tell, you must store the results ! of a RESHAPE command in an array that is exactly the dimensions ! given in the command. ! ! Modified: ! ! 09 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 = 8 ) X(3,204), the coordinates of the points. ! implicit none real ( kind = 8 ) x(3,204) real ( kind = 8 ) temp(3,50) temp(1:3,1:50) = 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 & /), (/ 3, 50 /) ) x(1:3,1:50) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & -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 & /), (/ 3, 50 /) ) x(1:3,51:100) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & -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 & /), (/ 3, 50 /) ) x(1:3,101:150) = temp(1:3,1:50) temp(1:3,1:50) = reshape ( (/ & 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 & /), (/ 3, 50 /) ) x(1:3,151:200) = temp(1:3,1:50) temp(1:3,1:4) = reshape ( (/ & -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, 4 /) ) x(1:3,201:204) = temp(1:3,1:4) return end subroutine design_20_216_3d ( x ) !*****************************************************************************80 ! !! DESIGN_20_216_3D returns a 20-design of 216 points in 3D. ! ! Discussion: ! ! The contortions necessary to get this routine to work arise because ! of the idiotic limit of 99 continuation statements, AND the ! inability of the RESHAPE command to deal with situations involving ! an array section. As far as I can tell, you must store the results ! of a RES