発言者: tDB
発言日: 2005 11/15 18:22
発言元: 222.158.93.151
> ありがとうございます。定義域が間違っていました。計算しなおすと
> r = ±sqrt( (4*sin(T)-cos(T)) / (sin(T)**3-cos(T)**3) )
> の定義域は、π/4 < T ≦ π+α、ただし α=arctan(1/4) と思うのですが合ってるでしょうか?
私もそのように思います。
が,
> T →π/4 のときに r →±∞ になるので、区間の左端を $pi/4+0.001 として
> 区間の右端を求めて次のように描くとエラーになりました。よろしくお願いします。
近似計算の面倒なところです。
> \def\At{X-X**3/3+X**5/5}
> \funcval\At{0.25}\aval
arctan(1/4) の計算としては,誤差が大きいでしょうか。
結果は
0.244987
一方,perl の逆三角関数 atan2 を用いた結果は
0.244979
もっともこれでもエラーが出ます。
丸め誤差恐るべし。
この値に対して,\Fx の根号の中を計算すると0ではなく
-0.000002
となってしまいます。
そこで
\calcval{atan2(1,4)}\aval
ではなく
\calcval{atan2(1,4)-0.000001}\aval
としてやることで回避します。
%----------------------------------------------
\documentclass{jarticle}
\usepackage{emathPp}
\begin{document}
\begin{zahyou}[ul=5mm](-8,8)(-8,8)
\calcval{atan2(1,4)-0.000001}\aval
\def\Ft{sqrt((4*sin(T)-cos(T))/((sin(T)-cos(T))*(1+sin(T)*cos(T))))}
\RGurafu(.001)\Ft{$pi/4+0.001}{$pi+\aval}
\RGurafu(.001){-\Ft}{$pi/4+0.001}{$pi+\aval}
\Hasen{(\xmin,\xmin)(\xmax,\xmax)}
\end{zahyou}%
\end{document}
▼関連発言
│
└◆3892:曲線のグラフについて [z] 11/09 16:49
├◆3894:--- [---] 11/09 17:47
├◆3895:Re:曲線のグラフについて [tDB] 11/09 20:27
│└◆3898:Re[2]:曲線のグラフについて [z] 11/10 01:15
├◆3951:Re:曲線のグラフについて [z] 11/15 01:21
│└◆3952:Re[2]:曲線のグラフについて [飯島 徹] 11/15 02:47
│ └◆3953:Re[3]:曲線のグラフについて [z] 11/15 03:03
└◆3954:Re:曲線のグラフについて [z] 11/15 07:11
└◆3956:Re[2]:曲線のグラフについて [tDB] 11/15 13:54
└◆3957:Re[3]:曲線のグラフについて [z] 11/15 16:34
└◆3958:Re[4]:曲線のグラフについて [tDB] 11/15 18:22
├◆3960:Re[5]:曲線のグラフについて [z] 11/15 22:00
└◆3961:emathPp.sty v 0.40 [tDB] 11/15 22:08
└◆3964:Re:emathPp.sty v 0.40 [z] 11/15 23:39<-last