// xsin.hpc -- Calculate the sine -*- C -*-

declare xsin(x) : s(0), n(0) {
	for (n = 0; n < 20; ++n)
		s += (-1)**n * x**(2*n + 1) / (2*n + 1)!;
	s;		// Return value
}