This version works for me:
r,z = var('r,z')
gro = 1.0
kro = 3.0
def f(r,z):
term1 = (1+(bessel_J(0, gro)/bessel_J(2,gro)))*(r/
kro)*(bessel_J(1, r)/bessel_J(1,gro))*cos(z)
term2 = -(bessel_J(0,r)/bessel_J(2, gro))*(r**2/kro**2)
return term1 + term2
p = contour_plot(f , (r, 0, 2.5), (z, 0, 4), 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(p)
Andrzej Chrzeszczyk
--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
