🚀 go-pugleaf

I seem to be having trouble with the macsyma fourier package.

When I try to obtain sine coefficients for with
		F: A*X
	and
		F: A*X + B
I get the same result using FOURSIN(F,X,L);

		       -2 A L     N
		B   := ------ (-1)
		 N	N %PI

However when I perform the integrals myself with
	(2/L) * DEFINT(F*SIN(N*%PI*X/L),X,0,L);
I get something different for F: A*X + B, namely

			-2                 N
		B   := ----- [(A L + B)(-1)  - B]
		 N     N %PI		

The problem is easily found when one tries FOURSIN(A,X,L) and gets

		B  := 0
		 N
which is clearly incorrect.
Any suggestions to correct this problem would be appreciated.

/* Jeffery G. Smith, BS-RHIT (AKA Doc. Insomnia, WMHD-FM)       *
 *    The Ohio State University, Graduate Physics Program       *
 *        3193 Smith Lab, Columbus, OH 43210  (614) 292-5321    *
 *    smithj@ohstpy.mps.ohio-state.edu                          */