您所在位置：主题：利用fzero循环求解y＝tan（y）
bjb05
[专家分：0] 发布于 2007-11-25 12:48:00
y=tan(y)在整个数轴上有无穷解，我只要求[0;60*pi]区间上的
下面是我的程序：
>>a=0.01;
>>for i=0:59;
>>X=fzero(@(y) y-tan(y),a);
>>a=X+pi;
>>t(i+1)=X;
>>end
不知道什么原因，只能求出3个解，到7.多时马上回到1.5几，下面是我的结果
-0.000000006571222 4.493409457909064 7.725251836937707
1.570796326794896 4.712388980384688 7.725251836937707
1.570796326794896 4.712388980384688…………
1.5几是pi/2,为啥就算到他了呢？
而且第二个解和后面的4.7几不一样，为什么会这样，tan（4.712388980384688）＝5.1019e+014也不等于4.712388980384688啊
请指教
下面是我的程序：
>>a=0.01;
>>for i=0:59;
>>X=fzero(@(y) y-tan(y),a);
>>a=X+pi;
>>t(i+1)=X;
>>end
不知道什么原因，只能求出3个解，到7.多时马上回到1.5几，下面是我的结果
-0.000000006571222 4.493409457909064 7.725251836937707
1.570796326794896 4.712388980384688 7.725251836937707
1.570796326794896 4.712388980384688…………
1.5几是pi/2,为啥就算到他了呢？
而且第二个解和后面的4.7几不一样，为什么会这样，tan（4.712388980384688）＝5.1019e+014也不等于4.712388980384688啊
请指教