问题描述
我需要使用区间算术来模拟几个模型,我发现最可行的软件包是:mpmath。但是我在绘制图形时遇到问题。在将其应用于模型之前,我进行了初始测试。有人可以帮助我吗?
这是我的代码:
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无法正常工作。
解决方法
这相当简单:
import numpy as np
import matplotlib.pyplot as plt
x = np.arange(0,5,0.01)
y = np.sin(x)
plt.plot(x,y)
plt.show()