Я написал код, чтобы использовать Римана для вычисления площади под определенной функцией.Это работает хорошо, но когда я устанавливаю границы на десятичное число, я получаю ошибку.Я думаю, что я должен сказать функции, которую я хочу, чтобы иметь возможность входить в поплавки, но я не понимаю, как.Я думаю, что это легко исправить, надеюсь, кто-то может мне помочь:
import numpy as np
import math
def func(x):
return x ** 2
def func1(x):
return x ** (x +0.5)
def func2(x):
return np.sin(x)
def func3(x):
return np.sin(x ** 2)
def riemann(f,a,b,N):
dx = (b-a)/N
x_midpoint = np.linspace(dx/2,b - dx/2,N)
riemann_sum = np.sum(f(x_midpoint) * dx)
return riemann_sum
res = riemann(func,0,1,10000)
res1 = riemann(func1,0,1,10000)
res2 = riemann(func2,0,1,10000)
res3 = riemann(func3,0,1,10000)
print("Midpoint Riemann Sum:",res)
print("Midpoint Riemann Sum:",res1)
print("Midpoint Riemann Sum:",res2)
print("Midpoint Riemann Sum:",res3)