128 views (last 30 days)

Show older comments

I need to solve the following equation to w(t). I just need the real part solution. w(t) = .....

All other variables are later given in a table so I can calculate different solutions.

I´m not able to solve this in a m.file?

It would be great if someone can help me. please

Thanks a lot

Walter Roberson
on 31 Jan 2021

syms rho_w w(t) p_u d_0 D_0 l h_0 t kappa m_St m_w rho_L c_w A_Q v(t) h_d lambda zeta n

syms sum_of_zeta

eqn = (rho_w / 2) * w(t)^2 + p_u == (rho_w / 2) * (d_0 / D_0)^4 * w(t)^2 + ((l - h_0) / (l - h_0 + (d_0 / D_0)^2 * w(t) * t))^kappa + (rho_w / (m_St + m_w - rho_w * (pi / 4) * d_0^2 * w(t) * t)) * (rho_w * (pi / 4) * d_0^2 * w(t)^2 - (rho_L / 2) * c_w * A_Q * v(t)^2) * (h_0 + h_d - (d_0 / D_0)^2 * w(t) * t) + (rho_w / 2) * (d_0 / D_0)^4 * w(t)^2 * (lambda * ((h_0 - (d_0 / D_0)^2 * w(t) * t) / D_0) + sum_of_zeta) ;

syms W V

eqnW = subs(eqn, w(t), W)

solw = solve(eqnW, W)

char(eqnW)

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

Start Hunting!