In [1]:
import  matplotlib.pyplot as plt
from numpy import *
In [7]:
ax1 = plt.subplot(5,2,1)
ax2 = plt.subplot(5,2,2)
ax3 = plt.subplot(5,2,10)

t = linspace(0,2*pi,300)
for c in linspace(0,1,10):
    ax1.plot(c*cos(t),c*sin(t))
    ax2.plot(t,c*t**2)
    ax3.plot(t,c*t**2)
plt.show()
In [8]:
x = 1.e300
x
Out[8]:
1e+300
In [9]:
x*x
Out[9]:
inf
In [11]:
x = linspace(-2e-15,2e-15,2000)
y = ((1+x)-1)/x
plt.plot(x,y,'.')
plt.show()
In [13]:
x = linspace(-2e-15,2e-15,2000)
y = log(1+x)/x
ya = 1 -x/2 + x**2/3 -x*3/4
plt.plot(x,y,'b.')
plt.plot(x,ya,'r')
plt.show()
In [17]:
# larger window
x = linspace(-2e-1,2e-1,2000)
y = log(1+x)/x
ya = 1 -x/2 + x**2/3 -x**3/4
plt.plot(x,y,'b.')
plt.plot(x,ya,'y',lw=3)
plt.show()
In [20]:
# larger window
%matplotlib notebook

x = linspace(-2e-1,2,2000)
y = log(1+x)/x
ya = 1 -x/2 + x**2/3 -x**3/4
yabetter =  1 -x/2 + x**2/3 -x**3/4 + x**4/5 - x**5/6 + x**6/7
plt.plot(x,y,'b.')
plt.plot(x,ya,'y',lw=3)
plt.plot(x,yabetter,'c',lw=3)
plt.show()