fun=@(x) 1+cos(x).*cosh(x);
xx=(0:9)*pi+pi/2;
z=zeros(size(xx));
for i=1:length(xx)
z(i)=fzero(fun,xx(i));
end
因為cos和cosh都是偶函數(shù),所以它的根正負(fù)對稱的
這個函數(shù)是個超越函數(shù),所以解析解不能用初等函數(shù)表示
我們可以重數(shù)值解法
因為cos(x)的根是pi/2+n*pi,可以預(yù)見函數(shù)的根在 pi/2+n*pi附近
以pi/2+n*pi為初始值,尋找附近的根就可以了
得到的z就是10個最靠近0點的正數(shù)根近似值,它們的相反數(shù)也是根
注意第一個根的誤差最小,后面的根的誤差會越來越大我輸入之后為什么沒有結(jié)果顯示呢？加一句disp(z)