ここはRやWolframによる解法のコード置き場。
亀レスとの照合用。