NumerovLin9[fcn_,x0_,xe_,ystart0_,ystart1_,n_]:=
Module[{a,b,c,h,x,m,y,f,f0,f1},
a={{0,0,0,0,0,0,0,0},
   {0,0,0,0,0,0,0,0},
   {-0.598998471886051805067313262522717,0.505248471886051805067313262522717,0,0,0,0,0,0},
   {0.611940860843265466836612918938395,-0.457598599630250945661536894626551,0.0019077387869854788249239756881560630,0,0,0,0,0},
   {5.427666026511682927591036891473420,0,-0.14260519141008493021164346639423363,-5.393130752150043699996669546434617,0,0,0,0},
   {-5.610135701858784405155111668019878,0,0.5506567074629963937612723050126895,5.636262548192080006678383019228344,-0.000842150804071192265661109162914162,0,0,0},
   {-0.136636677728244444569934203059859,0.561663202085575874213271190674143,0,-0.5341258126703869729688709410700263,-0.0010317928621100818145281444837034207,-0.00098002993594548597104901317166502217,0,0},
   {3.3958135995527209430054039199949006,0,0,0,0.77107979659239720331001070341560209,0.4029386965596218656312998325642190,-3.5698320927047400119467144559747217,0}};
b={0.0052416891284815813117699910152740341,0.078091150623339464541181279378704271,0.33333215391419725785050862450243255,0.33333215391419725785050862450243255,
   0.12238058164565142856713074479294128,0.12238058164565142856713074479294128,0,0.0052416891284815813117699910152740341};
c=-{1,0,1/4,-1/4,0.68401132004066650563615866841367152,-0.6840113200406665056361586684136715,1/3,-1};
h=SetPrecision[(xe-x0)/n,33];
s=Length[c];
x=x0+Range[0,n]*h;
m=Length[ystart0];
y=Array[Array[0&,m]&,n+1];
y[[1]]=ystart0;
y[[2]]=ystart1;
f=Array[Array[0&,m]&,s];
f[[2]]=fcn[x[[1]],y[[1]]];
Do[f0=f[[2]];
   f1=fcn[x[[k]],y[[k]]];
   f=Array[Array[0&,m]&,s];
   f[[1]]=f0;f[[2]]=f1;
   Do[f[[o]]=fcn[x[[k]]+c[[o]]*h,(1+c[[o]])*y[[k]]-c[[o]]*y[[k-1]]+h^2*a[[o]].f],{o,3,s}];
   y[[k+1]]=SetPrecision[2*y[[k]]-y[[k-1]]+h^2*b.f,33],{k,2,n}];
Return[{x,y}]];