Question

请教如何用matlab解隐性方程：y=((9.8*T^2)/2*pi)*tanh((2*pi*d)/y)

Answer 1

1.直接用solve函数

T=input('');

d=input('');

f=subs('((9.8*T^2)/2*pi)*tanh((2*pi*d)/y)-y',{'T','d'},{T,d});

a=solve(f,'y')

》T=2,d=3

结果

a=-32.323847772664794935753715077212.

2.Newton迭代法

设定了迭代的最大次数100,可以自己改动的

T=input(''); d=input('');

syms xk x a i

fx=input('Please input a f(x):');

gx=x-fx/diff(fx)

x0=input('Please input a initial number x0=');

p =input('Please input a precision (即|Xk-Xk-1|/|Xk|