"複数の状態が想定されるステップ、球面上に多数のセルが想定されるステップ、それぞれのセルに、いずれかの状態が付与されるステップ、このセルの状態及びこのセルの近傍に位置する複数のセルの状態に基づいて、当該セルの属性として、インサイド、アウトサイド及び境界のいずれかが付与されるステップ、多数のセルの属性に基づいてクレーターが想定されるステップ並びに上記クレーターの面積が拡大するように、それぞれのセルの属性及びこのセルの近傍に位置する複数のセルの属性に基づいて当該セルの属性が更新されるステップを含む、ゴルフボール表面の凹凸パターン設計方法であって、上記セルの状態の変更の要否の判定及びセルの状態の更新が、セルオートマトン法の反応拡散モデルによってなされる設計方法。上記クレーターが想定されるステップが、その属性がインサイド又は境界であるセルにリセスが割り当てられるステップ及びその属性がアウトサイドであるセルにランドが割り当てられるステップを含んでおり、複数のリセスの集合によってクレーターが想定される請求項1に記載の設計方法。上記クレーターが想定されるステップが、その属性がインサイドであるセルにリセスが割り当てられるステップ及びその属性が境界又はアウトサイドであるセルにランドが割り当てられるステップを含んでおり、複数のリセスの集合によってクレーターが想定される請求項1に記載の設計方法。上記クレーターの面積が算出されて属性の更新の要否が判定されるステップをさらに含む請求項1から3のいずれかに記載の設計方法。上記属性の更新が複数回繰り返される請求項1から4のいずれかに記載の設計方法。上記セルの属性の更新のステップにおいて、境界のセルに隣接するアウトサイドのセルの属性が境界に変更される請求項1から5のいずれかに記載の設計方法。上記セルへの状態の付与が、それぞれのセルの初期の状態が決定されるステップ、このセルの状態の変更の要否が、このセルの近傍に位置する複数のセルの状態に基づいて判定されるステップ及び上記判定に基づき、このセルの状態の更新がなされるステップによってなされる請求項1から6のいずれかに記載の設計方法。上記初期の状態が無作為に決定される請求項7に記載の設計方法。上記セルの状態の変更の要否の判定及びセルの状態の更新が3回以上繰り返される請求項7又は8に記載の設計方法。上記セルの数が5000個以上300000個以下である請求項7から9のいずれかに記載の設計方法。下記数式1で算出される値Eに基づいて、状態の変更の要否が判定される請求項10から10のいずれかに記載の設計方法。E=W1NR1W2NR1-R21この数式1において、W1は第一濃度であり、NR1は第一円に含まれこの第一円の中心に位置しない特定の状態にあるセルの数であり、W2は第二濃度であり、NR1-R2は第二円に含まれ第一円に含まれない特定の状態にあるセルの数である。第一濃度は正であり、第二濃度は負である。第一円は、当該セルを中心とする、指数半径がR1である円である。第二円は、当該セルを中心とする、指数半径がR2である円である。半径R2は、半径R1よりも大きい。それぞれのセルの初期の状態が分化又は未分化であり、下記数式1で算出される値Eが正であるとき、セルの状態が分化であればこの状態が維持され、セルの状態が未分化であればこの状態が分化に変更され、この値Eがゼロであるとき、セルの状態が維持され、この値Eが負であるとき、セルの状態が分化であればこの状態が未分化に変更され、セルの状態が未分化であればこの状態が維持される請求項11に記載の設計方法。E=W1NR1W2NR1-R21この数式1において、W1は第一濃度であり、NR1は第一円に含まれこの第一円の中心に位置しない分化のセルの合計数であり、W2は第二濃度であり、NR1-R2は第二円に含まれ第一円に含まれない分化のセルの合計数である。第一濃度は正であり、第二濃度は負である。第一円は、当該セルを中心とする、指数半径がR1である円である。第二円は、当該セルを中心とする、指数半径がR2である円である。半径R2は、半径R1よりも大きい。上記第一濃度W1が0.80以上1.20以下である請求項11又は12に記載の設計方法。上記第二濃度W2が-0.70以上-0.50以下である請求項11から13のいずれかに記載の設計方法。上記指数半径R1が2.20以上5.0以下である請求項11から14のいずれかに記載の設計方法。上記指数半径R2が3.0以上10.0以下である請求項11から15のいずれかに記載の設計方法。その表面に多数のクレーターを備えており、これらクレーターのパターンが請求項1から16のいずれかに記載の方法で設計されたゴルフボール。下記ステップ1から17によって得られる差dRの絶対値が2.5mm以下である請求項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以下である請求項18に記載のゴルフボール。上記変動幅Rhが3.0mm以下であり、上記変動幅Roが3.0mm以下である請求項18又は19に記載のゴルフボール。上記クレーターの総面積の、ゴルフボールの仮想球の表面積に対する比率が65以上である請求項17から20のいずれかに記載のゴルフボール。特開昭61-284264号公報に開示されたゴルフボールでは、ディンプルパターンに起因する不都合が、容積の相違によって解消されている。ディンプルパターンに起因する不都合が、ディンプルパターン自体の工夫で解消されているわけではない。このゴルフボールでは、ディンプルパターンが本来備えるポテンシャルが犠牲にされる。このゴルフボールの飛距離は、十分ではない。赤道の近傍の特異性、並びにこれに起因する不十分な対称性に関し、原因究明の努力がなされている。しかし、未だ原因は明らかでなく、改善のための普遍的な理論も確立されていない。本発明の目的は、占有率が高く、飛行性能及び空力的対称性に優れたゴルフボールの提供にある。本発明に係る凹凸パターン設計方法は、複数の状態が想定されるステップ、球面上に多数のセルが想定されるステップ、それぞれのセルに、いずれかの状態が付与されるステップ、このセルの状態及びこのセルの近傍に位置する複数のセルの状態に基づいて、当該セルの属性として、インサイド、アウトサイド及び境界のいずれかが付与されるステップ、多数のセルの属性に基づいてクレーターが想定されるステップ並びに上記クレーターの面積が拡大するように、それぞれのセルの属性及びこのセルの近傍に位置する複数のセルの属性に基づいて当該セルの属性が更新されるステップを含む。好ましくは、クレーターが想定されるステップは、その属性がインサイド又は境界であるセルにリセスが割り当てられるステップ及びその属性がアウトサイドであるセルにランドが割り当てられるステップを含む。複数のリセスの集合により、クレーターが想定される。クレーターが想定されるステップが、その属性がインサイドであるセルにリセスが割り当てられるステップ及びその属性が境界又はアウトサイドであるセルにランドが割り当てられるステップを含んでもよい。複数のリセスの集合により、クレーターが想定される。好ましくは、この設計方法は、クレーターの面積が算出されて属性の更新の要否が判定されるステップをさらに含む。好ましくは、この設計方法において、属性の更新が複数回繰り返される。好ましくは、セルの属性の更新のステップにおいて、境界のセルに隣接するアウトサイドのセルの属性が境界に変更される。好ましくは、セルへの状態の付与は、それぞれのセルの初期の状態が決定されるステップ、このセルの状態の変更の要否が、このセルの近傍に位置する複数のセルの状態に基づいて判定されるステップ及び上記判定に基づき、このセルの状態の更新がなされるステップによってなされる。好ましくは、初期の状態は、無作為に決定される。好ましくは、セルの状態の変更の要否の判定及びセルの状態の更新が3回以上繰り返された後に、それぞれのセルに、属性が付与される。好ましくは、セルの状態の変更の要否の判定及びセルの状態の更新は、セルオートマトン法によってなされる。好ましくは、この判定及び更新は、セルオートマトン法の反応拡散モデルによってなされる。好ましくは、セルの数は5000個以上300000個以下である。好ましくは、状態の変更の要否は、下記数式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以上5.0以下である。好ましくは、指数半径R2は、3.0以上10.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以下である。好ましくは、変動幅Rh及び変動幅Roのそれぞれは、3.0mm以下である。好ましくは、クレーターの総面積の、ゴルフボールの仮想球の表面積に対する比率は、65以上である。本発明に係るゴルフボールの表面には、多数のクレーターが形成される。これらクレーターは、乱流剥離を促進する。これらクレーターは、ゴルフボールの飛行性能に寄与する。これらクレーターはランダムに配置されるので、このパターンは方向性を有さない。このゴルフボールは、空力的対称性に優れる。クレーターの占有率は、大きい。以下、適宜図面が参照されつつ、好ましい実施形態に基づいて本発明が詳細に説明される。図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の中心には、判定の対象であるセル16aが示されている。図6にはさらに、第一円18及び第二円20が示されている。第一円18は、セル16aを中心とする、指数半径がR1である円である。第二円20は、セル16aを中心とする、指数半径がR2である円である。塗りつぶされた円で示されているのは、第一円18に含まれる、セル16a以外のセルである。塗りつぶされた四角形で示されているのは、第二円20に含まれ、第一円18に含まれないセルである。塗りつぶされた三角形で示されているのは、第二円20に含まれないセルである。この設計方法では、第一円18に含まれこの第一円18の中心に位置しない、特定の状態にあるセルの数NR1がカウントされる。好ましい態様によれば、その状態が分化であるセルの数がカウントされ、合計NR1が算出される。この設計方法ではさらに、第二円20に含まれ第一円18に含まれない、特定の状態にあるセルの数NR1-R2がカウントされる。好ましい態様によれば、その状態が分化であるセルの数がカウントされ、合計NR1-R2が算出される。この数NR1及びNR1-R2が下記数式1に代入され、値Eが算出される。この値Eに基づいて、セル16aの状態の変更の要否が決定される。E=W1NR1W2NR1-R21この判定に基づき、セル16aの状態の更新がなされるSTEP5。更新において、セル16a状態が変化することも、変化しないことも生じうる。好ましい態様によれば、値Eが正であるとき、セル16aの状態が分化であればこの状態が維持され、セル16aの状態が未分化であればこの状態が分化に変更される。この値Eがゼロであるとき、セル16aの状態が維持される。この値Eが負であるとき、セル16aの状態が分化であればこの状態が未分化に変更され、セル16aの状態が未分化であればこの状態が維持される。全てのセルについての第一回目の更新が完了したメッシュ12は、ステージ2にある。以下に、判定及び更新の計算例が示される。条件第一濃度W1:1.00第二濃度W2:-0.60第一円18に含まれる、その状態が分化であるセルの数当該セル16aを除く:8第二円20に含まれ第一円18に含まれない、その状態が分化であるセルの数:13計算例E=1.008-0.6013=0.2この場合、値Eが正なので、セル16aの状態が分化であればこの状態が維持され、セル16aの状態が未分化であればこの状態が分化に変更される。以下に、判定及び更新の他の計算例が示される。条件第一濃度W1:1.00第二濃度W2:-0.60第一円18に含まれる、その状態が分化であるセルの数当該セル16aを除く:5第二円20に含まれ第一円18に含まれない、その状態が分化であるセルの数:9計算例E=1.005-0.609=-0.4この場合、値Eが負なので、セル16aの状態が分化であればこの状態が未分化に変更され、セル16aの状態が未分化であればこの状態が維持される。この判定と更新とが、繰り返される。図3のフローチャートでは、繰り返し数はMである。M回の繰り返しが完了した後のメッシュ12は、ステージM1にある。ステージの進行に伴い、更新によって状態が変化するセルの数が減少する。図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以下がより好ましい。判定と更新とがM回繰り返されることにより、それぞれのセルの状態が確定する。この確定は、セルへの状態の付与である。図8は、状態の付与が完了した後のメッシュの一部が示された拡大図である。図8において、円で示されているのは分化のセルであり、四角形で示されているのは未分化のセルである。この状態に基づき、セルにiflagが付与される。まず暫定的に、全てのセルにiflagとして0が付与される。次に、状態が分化であるセルに関し、iflagが変更される。図8において符号16bで示されたセルは、6個のセル16c-16hと隣接している。本発明では、一方のセルを頂点とする三角形の他の頂点に他方のセルが存在するとき、一方のセルが他方のセルと隣接すると称される。これらのセル16c-16hの状態は、分化である。隣接する全てのセル16c-16hの状態が分化であるとき、当該セル16bのiflagが0から1に変更される。図8において符号16nで示されたセルは、6個のセル16h-16mと隣接している。セル16h、16i、16l、16mの状態は、分化である。セル16j、16kの状態は、未分化である。状態が未分化である1又は2以上のセルと隣接するとき、当該セル16nのiflagが0から2に変更される。状態が分化である全てのセルに関し、そのiflagが変更される。状態が未分化であるセルのiflagは、変更されない。このiflagに基づき、全てのセルに属性が付与されるSTEP6。属性の付与は、下記のルールに基づいてなされる。iflag:0属性:アウトサイドiflag:1属性:インサイドiflag:2属性:境界属性の付与が完了したメッシュ12は、第一フェーズにある。属性が境界である複数のセルを結ぶことにより、輪郭21第一輪郭が得られる。図8において、第一輪郭21が太線で示されている。それぞれのセルに、その属性に応じて、ランド10又はリセスが割り当てられるSTEP7。具体的には、属性がアウトサイドであるセルにランド10が割り当てられ、属性がインサイドであるセルにリセスが割り当てられ、属性が境界であるセルにリセスが割り当てられる。図9は、図8のIX-IX線に沿った断面図である。図9に示されるように、仮想球14の半径はRaである。図9には、第二の球22も示されている。この第二の球22は、仮想球14と同心である。第二の球22の半径Rbは、半径Raよりも小さい。属性がインサイドであるセル16g、16hは、仮想球14の半径方向に沿って、第二の球22の表面に移動させられる。属性が境界であるセル16m、16oも、仮想球14の半径方向に沿って、第二の球22の表面に移動させられる。セル16g、16h、16m、16oの移動距離は、Ra-Rbである。図9において三角形で示されているのが、移動後のセルである。半径Rbが半径Raよりも小さいので、この移動はリセスの割り当てに相当する。図9に示されるように、複数のリセスの連続より、クレーター8が形成される。クレーター8が、単一のリセスからなってもよい。未分化のセル四角形で示されるは、移動しない。移動のないことは、ランド10の割り当てに相当する。リセス及びランド10の割り当てにより、凹凸パターンが形成される。ランドのセルと境界のセルとが隣接する箇所では、ランドのセルから移動後の境界のセルに向かってスロープ24が形成される。スロープ24が円弧状であってもよい。1つのセルの移動距離が、他のセルの移動距離と異なってもよい。この場合、クレーター8は平坦でない底面を備えうる。この実施形態では、クレーター8は、スロープ24及び底面26からなる。図8に示された第一輪郭21は、底面26の外縁である。図9から明らかなように、クレーター8の外縁は、第一輪郭21よりも大きい。クレーター8の外縁は、第一輪郭21に依存して決まる。クレーター8が、他の断面形状を有してもよい。このクレーター8が形成されたのちの凹凸パターンが、図10に示されている。このパターンは、多数のクレーター8とランド10とからなる。このパターンの占有率が算出されるSTEP8。この算出では、クレーター8の輪郭で囲まれた面積が算出される。全てのクレーター8の面積が合計される。この合計の、仮想球14の表面積に対する比率が、占有率である。図4に示された多数の三角形が用いられ、占有率が近似的に算出されてもよい。近似的な算出では、クレーター8に含まる三角形の合計面積が、全ての三角形の合計面積で除される。得られた占有率に基づき、判定がなされるSTEP9。このステップでは、占有率が所定値以上であるか否かが判定される。図3に示された実施形態では、占有率Yが65以上であるか否かが判定される。占有率Yが65未満である場合、属性の更新がなされるSTEP10。以下、この更新の方法が詳説される。図11は、属性の付与が完了した後のメッシュの一部が示された拡大図である。符号16nで示されたセルは、第一輪郭21の上に存在する。このセル16nには、6個のセル16hから16mが隣接している。セル16hのiflagは1であり、その属性はインサイドである。属性がインサイドであるセルでは、iflagの変更はなされない。セル16i、16l、16mのiflagは2であり、その属性は境界である。その属性が境界であり、かつ属性が境界である他のセルと隣接するセルでは、iflagの変更はなされない。セル16j、16kのiflagは0であり、その属性はアウトサイドである。その属性がアウトサイドであり、かつ属性が境界である他のセルと隣接するセルでは、iflagが0から3に変更される。第一輪郭21の上に存在する全てのセルに関し、このセルと隣接するセルのiflagが決定される。このiflagに基づき、属性の更新STEP10がなされる。属性の更新は、下記のルールに基づいてなされる。iflag:0属性:アウトサイドiflag:1-2属性:インサイドiflag:3属性:境界属性の更新が1回なされた後のメッシュ12は、第二フェーズにある。属性が境界である複数のセルを結ぶことにより、輪郭28第二輪郭が得られる。それぞれのセルに、その属性に応じて、ランド10又はリセスが割り当てられるSTEP7。具体的には、属性がアウトサイドであるセルにランド10が割り当てられ、属性がインサイドであるセルにリセスが割り当てられ、属性が境界であるセルにリセスが割り当てられる。図12は、図11のXII-XII線に沿った断面図である。図12では、リセスの割り当てにより、クレーター8が形成されている。図9及び12の対比から明らかなように、属性の更新によりクレーター8が拡張される。換言すれば、属性の更新により占有率が大きくなる。このクレーター8が形成されたのちの凹凸パターンが、図13に示されている。図10及び13の対比から明らかなように、図13のパターンの占有率は、図10のそれよりも大きい。このパターンの占有率が算出されるSTEP8。得られた占有率に基づき、判定がなされるSTEP9。このステップでは、占有率が所定値以上であるか否かが判定される。図3に示された実施形態では、占有率Yが65以上であるか否かが判定される。以下同様に、占有率Yが65以上となるまで、属性の更新STEP10、リセス又はランド10の割り当てSTEP7、占有率の算出STEP8及び判定STEP9が繰り返される。第N回目の属性更新に先立ち、その属性がアウトサイドであり、かつ属性が境界である他のセルと隣接するセルにおいて、iflagが0からN2に変更される。第N回目の属性更新は、下記のルールに基づいてなされる。iflag:0属性:アウトサイドiflag:1からN1まで属性:インサイドiflag:N2属性:境界属性の更新がN回なされた後のメッシュ12は、第N1フェーズにある。図2に示されたゴルフボール2のパターンは、2回の属性の更新がなされて得られたものである。換言すれば、このパターンのメッシュ12は、第三フェーズにある。図2、10及び13の対比から明らかなように、このゴルフボール2の占有率は大きい。このゴルフボール2の占有率は、79である。大きな占有率を備えたゴルフボール2は、飛行性能に優れる。飛行性能の観点から、占有率は65以上が好ましく、75以上がより好ましく、80以上が特に好ましい。図2に示されたゴルフボール2では、クレーター8がランダムに配置されている。この凹凸パターンは、方向性を有しない。このゴルフボール2は、空力的対称性に優れる。図14は、3回の属性の更新がなされて得られたパターンが示された正面図である。換言すれば、このパターンのメッシュ12は、第四フェーズにある。このパターンでは、隣接するクレーター8同士が合体している。このパターンの占有率は極めて高い。このパターンでは、ランド10の面積が極めて小さい。ランド10の面積が小さなゴルフボール2がゴルフクラブで打撃されると、ランド10が摩滅しやすい。ランド10の耐久性の観点から、過大な占有率は好ましくない。ランド10の耐久性の観点から、占有率は95以下が好ましく、90以下がより好ましく、87以下が特に好ましい。隣接するクレーター8同士の合体を避ける目的で、他のクレーター8との間隔が小さいクレーター8が拡張されず、残余のクレーター8のみが拡張されてもよい。この設計方法では、その状態が未分化であるセルに近接しかつ自らの状態が分化であるセルに、属性として境界が付与される。その状態が分化であるセルに近接しかつ自らの状態が未分化であるセルに、属性として境界が付与されてもよい。この設計方法では、状態変更の要否判定に影響を与える事項は、下記の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は、パターンに影響を与える因子である。その幅が過大でないクレーター8が得られるとの観点から、セルの数は5000個以上が好ましく、9000個以上がより好ましく、20000個以上がさらに好ましく、40000個以上が特に好ましい。その幅が過小でないクレーター8が得られるとの観点から、セルの数は300000個以下が好ましい。幅が適正であるクレーター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は5.0以下が好ましく、4.6以下がより好ましい。指数半径R2は3.0以上が好ましく、3.5以上がより好ましい。指数半径R2は10.0以下が好ましく、8.0以下がより好ましい。ゴルフボール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に基づく評価方法が説明される。図15は、この評価方法が説明されるための模式図である。この評価方法では、第一回転軸Ax1が想定される。この第一回転軸Ax1は、ゴルフボール2の2つの極点Poを通過する。それぞれの極点Poは、ゴルフボール2の成形に用いられるモールドの最深部である。一方の極点Poは上型の最深部であり、他方の極点Poは下型の最深部である。ゴルフボール2は、第一回転軸Ax1を中心として回転する。この回転は、PH回転である。このゴルフボール2の仮想球14の表面に存在し、かつ第一回転軸Ax1と直交する大円GCが想定される。ゴルフボール2の回転のとき、この大円GCの周速が最も速い。さらに、ゴルフボール2の仮想球14の表面に存在し、第一回転軸Ax1と直交する2つの小円C1、C2が想定される。図16には、図15のゴルフボール2の一部の断面が模式的に示されている。図16の左右方向は、軸方向である。図16に示されるように、小円C1と大円GCとの中心角の絶対値は、30°である。図示されていないが、小円C2と大円GCとの中心角の絶対値も、30°である。これらの小円C1、C2により上記仮想球14が区画され、仮想球14の表面のうちこれら小円C1、C2に挟まれた領域が特定される。図16における点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は、図16に示された断面における、表面の形状に依存するパラメータである。図17には、ゴルフボール2の一部の断面が示されている。図17において紙面垂直方向が、軸方向である。図17において符号βで示されているのは、ゴルフボール2の回転角度である。0°以上360°未満の範囲において、0.25°刻みに、回転角度βが設定される。それぞれの回転角度ごとに、総長さL2が算出される。この結果、回転方向に沿って1440の総長さL2が得られる。換言すれば、ゴルフボール2の1回転によって所定箇所に刻々と出現する表面の形状に依存するパラメータに関するデータ群が、算出される。このデータ群は、30240個の長さL1に基づいて算出されたものである。図2に示されたゴルフボール2のデータ群がプロットされたグラフが、図18に示されている。このグラフでは、横軸は回転角度βであり、縦軸は総長さ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のデータ群がプロットされたグラフが、図19に示されている。このグラフでは、横軸は回転角度βであり、縦軸は総長さ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.8mm以下がより好ましい。PH回転時及びPOP回転時のいずれにおいても大きな飛距離が得られうるとの観点から、変動幅Rh及び変動幅Roの両方が3.0mm以下であることが好ましく、2.8mm以下であることがより好ましい。この設計方法では、セルへの状態の付与がセルオートマトン法によってなされている。他の方法によってセルへの状態の付与がなされてもよい。以下、実施例によって本発明の効果が明らかにされるが、この実施例の記載に基づいて本発明が限定的に解釈されるべきではない。[実施例1]図3に示された方法により、図2に示された凹凸パターンを設計した。パターンに影響を与える因子の詳細は、以下の通りである。セルの数:88266個第一濃度W1:1.00第二濃度W2:-0.60指数半径R1:4.6指数半径R2:8.0セルの属性の更新:2回第三フェーズクレーターの深さRa-Rb:0.1352mm[実施例2-3及び比較例1]属性の更新の回数とクレーター深さとを下記の表1に示される通りとした他は実施例1と同様にして、実施例2-3及び比較例1の凹凸パターンを設計した。[比較例2]図20に示されたディンプルパターンを設計した。図20には、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[実施例2-3及び比較例1]属性の更新の回数とクレーター深さとを下記の表1に示される通りとした他は実施例1と同様にして、実施例2-3及び比較例1の凹凸パターンを設計した。[実施例4]図3に示された方法により、図22に示された凹凸パターンを設計した。パターンに影響を与える因子の詳細は、以下の通りである。セルの数:157045個第一濃度W1:1.00第二濃度W2:-0.60指数半径R1:4.6指数半径R2:8.0セルの属性の更新:1回第二フェーズクレーターの深さRa-Rb:0.1814mm[実施例5-6及び比較例3]属性の更新の回数とクレーター深さとを下記の表2に示される通りとした他は実施例4と同様にして、実施例5-6及び比較例3の凹凸パターンを設計した。[実施例7]図3に示された方法により、図26に示された凹凸パターンを設計した。パターンに影響を与える因子の詳細は、以下の通りである。セルの数:279329個第一濃度W1:1.00第二濃度W2:-0.60指数半径R1:4.6指数半径R2:8.0セルの属性の更新:1回第二フェーズクレーターの深さRa-Rb:0.1814mm[実施例8-9及び比較例4]属性の更新の回数とクレーター深さとを下記の表3に示される通りとした他は実施例7と同様にして、実施例8-9及び比較例4の凹凸パターンを設計した。[評価]前述の手法により、各パターンの差dRを算出した。結果の詳細が、下記の表1から3に示されている。表1から3に示されるように、各実施例のRh及びRoは小さい。しかも、各実施例の差dRの絶対値は、比較例2のそれよりも小さい。この評価結果から、本発明の優位性は明らかである。"