t1=@(x,y) ((x^2-2.10.*y+x.*sqrt(x^2-2.*2.10.*y))./2)^0.25;
t2=@(x,y) ((x^2-2.10.*y-x.*sqrt(x^2-2.*2.10.*y))./2)^0.25;
r=@(x,y) sqrt(x^2-2.*2.10.*y);
q1=@(x,y) -(x-r(x,y)).*(cos(t1(x,y).*y)+(exp(t1(x,y).*y)+exp(-t1(x,y).*y))./2)./(4.*r(x,y))+(x+r(x,y)).*(cos(t2(x,y).*y)+(exp(t2(x,y).*y)+exp(-t2(x,y).*y))./2)./(4.*r(x,y))+2.10.*y.*(-cos(t1(x,y).*y)+(exp(t1(x,y).*y)+exp(-t1(x,y).*y))./2)./(4.*r(x,y).*(t1(x,y))^2)-2.10.*y.*(-cos(t2(x,y).*y)+(exp(t2(x,y).*y)+exp(-t2(x,y).*y))./2)./(4.*r(x,y).*(t2(x,y))^2);
q2=@(x,y) -(x-r(x,y)).*(sin(t1(x,y).*y)+(exp(t1(x,y).*y)-exp(-t1(x,y).*y))./2)./(4.*r(x,y).*t1(x,y))+(x+r(x,y)).*(sin(t2(x,y).*y)+(exp(t2(x,y).*y)-exp(-t2(x,y).*y))./2)./(4.*r(x,y).*t2(x,y))+2.10.*y.*(-sin(t1(x,y).*y)+(exp(t1(x,y).*y)-exp(-t1(x,y).*y))./2)./(4.*r(x,y).*(t1(x,y))^3)-2.10.*y.*(-sin(t2(x,y).*y)&