Someone asked about the reason for looking for the roots. In case the
original question didn't get answered, or the problem has been changed,
the problem arises in getting the eigenfunctions (natural mode shapes)
for a cantilevered, built in vibrating uniform rod (or beam). After the
first two the rest come very close to (2*n+1)pi/2 because cosh(x) gets
so big that only when cos(x) is virtually zero does the solution appear.
Mike.
=========
On 30/03/2014 21:55, Paul CARRICO wrote:
Thanks ... indeed much more stable ..
Paul
-----Message d'origine-----
De : users [mailto:[email protected]] De la part de Rafael Guerra
Envoyé : dimanche 30 mars 2014 23:34
À : 'International users mailing list for Scilab.'
Objet : [Scilab-users] RE(2): finding roots
Hi Paul,
The new code here below does not show any such problem for all N solutions I
tried.
Note that your problem: cos(x) .* cosh(x) + 1 = 0; is equivalent to: cos(x) +
sech(x) = 0.
The latter form seems to be numerically more stable.
Regards,
Rafael
.............
_______________________________________________
users mailing list
[email protected]
http://lists.scilab.org/mailman/listinfo/users
_______________________________________________
users mailing list
[email protected]
http://lists.scilab.org/mailman/listinfo/users
