I have a contour plot that repeats itself. All I am only interested in
the contour that is centred around (1,0) to show up on the plot.
nothing else. Here is the code I used to plot the contour.
r,z,gro,kro = var('r,z,gro,kro')
gro = 2.405
kro = 3
def f(r,z): term1 = (1+(bessel_J(0, gro)/bessel_J(2,gro)))*(r/
1)*(bessel_J(1, gro*r)/bessel_J(1,gro))*cos(kro*z); term2 = -
(bessel_J(0,gro)/bessel_J(2, gro))*(r**2/1**2); return term1 + term2
cp = contour_plot(f , (r, 0, 2.6), (z, -1.1, 1.1), axes_labels=['$r
$','$z$'], contours = [0.0, 0.2, 0.4, 0.6, 0.8], labels = True,
label_inline= True, label_colors = 'black', label_inline_spacing =
0.01, fill=False, aspect_ratio=1, legend_label='Polodial')
show(cp)
--
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
