"複数の状態が想定されるステップ、球面上に多数のセルが想定されるステップ、それぞれのセルの初期の状態が決定されるステップ、このセルの状態の変更の要否が、このセルの近傍に位置する複数のセルの状態に基づいて判定されるステップ、上記判定に基づき、このセルの状態の更新がなされるステップ及びこのセルに、その状態に応じてランド又はリセスが割り当てられるステップを含む、ゴルフボール表面の凹凸パターン設計方法であって、上記判定及び更新が、セルオートマトン法の反応拡散モデルによってなされる設計方法。上記初期の状態が無作為に決定される請求項1に記載の設計方法。上記判定及び更新が3回以上繰り返された後に、それぞれのセルに、その状態に応じてランド又はリセスが割り当てられる請求項1又は2に記載の設計方法。上記セルの数が5000個以上100000個以下である請求項1から3のいずれかに記載の設計方法。下記数式1で算出される値Eに基づいて、状態の変更の要否が判定される請求項1から4のいずれかに記載の設計方法。E=W1NR1W2NR1-R21この数式1において、W1は第一濃度であり、NR1は第一円に含まれこの第一円の中心に位置しない特定の状態にあるセルの数であり、W2は第二濃度であり、NR1-R2は第二円に含まれ第一円に含まれない特定の状態にあるセルの数である。第一濃度は正であり、第二濃度は負である。第一円は、当該セルを中心とする、指数半径がR1である円である。第二円は、当該セルを中心とする、指数半径がR2である円である。半径R2は、半径R1よりも大きい。それぞれのセルの初期の状態が分化又は未分化であり、下記数式1で算出される値Eが正であるとき、セルの状態が分化であればこの状態が維持され、セルの状態が未分化であればこの状態が分化に変更され、この値Eがゼロであるとき、セルの状態が維持され、この値Eが負であるとき、セルの状態が分化であればこの状態が未分化に変更され、セルの状態が未分化であればこの状態が維持される請求項5に記載の設計方法。E=W1NR1W2NR1-R21この数式1において、W1は第一濃度であり、NR1は第一円に含まれこの第一円の中心に位置しない分化のセルの合計数であり、W2は第二濃度であり、NR1-R2は第二円に含まれ第一円に含まれない分化のセルの合計数である。第一濃度は正であり、第二濃度は負である。第一円は、当該セルを中心とする、指数半径がR1である円である。第二円は、当該セルを中心とする、指数半径がR2である円である。半径R2は、半径R1よりも大きい。上記第一濃度W1が0.80以上1.20以下である請求項5又は6に記載の設計方法。上記第二濃度W2が-0.70以上-0.50以下である請求項5から7のいずれかに記載の設計方法。上記指数半径R1が2.20以上2.60以下である請求項5から8のいずれかに記載の設計方法。上記指数半径R2が3.0以上5.0以下である請求項5から9のいずれかに記載の設計方法。その表面に凹凸パターンを備えており、この凹凸パターンが請求項1から10のいずれかに記載の方法で設計されたゴルフボール。下記ステップ1から17によって得られる差dRの絶対値が2.5mm以下である請求項11に記載のゴルフボール。1ゴルフボールの両極を結ぶ線が、第一回転軸に想定されるステップ2ゴルフボールの仮想球の表面に存在し、かつ上記第一回転軸と直交する大円が想定されるステップ3ゴルフボールの仮想球の表面に存在し、上記第一回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ4これらの小円により上記仮想球が区画され、仮想球の表面のうちこれら小円に挟まれた領域が特定されるステップ5上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ6それぞれの点から上記第一回転軸に下ろした垂線の長さL1が算出されるステップ7軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ8回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Rhが算出されるステップ9上記ステップ1で想定された第一回転軸に直交する第二回転軸が想定されるステップ10ゴルフボールの仮想球の表面に存在し、かつ上記第二回転軸と直交する大円が想定されるステップ11ゴルフボールの仮想球の表面に存在し、上記第二回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ12これらの小円により上記仮想球が区画され、仮想球の表面のうちこれら小円に挟まれた領域が特定されるステップ13上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ14それぞれの点から上記第二回転軸に下ろした垂線の長さL1が算出されるステップ15軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ16回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Roが算出されるステップ17上記変動幅Rh及びRoの差dRが算出されるステップ上記差dRの絶対値が1.0mm以下である請求項12に記載のゴルフボール。特開昭61-284264号公報に開示されたゴルフボールでは、ディンプルパターンに起因する不都合が、容積の相違によって解消されている。ディンプルパターンに起因する不都合が、ディンプルパターン自体の工夫で解消されているわけではない。このゴルフボールでは、ディンプルパターンが本来備えるポテンシャルが犠牲にされる。このゴルフボールの飛距離は、十分ではない。赤道の近傍の特異性、並びにこれに起因する不十分な対称性に関し、原因究明の努力がなされている。しかし、未だ原因は明らかでなく、改善のための普遍的な理論も確立されていない。本発明の目的は、空力的対称性に優れたゴルフボールの提供にある。本発明者等は、ゴルフボール表面の凹凸パターンに工夫を施すことで、空力的対称性に優れたゴルフボールが得られることを見出した。本発明に係る凹凸パターン設計方法は、複数の状態が想定されるステップ、球面上に多数のセルが想定されるステップ、それぞれのセルの初期の状態が決定されるステップ、このセルの状態の変更の要否が、このセルの近傍に位置する複数のセルの状態に基づいて判定されるステップ、上記判定に基づき、このセルの状態の更新がなされるステップ及びこのセルに、その状態に応じてランド又はリセスが割り当てられるステップを含む。好ましくは、初期の状態は、無作為に決定される。好ましくは、判定及び更新が3回以上繰り返された後に、それぞれのセルに、その状態に応じてランド又はリセスが割り当てられる。好ましくは、判定及び更新は、セルオートマトン法によってなされる。好ましくは、判定及び更新は、セルオートマトン法の反応拡散モデルによってなされる。好ましくは、セルの数は5000個以上100000個以下である。好ましくは、状態の変更の要否は、下記数式1で算出される値Eに基づいて判定される。E=W1NR1W2NR1-R21この数式1において、W1は第一濃度であり、NR1は第一円に含まれこの第一円の中心に位置しない特定の状態のセルの数であり、W2は第二濃度であり、NR1-R2は第二円に含まれ第一円に含まれない特定の状態のセルの数である。第一濃度は正であり、第二濃度は負である。第一円は、当該セルを中心とする、指数半径がR1である円である。第二円は、当該セルを中心とする、指数半径がR2である円である。半径R2は、半径R1よりも大きい。好ましくは、それぞれのセルの初期の状態は、分化又は未分化である。下記数式1で算出される値Eが正であるとき、セルの状態が分化であればこの状態が維持され、セルの状態が未分化であればこの状態が分化に変更される。この値Eがゼロであるとき、セルの状態が維持される。この値Eが負であるとき、セルの状態が分化であればこの状態が未分化に変更され、セルの状態が未分化であればこの状態が維持される。E=W1NR1W2NR1-R21この数式1において、W1は第一濃度であり、NR1は第一円に含まれこの第一円の中心に位置しない分化のセルの合計数であり、W2は第二濃度であり、NR1-R2は第二円に含まれ第一円に含まれない分化のセルの合計数である。第一濃度は正であり、第二濃度は負である。第一円は、当該セルを中心とする、指数半径がR1である円である。第二円は、当該セルを中心とする、指数半径がR2である円である。半径R2は、半径R1よりも大きい。好ましくは、第一濃度W1は、0.80以上1.20以下である。好ましくは、第二濃度W2は、-0.70以上-0.50以下である。好ましくは、指数半径R1は、2.20以上2.60以下である。好ましくは、指数半径R2は、3.0以上5.0以下である。本発明に係るゴルフボールは、その表面に凹凸パターンを備える。この凹凸パターンは、前述の方法によって設計される。好ましくは、このゴルフボールの差dRの絶対値は、2.5mm以下である。差dRは、下記ステップ1から17によって得られる。1ゴルフボールの両極を結ぶ線が、第一回転軸に想定されるステップ2ゴルフボールの仮想球の表面に存在し、かつ上記第一回転軸と直交する大円が想定されるステップ3ゴルフボールの仮想球の表面に存在し、上記第一回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ4これらの小円により上記仮想球が区画され、仮想球の表面のうちこれら小円に挟まれた領域が特定されるステップ5上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ6それぞれの点から上記第一回転軸に下ろした垂線の長さL1が算出されるステップ7軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ8回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Rhが算出されるステップ9上記ステップ1で想定された第一回転軸に直交する第二回転軸が想定されるステップ10ゴルフボールの仮想球の表面に存在し、かつ上記第二回転軸と直交する大円が想定されるステップ11ゴルフボールの仮想球の表面に存在し、上記第二回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ12これらの小円により上記仮想球が区画され、仮想球の表面のうちこれら小円に挟まれた領域が特定されるステップ13上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ14それぞれの点から上記第二回転軸に下ろした垂線の長さL1が算出されるステップ15軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ16回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Roが算出されるステップ17上記変動幅Rh及びRoの差dRが算出されるステップ好ましくは、差dRの絶対値は、1.0mm以下である。本発明に係るゴルフボールの表面には、リセスによって多数のクレーターが形成される。これらクレーターは、乱流剥離を促進する。これらクレーターは、ゴルフボールの飛行性能に寄与する。これらクレーターはランダムに配置されるので、このパターンは方向性を有さない。このゴルフボールは、空力的対称性に優れる。以下、適宜図面が参照されつつ、好ましい実施形態に基づいて本発明が詳細に説明される。図1は、本発明の一実施形態に係るゴルフボール2が示された模式的断面図である。このゴルフボール2は、球状のコア4と、カバー6とを備えている。カバー6の表面には、多数のクレーター8が形成されている。ゴルフボール2の表面のうちクレーター8以外の部分は、ランド10である。このゴルフボール2は、カバー6の外側にペイント層及びマーク層を備えているが、これらの層の図示は省略されている。コア4とカバー6との間に、中間層が設けられてもよい。このゴルフボール2の直径は、40mm以上45mm以下である。米国ゴルフ協会USGAの規格が満たされるとの観点から、直径は42.67mm以上がより好ましい。空気抵抗抑制の観点から、直径は44mm以下がより好ましく、42.80mm以下が特に好ましい。このゴルフボール2の質量は、40g以上50g以下である。大きな慣性が得られるとの観点から、質量は44g以上がより好ましく、45.00g以上が特に好ましい。USGAの規格が満たされるとの観点から、質量は45.93g以下がより好ましい。コア4は、ゴム組成物が架橋されることによって形成されている。ゴム組成物の基材ゴムとしては、ポリブタジエン、ポリイソプレン、スチレン-ブタジエン共重合体、エチレン-プロピレン-ジエン共重合体及び天然ゴムが例示される。2種以上のゴムが併用されてもよい。反発性能の観点からポリブタジエンが好ましく、特にハイシスポリブタジエンが好ましい。コア4の架橋には、共架橋剤が用いられうる。反発性能の観点から好ましい共架橋剤は、アクリル酸亜鉛、アクリル酸マグネシウム、メタクリル酸亜鉛及びメタクリル酸マグネシウムである。ゴム組成物には、共架橋剤と共に有機過酸化物が配合されるのが好ましい。好適な有機過酸化物としては、ジクミルパーオキサイド、1,1-ビスt-ブチルパーオキシ-3,3,5-トリメチルシクロヘキサン、2,5-ジメチル-2,5-ジt-ブチルパーオキシヘキサン及びジ-t-ブチルパーオキサイドが挙げられる。コア4のゴム組成物には、硫黄化合物、充填剤、老化防止剤、着色剤、可塑剤、分散剤等の各種添加剤が、必要に応じて適量配合される。ゴム組成物に、架橋ゴム粉末又は合成樹脂粉末が配合されてもよい。コア4の直径は30.0mm以上、特には38.0mm以上である。コア4の直径は42.0mm以下、特には41.5mm以下である。コア4が2以上の層から構成されてもよい。カバー6に好適なポリマーは、アイオノマー樹脂である。好ましいアイオノマー樹脂としては、α-オレフィンと炭素数が3以上8以下のα,β-不飽和カルボン酸との二元共重合体が挙げられる。好ましい他のアイオノマー樹脂としては、α-オレフィンと炭素数が3以上8以下のα,β-不飽和カルボン酸と炭素数が2以上22以下のα,β-不飽和カルボン酸エステルとの三元共重合体が挙げられる。二元共重合体及び三元共重合体において、好ましいα-オレフィンはエチレン及びプロピレンであり、好ましいα,β-不飽和カルボン酸はアクリル酸及びメタクリル酸である。二元共重合体及び三元共重合体において、カルボキシル基の一部は金属イオンで中和されている。中和のための金属イオンとしては、ナトリウムイオン、カリウムイオン、リチウムイオン、亜鉛イオン、カルシウムイオン、マグネシウムイオン、アルミニウムイオン及びネオジムイオンが例示される。アイオノマー樹脂に代えて、又はアイオノマー樹脂と共に、他のポリマーが用いられてもよい。他のポリマーとしては、熱可塑性ポリウレタンエラストマー、熱可塑性スチレンエラストマー、熱可塑性ポリアミドエラストマー、熱可塑性ポリエステルエラストマー及び熱可塑性ポリオレフィンエラストマーが例示される。カバー6には、必要に応じ、二酸化チタンのような着色剤、硫酸バリウムのような充填剤、分散剤、酸化防止剤、紫外線吸収剤、光安定剤、蛍光剤、蛍光増白剤等が適量配合される。比重調整の目的で、カバー6にタングステン、モリブデン等の高比重金属の粉末が配合されてもよい。カバー6の厚みは0.3mm以上、特には0.5mm以上である。カバー6の厚みは2.5mm以下、特には2.2mm以下である。カバー6の比重は0.90以上、特には0.95以上である。カバー6の比重は1.10以下、特には1.05以下である。カバー6が2以上の層から構成されてもよい。図2は、図1のゴルフボール2が示された拡大正面図である。図2から明らかなように、多数のクレーター8がランダムに配置されている。これらクレーター8とランド10とにより、ゴルフボール2の表面に凹凸パターンが形成されている。この凹凸パターンの設計には、セルオートマトンCellularAutomaton法が用いられている。セルオートマトン法は、計算可能性理論、数学、理論生物学等の分野で、広く利用されている。セルオートマトン法のモデルは、多数のセルと単純な規則とからなる。このモデルにより、生命現象、結晶の成長、乱流等の自然現象が模擬されうる。このモデルでは、それぞれのセルが状態を有する。この状態は、ステージの進行に応じ、他の状態に変化しうる。あるセルのステージt1における状態は、ステージtにおける当該セルの状態及びこのセルの近傍にある複数のセルの状態によって決定される。決定は、ある規則に従ってなされる。全てのセルに、この規則が等しく適用される。この凹凸パターンの設計には、セルオートマトン法の反応拡散モデルが適している。このモデルは、獣、鳥、魚、昆虫等の体表面の模様の模擬に用いられている。このモデルでは、複数の状態が想定される。状態の数は、通常は2以上8以下である。それぞれのセルにおいて、初期の状態が決定される。ステージの進行により、規則に基づいて、状態が更新される。更新により、その状態が変化するセルが存在する。更新により、状態が変化しないセルも存在する。セルオートマトン法は、セルオートマトン法複雑系の自己組織化と超並列処理加藤泰義ら著、森北出版株式会社発行の第25-28頁に開示されている。本発明に係る設計方法の特徴は、セルの状態が、当該セルの近傍に位置する他のセルの影響下にて更新される点にある。この更新により、多数のクレーター8がランダムに配置された凹凸パターンが得られる。この特徴が維持される限り、如何なるモデルも用いられうる。本発明に係る設計方法は、効率の観点から、コンピュータとソフトウエアとが用いられて実施されることが好ましい。もちろん、手計算でも本発明は実施されうる。本発明の本質がコンピュータソフトウエアにあるわけではない。以下、セルオートマトン法の反応拡散モデルを用いた設計方法が詳説される。図3は、本発明の一実施形態に係る設計方法が示されたフローチャートである。図4は、図3の設計方法に用いられるメッシュ12が示された正面図である。このメッシュ12の形成には、球14が仮想されるSTEP1。この仮想球14の直径は、ゴルフボール2の直径と同一である。この仮想球14の表面が、多数の三角形に分割されるSTEP2。分割は、前進先端法advancingfrontmethodに基づいてなされている。前進先端法が、大学院情報理工学3計算力学伊藤耿一編、講談社発行の第195-197頁に開示されている。このメッシュ12において、三角形の数は176528個であり、頂点の数は88266個である。それぞれの頂点は、セル又はセルの中心と定義される。このメッシュ12では、セルの数は88266個である。他の手法によって仮想球14が分割されてもよい。この設計方法では、分化及び未分化の、2つの状態が想定される。それぞれのセルにおいて、いずれかの状態初期の状態が決定されるSTEP3。決定は、好ましくは、無作為になされる。無作為な決定には、乱数と剰余系とが用いられる。状態の数が2なので、基数が2である剰余系が用いられる。具体的には、コンピュータによって0以上1未満であり、下5桁の乱数が発生させられる。この乱数に100000が乗され、この積がさらに2で除される。この商の余りは、1又は0である。この余りに基づいて、当該セルの状態が決定される。例えば、余りが1の場合に分化が選定され、余りが0の場合に未分化が選定される。全てのセルについて、決定がなされる。決定の後のメッシュ12は、ステージ1にある。それぞれのセルにおいて、状態の変更の要否が判定されるSTEP4。判定は、規則に従ってなされる。図5は、この規則が説明されるためのグラフである。このグラフにおいて、縦軸は濃度であり、横軸は指数半径である。指数半径は、当該セルからの距離が基準値で除された値である。この基準値は、当該セルに最も近いセルと、当該セルとの距離である。濃度W1は正であり、濃度W2は負である。濃度W1の絶対値は、濃度W2の絶対値よりも大きい。指数半径R2は、指数半径R1よりも大きい。指数半径が0を超えてR1以下のエリアでは、濃度はW1である。指数半径がR1を超えてR2以下のエリアでは、濃度はW2である。図6は、図4のメッシュ12の一部が示された拡大図である。便宜上、図6では、メッシュ12が二次元で画かれている。図6の中心には、判定の対象であるセル16が示されている。図6にはさらに、第一円18及び第二円20が示されている。第一円18は、セル16を中心とする、指数半径がR1である円である。第二円20は、セル16を中心とする、指数半径がR2である円である。塗りつぶされた円で示されているのは、第一円18に含まれる、セル16以外のセルである。塗りつぶされた四角形で示されているのは、第二円20に含まれ、第一円18に含まれないセルである。塗りつぶされた三角形で示されているのは、第二円20に含まれないセルである。この設計方法では、第一円18に含まれこの第一円18の中心に位置しない、特定の状態にあるセルの数NR1がカウントされる。好ましい態様によれば、その状態が分化であるセルの数がカウントされ、合計NR1が算出される。この設計方法ではさらに、第二円20に含まれ第一円18に含まれない、特定の状態にあるセルの数NR1-R2がカウントされる。好ましい態様によれば、その状態が分化であるセルの数がカウントされ、合計NR1-R2が算出される。この数NR1及びNR1-R2が下記数式1に代入され、値Eが算出される。この値Eに基づいて、セル16の状態の変更の要否が決定される。E=W1NR1W2NR1-R21この判定に基づき、セル16の状態の更新がなされるSTEP5。更新において、セル16状態が変化することも、変化しないことも生じうる。好ましい態様によれば、値Eが正であるとき、セル16の状態が分化であればこの状態が維持され、セル16の状態が未分化であればこの状態が分化に変更される。この値Eがゼロであるとき、セル16の状態が維持される。この値Eが負であるとき、セル16の状態が分化であればこの状態が未分化に変更され、セル16の状態が未分化であればこの状態が維持される。全てのセルについての第一回目の更新が完了したメッシュ12は、ステージ2にある。以下に、判定及び更新の計算例が示される。条件第一濃度W1:1.00第二濃度W2:-0.60第一円18に含まれる、その状態が分化であるセルの数当該セル16を除く:8第二円20に含まれ第一円18に含まれない、その状態が分化であるセルの数:13計算例E=1.008-0.6013=0.2この場合、値Eが正なので、セル16の状態が分化であればこの状態が維持され、セル16の状態が未分化であればこの状態が分化に変更される。以下に、判定及び更新の他の計算例が示される。条件第一濃度W1:1.00第二濃度W2:-0.60第一円18に含まれる、その状態が分化であるセルの数当該セル16を除く:5第二円20に含まれ第一円18に含まれない、その状態が分化であるセルの数:9計算例E=1.005-0.609=-0.4この場合、値Eが負なので、セル16の状態が分化であればこの状態が未分化に変更され、セル16の状態が未分化であればこの状態が維持される。この判定と更新とが、繰り返される。図3のフローチャートでは、繰り返し数はNである。N回の繰り返しが完了した後のメッシュ12は、ステージN1にある。ステージの進行に伴い、更新によって状態が変化するセルの数が減少する。図7は、繰り返しによる状態の変化が示された写真である。なお、図7の写真は、図2に示されたゴルフボール2のものではない。図7では、便宜上、分化のセルが黒に着色され、未分化のセルが白に着色されている。各写真の詳細は、下記の通りである。ステージ1繰り返し数:0初期状態ステージ2繰り返し数:1ステージ3繰り返し数:2ステージ4繰り返し数:3ステージ5繰り返し数:4ステージ6繰り返し数:5ステージ31繰り返し数:30図7から明らかなように、繰り返し数が小さいステージでは、更新によるパーターンの変化が激しい。多数回の更新がなされることで、パターンが収束する。図7から明らかなように、繰り返し数は3以上が好ましく、5以上がより好ましい。繰り返し数が過大であると、コンピュータの負荷が大きい。この観点から、繰り返し数は30以下が好ましく、10以下がより好ましい。繰り返しの完了後、凹凸パターンが形成されるSTEP6。この凹凸パターンの形成では、それぞれのセルに、その状態に応じて、ランド10又はリセスが割り当てられる。具体的には、分化のセルにはリセスが割り当てられ、未分化のセルにはランド10が割り当てられる。分化のセルにランド10が割り当てられ、未分化のセルにリセスが割り当てられてもよい。以下、ランド10及びリセスの割り当て方法が説明される。図8は、更新が完了したのちのメッシュの一部が示された拡大図である。図9は、図8のIX-IX線に沿った断面図である。図8及び図9において、円で示されているのは分化のセルであり、四角形で示されているのは未分化のセルである。図9に示されるように、仮想球14の半径はRaである。図9には、第二の球22も示されている。この第二の球22は、仮想球14と同心である。第二の球22の半径Rbは、半径Raよりも小さい。分化のセル円で示されるは、仮想球14の半径方向に沿って、第二の球22の表面に移動させられる。移動距離は、Ra-Rbである。図9において三角形で示されているのが、移動後のセルである。半径Rbが半径Raよりも小さいので、この移動はリセスの割り当てに相当する。図9に示されるように、複数のリセスの連続より、クレーター8が形成される。クレーター8が、単一のリセスからなってもよい。未分化のセル四角形で示されるは、移動しない。移動のないことは、ランド10の割り当てに相当する。リセス及びランド10の割り当てにより、凹凸パターンが形成される。分化のセルと未分化のセルとが隣接する箇所では、未分化のセルから移動後の分化のセルに向かってスロープ24が形成される。スロープ24が円弧状であってもよい。1つのセルの移動距離が、他のセルの移動距離と異なってもよい。この場合、クレーター8は平坦でない底面を備えうる。図2に示されたゴルフボール2では、クレーター8がランダムに配置されている。この凹凸パターンは、方向性を有しない。このゴルフボール2は、空力的対称性に優れる。この設計方法では、状態変更の要否判定に影響を与える事項は、下記の2つである。I当該セル16の状態II当該セル16の近傍に位置する他のセルの状態この設計方法では、以下の仮定に基づき、判定がなされる。1分化のセルは、当該セル16の分化を促す活性化物質を生成する。2分化のセルは、当該セル16を分化状態から未分化の状態に戻す阻害物質を生成する。3未分化のセルは、当該セル16に影響を与えない。4活性化物質は大きな影響力をもち、阻害物質は小さな影響力を持つ。5活性化物質は遠方まで拡散せず、阻害物質は遠方まで拡散する。図5に示されたグラフでは、濃度W1の絶対値は、濃度W2の絶対値よりも大きい。これには、上記4の仮定が反映されている。図5に示されたグラフでは、指数半径R2は、指数半径R1よりも大きい。これには、上記5の仮定が反映されている。活性化物質の影響と、阻害物質の影響とは、互いに打ち消し合う。その状態が未分化であるセル16に働く活性化物質の影響が、このセル16に働く阻害物質の影響よりも大きい場合、このセル16が分化する。その状態が分化であるセル16に働く阻害物質の影響が、このセル16に働く活性化物質の影響よりも大きい場合、このセル16の状態が未分化に変化する。セルの数、第一濃度W1、第二濃度W2、指数半径R1及び指数半径R2は、パターンに影響を与える因子である。図10には、11個のパターンが示されている。便宜上、分化のセルが黒に着色され、未分化のセルが白に着色されている。それぞれのパターンの因子が、下記の表1に示されている。その幅が過大でないクレーター8が得られるとの観点から、セルの数は5000個以上が好ましく、9000個以上がより好ましく、20000個以上がさらに好ましく、40000個以上が特に好ましい。その幅が過小でないクレーター8が得られるとの観点から、セルの数は100000個以下が好ましく、90000個以下がより好ましい。幅が適正であるクレーター8は、乱流剥離を促進する。第一濃度W1は0.80以上が好ましく、0.95以上がより好ましい。第一濃度W1は1.20以下が好ましく、1.05以下がより好ましい。第二濃度W2は-0.70以上が好ましく、-0.65以上がより好ましい。第二濃度W2は-0.50以下が好ましく、-0.55以下がより好ましい。指数半径R1は2.20以上が好ましく、2.45以上がより好ましい。指数半径R1は2.60以下が好ましく、2.55以下がより好ましい。指数半径R2は3.0以上が好ましく、3.5以上がより好ましい。指数半径R2は5.0以下が好ましく、4.5以下がより好ましい。本発明では、全てのクレーター8の面積の合計の、仮想球14の表面積に対する比率は、占有率と称される。乱流剥離の促進の観点から、占有率は55以上が好ましく、65以上がより好ましく、82以上が特に好ましい。占有率は、95以下が好ましい。ゴルフボール2のホップが抑制されるとの観点から、クレーター8の深さRa-Rbは0.05mm以上が好ましく、0.08mm以上がより好ましく、0.10mm以上が特に好ましい。ゴルフボール2のドロップが抑制されるとの観点から、深さRa-Rbは0.60mm以下が好ましく、0.45mm以下がより好ましく、0.40mm以下が特に好ましい。本発明においてクレーターの容積とは、仮想球14とクレーター8の表面とに囲まれた部分の容積を意味する。全てのクレーター8の容積の和総容積は、ゴルフボール2のホップが抑制されるとの観点から400mm3以上が好ましく、450mm3以上がより好ましく、500mm3以上が特に好ましい。ゴルフボール2のドロップが抑制されるとの観点から、総容積は700mm3以下が好ましく、650mm3以下がより好ましく、600mm3以下が特に好ましい。好ましくは、このゴルフボール2の差dRの絶対値は、2.5mm以下である。この絶対値は、ゴルフボール2の空力的対称性と相関するパラメータである。この絶対値が小さいほど、PH回転時の弾道とPOP回転時の弾道との差が小さい。この観点から、絶対値は1.0mm以下がより好ましく、0.80mm以下が特に好ましい。以下、差dRに基づく評価方法が説明される。図11は、この評価方法が説明されるための模式図である。この評価方法では、第一回転軸Ax1が想定される。この第一回転軸Ax1は、ゴルフボール2の2つの極点Poを通過する。それぞれの極点Poは、ゴルフボール2の成形に用いられるモールドの最深部である。一方の極点Poは上型の最深部であり、他方の極点Poは下型の最深部である。ゴルフボール2は、第一回転軸Ax1を中心として回転する。この回転は、PH回転である。このゴルフボール2の仮想球14の表面に存在し、かつ第一回転軸Ax1と直交する大円GCが想定される。ゴルフボール2の回転のとき、この大円GCの周速が最も速い。さらに、ゴルフボール2の仮想球14の表面に存在し、第一回転軸Ax1と直交する2つの小円C1、C2が想定される。図12には、図11のゴルフボール2の一部の断面が模式的に示されている。図12の左右方向は、軸方向である。図12に示されるように、小円C1と大円GCとの中心角の絶対値は、30°である。図示されていないが、小円C2と大円GCとの中心角の絶対値も、30°である。これらの小円C1、C2により上記仮想球14が区画され、仮想球14の表面のうちこれら小円C1、C2に挟まれた領域が特定される。図12における点Pαは、ゴルフボール2の表面に位置し、かつ大円GCとの中心角がα°degreeである点である。点Fαは、点Pαから第一回転軸Ax1に下ろした垂線Peαの足である。矢印L1αで示されているのは、垂線Peαの長さである。換言すれば、長さL1αは、点Pαと第一回転軸Ax1との距離である。1つの断面において、21個の点Pαに関し、長さL1αが算出される。具体的には、-30°、-27°、-24°、-21°、-18°、-15°、-12°、-9°、-6°、-3°、0°、3°、6°、9°、12°、15°、18°、21°、24°、27°及び30°の角度αに関し、長さL1αが算出される。21個の長さL1αが合計され、総長さL2mmが得られる。総長さL2は、図12に示された断面における、表面の形状に依存するパラメータである。図13には、ゴルフボール2の一部の断面が示されている。図13において紙面垂直方向が、軸方向である。図13において符号βで示されているのは、ゴルフボール2の回転角度である。0°以上360°未満の範囲において、0.25°刻みに、回転角度βが設定される。それぞれの回転角度ごとに、総長さL2が算出される。この結果、回転方向に沿って1440の総長さL2が得られる。換言すれば、ゴルフボール2の1回転によって所定箇所に刻々と出現する表面の形状に依存するパラメータに関するデータ群が、算出される。このデータ群は、30240個の長さL1に基づいて算出されたものである。図2に示されたゴルフボール2のデータ群がプロットされたグラフが、図14に示されている。このグラフでは、横軸は回転角度βであり、縦軸は総長さL2である。このグラフから、総長さL2の最大値と最小値とが決定される。この最大値から最小値が減じられ、変動幅Rhが算出される。変動幅Rhは、PH回転における空力特性を表す数値である。さらに、第一回転軸Ax1と直交する第二回転軸Ax2が決定される。第二回転軸Ax2を中心としたゴルフボール2の回転は、POP回転である。PH回転と同様、POP回転についても、大円GCと2つの小円C1、C2が想定される。小円C1と大円GCとの中心角の絶対値は、30°である。小円C2と大円GCとの中心角の絶対値も、30°である。仮想球14の表面のうちこれら小円に挟まれた領域において、1440の総長さL2が算出される。換言すれば、ゴルフボール2の1回転によって所定箇所に刻々と出現する表面の形状に依存するパラメータに関するデータ群が、算出される。図2に示されたゴルフボール2のデータ群がプロットされたグラフが、図15に示されている。このグラフでは、横軸は回転角度βであり、縦軸は総長さL2である。このグラフから、総長さL2の最大値と最小値とが決定される。この最大値から最小値が減じられ、変動幅Roが算出される。変動幅Roは、POP回転における空力特性を表す数値である。変動幅Rhから変動幅Roが減じられ、差dRが算出される。差dRは、ゴルフボール2の空力的対称性を表す指標である。本発明者が得た知見によれば、差dRが小さいゴルフボール2は、空力的対称性に優れる。その理由は、PH回転時の表面形状とPOP回転時の表面形状との類似性が高いためと推測される。第一回転軸Ax1と直交する直線は、無数に存在する。従って、大円GCも無数に存在する。クレーター8に含まれる部分が最長である大円GCが選択され、変動幅Ro及び差dRが算出される。これに代えて、無作為に抽出された20個の大円GCに基づき、20個の変動幅Roが算出されてもよい。この場合、20個のデータの中の最大値に基づいて、差dRが算出される。変動幅Rhが小さいほど、PH回転時に大きな飛距離が得られうる。この観点から、変動幅Rhは3.0mm以下が好ましく、2.8mm以下がより好ましい。変動幅Roが小さいほど、POP回転時に大きな飛距離が得られうる。この観点から、変動幅Roは3.0mm以下が好ましく、2.7mm以下がより好ましい。以下、実施例によって本発明の効果が明らかにされるが、この実施例の記載に基づいて本発明が限定的に解釈されるべきではない。[実施例1]図3に示された方法により、図2に示された凹凸パターンを設計した。パターンに影響を与える因子の詳細は、以下の通りである。セルの数:47415個第一濃度W1:1.00第二濃度W2:-0.60指数半径R1:2.5指数半径R2:4.0クレーターの深さRa-Rb:0.213mm[実施例2から4]因子の値を下記の表2に示される通りとした他は実施例1と同様にして、実施例2から4の凹凸パターンを設計した。実施例2の凹凸パターンが、図16に示されている。実施例3の凹凸パターンが、図17に示されている。実施例4の凹凸パターンが、図18に示されている。[比較例1]図19に示されたディンプルパターンを設計した。図19には、1つのユニットに関し、ディンプルの種類が符号で示されている。このユニットは、球面が10分割されることで得られる。このユニットのパターンが、球面全体に展開される。このディンプルパターンは、直径が4.00mmであるディンプルAと、直径が3.70mmであるディンプルBと、直径が3.40mmであるディンプルCと、直径が3.20mmであるディンプルDとを備えている。それぞれのディンプルの断面形状は、円弧である。ディンプルの詳細は、下記の通りである。種類数直径深さ容積A1204.000.1841.737B1523.700.1841.414C603.400.1841.137D603.200.1840.977[評価]前述の手法により、各パターンの差dRを算出した。結果の詳細が、下記の表3に示されている。表3に示されるように、各実施例の差dRの絶対値は、比較例1のそれよりも小さい。この評価結果から、本発明の優位性は明らかである。"