algebrite/sources/sin.coffee

Computer Algebra System in Coffeescript

# Sine function of numerical and symbolic arguments



Eval_sin = ->
	#console.log "sin ---- "
	push(cadr(p1))
	Eval()
	sine()
	#console.log "sin end ---- "

sine = ->
	#console.log "sine ---- "
	save()
	p1 = pop()
	if (car(p1) == symbol(ADD))
		sine_of_angle_sum()
	else
		sine_of_angle()
	restore()
	#console.log "sine end ---- "

# Use angle sum formula for special angles.

#define A p3
#define B p4

sine_of_angle_sum = ->
	#console.log "sin of angle sum ---- "
	p2 = cdr(p1)
	while (iscons(p2))
		p4 = car(p2); # p4 is B
		if (isnpi(p4)) # p4 is B
			push(p1)
			push(p4); # p4 is B
			subtract()
			p3 = pop(); # p3 is A
			push(p3); # p3 is A
			sine()
			push(p4); # p4 is B
			cosine()
			multiply()
			push(p3); # p3 is A
			cosine()
			push(p4); # p4 is B
			sine()
			multiply()
			add()
			#console.log "sin of angle sum end ---- "
			return
		p2 = cdr(p2)
	sine_of_angle()
	#console.log "sin of angle sum end ---- "

sine_of_angle = ->

	if (car(p1) == symbol(ARCSIN))
		push(cadr(p1))
		return

	if (isdouble(p1))
		d = Math.sin(p1.d)
		if (Math.abs(d) < 1e-10)
			d = 0.0
		push_double(d)
		return

	# sine function is antisymmetric, sin(-x) = -sin(x)

	if (isnegative(p1))
		push(p1)
		negate()
		sine()
		negate()
		return

	# sin(arctan(x)) = x / sqrt(1 + x^2)

	# see p. 173 of the CRC Handbook of Mathematical Sciences

	if (car(p1) == symbol(ARCTAN))
		push(cadr(p1))
		push_integer(1)
		push(cadr(p1))
		push_integer(2)
		power()
		add()
		push_rational(-1, 2)
		power()
		multiply()
		return

	# multiply by 180/pi

	push(p1)
	push_integer(180)
	multiply()
	push_symbol(PI)
	divide()

	n = pop_integer()

	if (n < 0 || n == 0x80000000)
		push(symbol(SIN))
		push(p1)
		list(2)
		return

	switch (n % 360)
		when 0, 180
			push_integer(0)
		when 30, 150
			push_rational(1, 2)
		when 210, 330
			push_rational(-1, 2)
		when 45, 135
			push_rational(1, 2)
			push_integer(2)
			push_rational(1, 2)
			power()
			multiply()
		when 225, 315
			push_rational(-1, 2)
			push_integer(2)
			push_rational(1, 2)
			power()
			multiply()
		when 60, 120
			push_rational(1, 2)
			push_integer(3)
			push_rational(1, 2)
			power()
			multiply()
		when 240, 300
			push_rational(-1, 2)
			push_integer(3)
			push_rational(1, 2)
			power()
			multiply()
		when 90
			push_integer(1)
		when 270
			push_integer(-1)
		else
			push(symbol(SIN))
			push(p1)
			list(2)