Beispiel #1
0
 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
Beispiel #2
0
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