Dmitry wrote: > > Fellows, > if somebody uses float point, could you check that everything works Ok. Hi Dimtry, This might be easier to debug. Obviously it only works for 0 <= x <= pi/2
double cosine(double x)
{
double y;
y = 1 + x*(-0.4999999963 + x*(0.416666418 + x*(-0.0013888397 +
x*(0.0000247609 +
x * -0.0000002605))));
return y;
}
But why can't I find sqrt when I
#include <math.h>
double sine(double x)
{
return sqrt(1 - cosine(x)*cosine(x));
}
Garst
