%% Plot Jacobi Zeta Function
% Plot real and imaginary values of the Jacobi zeta function using fcontour. Set u on the x-axis and m on the y-axis by using the symbolic function f with the variable order (u,m). Fill plot contours by setting Fill to on.

syms f(u,m)
f(u,m) = jacobiZeta(u,m);

subplot(2,2,1)
fcontour(real(f),'Fill','on')
title('Real Values of Jacobi Zeta')
xlabel('u')
ylabel('m')

subplot(2,2,2)
fcontour(imag(f),'Fill','on')
title('Imaginary Values of Jacobi Zeta')
xlabel('u')
ylabel('m')

% Copyright 2012 The MathWorks, Inc.