这是我为使用牛顿法制作分形而编写的一个小脚本。importnumpyasnpimportmatplotlib.pyplotaspltf=np.poly1d([1,0,0,-1])#x^3-1fp=np.polyder(f)defnewton(i,guess):ifabs(f(guess))>.00001:returnnewton(i+1,guess-f(guess)/fp(guess))else:returnipic=[]foryinnp.linspace(-10,10,1000):pic.append([newton(0,x+y*1j)forxinnp.linspace(-10,10