The R Project SVN R

Rev

Rev 81560 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 81560 Rev 86647
Line 144... Line 144...
144 
#endif
 144 
#endif
145 
 
 145 
 
146 
#define n_max (100)
 146 
#define n_max (100)
147 
 
 147 
 
148 
// Compute  d_n(x) = (d/dx)^n cot(x)  ; cot(x) := cos(x) / sin(x)
 148 
// Compute  d_n(x) = (d/dx)^n cot(x)  ; cot(x) := cos(x) / sin(x)
149 
double d_n_cot(double x, int n)
 149 
static double d_n_cot(double x, int n)
150 
{
 150 
{
151 
    if (n == 0)
 151 
    if (n == 0)
152 
	return cos(x)/sin(x);
 152 
	return cos(x)/sin(x);
153 
    else if (n == 1) // -1/sin^2
 153 
    else if (n == 1) // -1/sin^2
154 
	return -1/R_pow_di(sin(x), 2);
 154 
	return -1/R_pow_di(sin(x), 2);