function [wy]=wartwl(a,x,y,z)
wart=zeros(3,3);
ww=zeros(1,3);
wart(1,1)=(a(1,1)+a(1,4)*y+a(1,5)*z+a(1,7)*y*z)-(3*x^2+y^2+z^2);
wart(1,2)=(a(1,2)+a(1,4)*x+a(1,6)*z+a(1,7)*x*z)-2*x*y;
wart(1,3)=(a(1,3)+a(1,5)*x+a(1,6)*y+a(1,7)*y*x)-2*x*z;
wart(2,1)=(a(2,1)+a(2,4)*y+a(2,5)*z+a(2,7)*y*z)-2*x*y;
wart(2,2)=(a(2,2)+a(2,4)*x+a(2,6)*z+a(2,7)*x*z)-(3*y^2+x^2+z^2);
wart(2,3)=(a(2,3)+a(2,5)*x+a(2,6)*y+a(2,7)*y*x)-2*z*y;
wart(3,1)=(a(3,1)+a(3,4)*y+a(3,5)*z+a(3,7)*y*z)-2*x*y;
wart(3,2)=(a(3,2)+a(3,4)*x+a(3,6)*z+a(3,7)*x*z)-2*y*z;
wart(3,3)=(a(3,3)+a(3,5)*x+a(3,6)*y+a(3,7)*y*x)-(3*z^2+y^2+z^2);
wy=eig(wart)';
%disp('------------------------------')
%disp(eig(wart)')
%disp('------------------------------')
