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);
1 个评论
Pawel Ladosz
2016-9-29
The code runs just fine, so not sure what you mean by problem? You not getting right answer?
回答（2 个）
Zakir Khan
2021-8-10
0 个投票
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
2021-8-10
0 个投票
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
类别
在 帮助中心 和 File Exchange 中查找有关 Vibration Analysis 的更多信息
2016-9-28
2021-8-10
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
