what is my problem in numerov method code?

this is my hole code for vibrating string

clc
  clear
  disp 'numerov method for finding zero of function'
  k=1;
  dk=0.5;
  tol=1e-08;
  phim=0;
  phiz=0.01;
  phipold=funcnuv(k,phiz,phim);

        for i=1:90
            phip=funcnuv(k,phiz,phim);
            phim=phiz;
            phiz=phip;
            k=k+dk;
                  if(phipold*phip<0)
                    k=k-dk;
                    dk=dk/2;
                  end
                  if (abs(dk)<tol)
                      break
                  end
        end

    disp 'phip=' ; disp(phip)
    disp 'k='; disp(k)
    disp 'i='; disp(i)

and this is function that i have used

function phip=funcnuv(k,phiz,phim)
h=0.01;
const=(k*h)*(k*h)/12;
phip=(2*(1-5*const)*phiz-(1+const)*phim)/(1+const);

[Pawel Ladosz]

The code runs just fine, so not sure what you mean by problem? You not getting right answer?

Answer 1

[Zakir Khan]

Detail is given here: https://www.iue.tuwien.ac.at/phd/brech/ch_2_2.htm, I implemented the code. https://www.fiverr.com/users/tala1971

[Zakir Khan]

wonderful

Answer 2

[Zakir Khan]

Here it is given at the link: https://fr.mathworks.com/matlabcentral/fileexchange/23384-gaas-gaalas-quantum-wells which I tested it my self. https://www.upwork.com/freelancers/~015508901db48c1b4c