a1 =        1.38;

a2 =      0.3108;

a3 =      0.6671;

a4 =      0.3192;

a5 =      0.9714;

a6 =      0.4783;

a7 =      -9.493;

a8 =      0.1478;

c1 =       1.094;

c2 =    0.001242;

c3 =      0.7892;

c4 =       1.287;

c5 =      0.2025;

c6 =      0.6889;

c7 =     0.03396;

c8 =      0.8076;


f2 = chebfun(@(x) (a1*exp(-((x)/c1)^2) + a2*abs(x)*exp(-((x)/c2)^2) + a3*x^2*exp(-((x)/c3)^2) + a4*abs(x)^3*exp(-((x)/c4)^2) + a5*x^4*exp(-((x)/c5)^2) + a6*abs(x)^5*exp(-((x)/c6)^2) + a7*x^6*exp(-((x)/c7)^2) + a8*abs(x)^7*exp(-((x)/c8)^2)), [-50,50]);

L = chebop(@(x,u) diff(u,2)/2-x*diff(u,1)/2-f2*u,[-50,50]);

L.bc = 'neumann';

[V,D] = eigs(L);

diag(D)

Warning: Function failed to evaluate on array inputs.
Vectorizing the function may speed up its evaluation
and avoid the need to loop over array elements.
Use 'vectorize' flag in the CHEBFUN constructor call
to avoid this warning message.
> In chebfun>vectorCheck (line 913)
  In chebfun>parseInputs/parseOp (line 764)
  In chebfun>parseInputs (line 755)
  In chebfun (line 203)
  In main (line 216)
Warning: Function not resolved using 65537 pts. Have you tried
'splitting on'?
> In chebfun/constructor>constructorNoSplit (line 120)
  In chebfun/constructor (line 63)
  In chebfun (line 219)
  In main (line 216)

ans =

   -2.8005
   -2.3387
   -1.8797
   -1.5001
   -1.0001
   -0.8883