function z=f(x,y)
global k l abj
z=p_jacobi01(x,k)*p_jacobi01(y,l);
