@morotesovtannu function ret = plane_jacob(theta, mu, lambda)
sin_theta = sin(theta);
cos_theta = cos(theta);
if abs(cos_theta) > 1.0e-10
ret = [-sin_theta/cos_theta^2 -mu*cos_theta*sin_theta/cos_theta^3;...
-sin_theta/cos_theta^2 -lambda*cos_theta*sin_theta/cos_theta^3];
else
ret = [-1/sin_theta^2 -mu*sin_theta/sin_theta^3;...
-1/sin_theta^2 -lambda*sin_theta/sin_theta^3];
end