- 315 名前：デフォルトの名無しさん [2009/02/16(月) 19:15:48 ]
- 丸投げっぽくてアレだけどだれか教えてくれない？
余弦をマクローリン展開してるつもりなんだけどどうしても
マイナス方向に誤差がでる。原因は式？
それとも浮動小数点の誤差？
template<class type> type Cos(type n)
{
type r=0,pow=1.0;
unsigned long factorial=1,i=0;
n*=n;
while(i<18)
{
type x,y;
x=pow/factorial;
pow*=n;
factorial*=++i;
factorial*=++i;
y=pow/factorial;
pow*=n;
factorial*=++i;
factorial*=++i;
r+=x-y;
}
return r;
}
