matlab求频率方程的超越方程y1=1+cos(x).*cosh(x);y2=2.4*x.*(sin(x).*cosh(x)-cos(x).*sinh(x));x的取值范围是（0,30）.求两条曲线的交点坐标.因为是多解,希望能给出具体的程序算法.不要用作图法,因为作图法

matlab求频率方程的超越方程y1=1+cos(x).*cosh(x);y2=2.4*x.*(sin(x).*cosh(x)-cos(x).*sinh(x));x的取值范围是（0,30）.求两条曲线的交点坐标.因为是多解,希望能给出具体的程序算法.不要用作图法,因为作图法
matlab求频率方程的超越方程
y1=1+cos(x).*cosh(x);
y2=2.4*x.*(sin(x).*cosh(x)-cos(x).*sinh(x));
x的取值范围是（0,30）.
求两条曲线的交点坐标.因为是多解,希望能给出具体的程序算法.不要用作图法,因为作图法精度太低了,希望能得到小数点后至少六位的数值,双精度的就更好了.小弟摸索了一段时间,始终得不到数值.真心希望论坛里的老师们帮忙解除疑惑,

matlab求频率方程的超越方程y1=1+cos(x).*cosh(x);y2=2.4*x.*(sin(x).*cosh(x)-cos(x).*sinh(x));x的取值范围是（0,30）.求两条曲线的交点坐标.因为是多解,希望能给出具体的程序算法.不要用作图法,因为作图法
function f=fx(x)
f(1)=1+cos(x(1))*cosh(x(1))-x(2);
f(2)=2.4*x(1)*(sin(x(1))*cosh(x(1))-cos(x(1))*sinh(x(1)))-x(2);
>> [x,f]=fsolve(@fx,[1,30])
x =
1.0328 1.8109
f =
1.0e-010 *
-0.0955 0.9024
>> [x,f]=fsolve(@fx,[4,30])
x =
3.9738 -16.9100
f =
1.0e-013 *
0.0355 0.1776
>>[x,f]=fsolve(@fx,[9,30])
x =
7.0972 415.8972
f =
1.0e-011 *
0 -0.3013