Wrong answer in below code?

3 次查看（过去 30 天）

Raj Arora 2021-7-2

编辑： G A 2021-7-2

bw = 36;

angle = 60;

syms ds

tw = bw + 2*(ds / tand(angle));

nm = 0.025;

Area = 0.5*(bw + tw)*ds;

P = bw + 2*(ds/sind(angle));

S = 0.00001;

A = (((Area)^(5/3))/((P)^(2/3))) == ((200*nm)/(S^0.5))

height =(vpasolve(A,ds));

OUTPUT: - 49.95106 - 10.88106i

The output is wrong, the correct answer for the above code should be 10.19

Can anyone please tell me why I am not getting the actual answer as 10.19

Walter Roberson 2021-7-2

height = vpasolve(A, ds, 1);

Walter Roberson 2021-7-2

Most any positive initial value should work

G A 2021-7-2

编辑：G A 2021-7-2

Thanks, Walter - I understood postfactum that nm was interpreted as the initial value because it cannot be an argument.

在 Help Center 和 File Exchange 中查找有关 Newton-Raphson Method 的更多信息

Find the treasures in MATLAB Central and discover how the community can help you!

Translated by