Tes

1 minute read

test for jupyter

def trapezoidal(f, a, b, n):
    h = float(b - a) / n
    s = 0.0
    s += f(a)/2.0
    for i in range(1, n):
        s += f(a + i*h)
    s += f(b)/2.0
    return s * h
print( trapezoidal(lambda x:x**2, 5, 10, 100))
291.66875000000005
x = (lambda x:x**2)
# a = [.2,.3,.4,.5,.6,.7,.8,.9]
# b = [.2,.35,.34,.45,.46,.57,.68,.79]
# print( trapezoidal(lambda x:x**2, a, b, len(a)))
# h = (list(set(a)-set(b)))
import numpy as np
# a = 0
# b = np.pi
# n = 11
# h = (b - a) / (n - 1)
# x = np.linspace(a, b, n)
# f = np.sin(x)
# I_trap = (h/2)*(f[0] + 2 * sum(f[1:n-1]) + f[n-1])
# err_trap = 2 - I_trap
# print(I_trap)
# print(err_trap)
# e = 2.718281828459045
# X = 0.1
# (e**(X*1j)).imag   # sin(X)
# 0.09983341664682815
# (e**(X*1j)).real   # cos(X)
# 0.9950041652780258