File:Stability region for Euler method.svg
Jump to navigation
Jump to search
Size of this PNG preview of this SVG file: 720 × 540 pixels. Other resolutions: 320 × 240 pixels | 640 × 480 pixels | 1,024 × 768 pixels | 1,280 × 960 pixels | 2,560 × 1,920 pixels.
Original file (SVG file, nominally 720 × 540 pixels, file size: 23 KB)
Summary
| Description |
English: Stability region for Euler method |
| Date | |
| Source | Own work |
| Author | Helmut Podhaisky |
| Other versions | Euler method -- Trapzoidal method -- Adams Bashforth |
| SVG development |  This plot was created with Matplotlib. |
| Source code | Python codefrom pylab import linspace,pi,exp,real,imag,axis
import matplotlib.path as mpath
import matplotlib.patches as mpatches
import matplotlib.pyplot as plt
w=exp(1j*linspace(0,2*pi,200))
Path=mpath.Path
fig = plt.figure()
for name,z in zip(['Euler','trapezoidal', 'AB2'],\
[ w-1, [5j,-5+5j,-5-5j,-5j,5j],2*(w**2-w)/(3*w-1)]):
fig.clf()
ax=fig.add_subplot(111)
verts=map(lambda z:(real(z),imag(z)), z)
codes=[Path.MOVETO]+[Path.LINETO]*(len(verts)-2)+[Path.CLOSEPOLY]
path=mpath.Path(verts,codes)
patch=mpatches.PathPatch(path,facecolor=[1, 0.5, 0.8],edgecolor='black',alpha=1)
ax.add_patch(patch)
ax.plot([-3,2],[0,0],'--',color='black')
ax.plot([0,0],[-2,2],'--',color='black')
ax.set_xlim(-3,2)
ax.set_ylim(-2,2)
fig.savefig("Stability region for "+name+" method.svg")
|
Licensing
|
I, the copyright holder of this work, release this work into the public domain. This applies worldwide. In some countries this may not be legally possible; if so: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law. |
File history
Click on a date/time to view the file as it appeared at that time.
| Date/Time | Thumbnail | Dimensions | User | Comment | |
|---|---|---|---|---|---|
| current | 11:46, 26 August 2010 | | 720 × 540 (23 KB) | wikimediacommons>Helmut.podhaisky | no |
File usage
There are no pages that use this file.
