阪本ひろむ氏から（その５６）についての追加コメントが届いたので紹介する．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

　すでに，これまでの方法でうまく描画できない理由として浮動小数点のアンダーフローの話をしたが，

　　ｒ^(n/2)＝ｃｏｓ（ｎ／２・ｔ）

の場合，他に事情がある．

　プロットする区間はｔが［０，２π］である．プロット関数はこの区間をｔ0，ｔ1，ｔ2，・・・と適当に分割し，その後，点ｒ（ｔ0），ｒ（ｔ1），ｒ（ｔ2），・・・を結合して曲線を描く．

　ところが，この場合，ｒ＝０となるｔの付近でのｒの変化は非常に大きく

　　ｄｒ／ｄｔ＝±∞

である．よって，区間をどのように細分しても，これらのｔjに対して，ｒ（ｔj）＝０となる点がないのである．

　ｔj-1までｒは実数値，ｔjで虚数となる場合，描画はｔj-1のところで途切れてしまう．ｔj-1までｒは虚数値，ｔjで実数となる場合，描画はｔjから始まる．よって，ｒ＝０の付近での描画はできない．

　そこで，

　　ｒ＝ｃｏｓ（ｎ／２・ｔ）^2/n　　　ｃｏｓ（ｎ／２・ｔ）≧０のとき

　　ｒ＝０　　　ｃｏｓ（ｎ／２・ｔ）＜０のとき

なる条件を加えて，ｒ＝０までの描画をせざるを得ないのである．

　ところで，今回のケースでは，ある区間でｒ＝０という条件を設定できた．そうでない場合，たとえば，

　　ｒ＝ｃｏｓ（ｎ／２・ｔ）^2/n　　　ｃｏｓ（ｎ／２・ｔ）≧０のとき

　　ｒ＝（−ｃｏｓ（ｎ／２・ｔ））^2/n　　　ｃｏｓ（ｎ／２・ｔ）＜０のとき

とした場合，この手はもはや通用しない（ｎ／２が偶数の場合，この式が正解）．やはりｒ＝０付近での描画はうまくできない．このような場合は

［１］ｒ≧０となる区間をあらかじめ設定しておく．

［２］それ以外のところで，ｒ＝０となるように関数を定義．

［３］これらのプロット図を重ね合わせるという方法をとることになる．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝