我需要用区间算法模拟几个模型,我找到的最可行的包是:mpmath。但是我在绘制图形时遇到问题。在将其应用于模型之前,我进行了初步测试。有谁能够帮我?
- 另一个问题是我总是需要一个 for 来创建我的区间变量,这大大增加了计算成本。会有另一种选择吗?
这是我的代码:
import mpmath as mp
import math as mt
from mpmath import *
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
iv.dps = 10
iv.pretty = True
X = np.linspace(-np.pi, np.pi, 10, endpoint=True)
a=iv.mpf(X[1])
b=[]
for k in range(len(X)):
b = np.append(b,iv.mpf(X[k]) )
C=[]
for k in range(len(X)):
C = np.append(C, iv.sin(b[k]))
print(C)
- 我需要绘制罪恶,而 mp.plot 不起作用。