function F = acrit(a,x,D2);
 
A  = -D2+diag(x.^4+i*a*x);
E = sort(eig(A));
 
F = abs(E(1)-E(2));
 
 

