When I evaluate the code shown below, it gives an error:
NDSolve::ndsz: At x == -0.00001, step size is effectively zero; singularity or stiff system suspected.
W := (2*Pi*3*10^14)
CC := 3*10^8
vt := ((2*T/me)^(1/2))
me := 9.11*10^-31
T := 1.6*10^-15
Wc := 1.6*10^-19/me
L := 0.00001
wp := W*Sqrt[(1 - x/L)]
Om := Wc
e0 := 8.85*10^-12
LnLumbda := 10
k := W/CC (1 - (1 - x/L)/(1 + (Wc)/W))^(1/2)
R := ((1 - x/L)*W^2*(1.6*10^-19)^2*Pi*Sqrt[me]*LnLumbda)/( T^(3/2)*e0)
q := (W + 2 I R - Om)/vt
F := q/k
sol =
NDSolve[
{(y''[x] +
(W^2/ (CC^2) + ((W (wp)^2)/((CC^2) vt k ))
((E^-F^2 (-π Erfi[F] + Log[-(1/F)] + Log[F]))/Sqrt[π])) y[x]) == 0,
y'[-L] == 0, y[-L] == 10},
y, {x, -L, L}]
I will appreciate any help.