def __call__(self, x):
from scitools.std import linspace, sum, zeros,iseq
f, a, n = self.f, self.a, self.n
h = (x-a)/float(n)
I = 0.5*f(a)
ar = zeros(len(linspace(1,n-1, n-1)))
for i in iseq(1, n-1):
ar[i-1] = f(a + i*h)
I = sum(ar)
I += 0.5*f(x)
I *= h
return I
def __call__(self, x):
from scitools.std import linspace, sum, zeros, iseq
f, a, n = self.f, self.a, self.n
h = (x - a) / float(n)
I = 0.5 * f(a)
ar = zeros(len(linspace(1, n - 1, n - 1)))
for i in iseq(1, n - 1):
ar[i - 1] = f(a + i * h)
I = sum(ar)
I += 0.5 * f(x)
I *= h
return I
def __call__(self, x):
from scitools.std import iseq
f, a, n = self.f, self.a, self.n
h = (x-a)/float(n)
I = 0.5*f(a)
l = []
for i in iseq(1, n-1):
I = f(a + i*h)
l.append(I)
I = sum(l)
I += 0.5*f(x)
I *= h
return I
def __call__(self, x):
from scitools.std import iseq
f, a, n = self.f, self.a, self.n
h = (x - a) / float(n)
I = 0.5 * f(a)
l = []
for i in iseq(1, n - 1):
I = f(a + i * h)
l.append(I)
I = sum(l)
I += 0.5 * f(x)
I *= h
return I
