def integrand(x): # . . . eller , om du er fancy : integrand = lambda x : f ( x ) *np . cos (n*x ) return f(x) * np.cos(n * x) integral = s.simpson(-np.pi, np.pi, 1000, integrand) if n == 0: return integral / ( 2.0 * np.pi ) else: return integral / np.pi
def sinprojection(n, f): def integrand(x): return f(x) * np.sin(n * x) return s.simpson(-np.pi, np.pi, 1000, integrand) / np.pi