Thank you so much, for the answer. I have to get the solutions to plot
them. The solutions represent the critical speed of a shaft and there's no
other way. Other students are using other software for their project like
Maple or Mathematica, which they cause less problem. I was using python
Thank you so much, for the answer.
Il giorno martedì 11 giugno 2019 12:49:36 UTC+2, Shekhar Prasad Rajak ha
scritto:
>
> Also, using solveset we can get the solution pretty fast :
>
> ```
> solveset(omega_nf - 942.5 , J_u)
> {0.00235331614197391}
>
> ```
>
> Regards,
> Shekhar
>
> On Thursday,
Thank you so much, for the answer.
Il giorno sabato 8 giugno 2019 16:22:00 UTC+2, Chris Smith ha scritto:
>
> With check=False I get the solutions
> ```
> >>> solve(omega_nf-942.5,check=False)
> [-0.00027375075000, 0.00235331614197392]
> ```
> And bisect agrees with positive root
> ```
> >>>
Thank you so much. I'm using sympy 1.4. I'll try it.
Il giorno mercoledì 5 giugno 2019 23:20:27 UTC+2, Aaron Meurer ha scritto:
>
> What version of SymPy are you using? For me in 1.4, solve(omega_nf -
> 942.5, J_u, dict=True) returns [{J_u: 0.00235331614197392}]
>
> In general, solve() only
Also, using solveset we can get the solution pretty fast :
```
solveset(omega_nf - 942.5 , J_u)
{0.00235331614197391}
```
Regards,
Shekhar
On Thursday, 6 June 2019 02:50:27 UTC+5:30, Aaron Meurer wrote:
>
> What version of SymPy are you using? For me in 1.4, solve(omega_nf -
> 942.5, J_u,
With check=False I get the solutions
```
>>> solve(omega_nf-942.5,check=False)
[-0.00027375075000, 0.00235331614197392]
```
And bisect agrees with positive root
```
>>> nsolve(omega_nf-942.5,(.002,.003),solver='bisect')
0.00235331614197390
```
On Wednesday, June 5, 2019 at 4:20:27 PM UTC-5,
SymPy will find the solution eventually I think.
This is a polynomial of order 8 having only even powers so with a
substitution omega**2 -> x it's a quartic in x with complicated
symbolic coefficients. The general formula for a quartic is horrendous
and in this case your coefficients are already
Hi everybody, I'm tring to solve this equation without succes: omega_nf_eq
= 0
import sympy as sym
m,J_d,J_p,y,Y,omega,Omega,phi,Phi,z,Z,theta,Theta,k_yy,k_zz,k_phiphi,k_yphi,k_ztheta,k_thetatheta,plane_xy1,plane_xy2,plane_xz1,plane_xz2
=
What version of SymPy are you using? For me in 1.4, solve(omega_nf -
942.5, J_u, dict=True) returns [{J_u: 0.00235331614197392}]
In general, solve() only returns closed-form solutions, so if it
doesn't return a solution, it may just mean that it couldn't find one
in closed-form. If you know that
Probably I'm not using sympy correctly, here I have this equation :
omega_nf = sqrt(2)*sqrt(87791997.5351708 - 12563210.5479217*sqrt(-
0.00144380926150678*J_u + 48.4180817181289*(J_u + 0.00027375075)**2 -
3.58991000729848e-7)/(J_u + 0.00027375075) + 2392.05862861605/(J_u +
0.00027375075))/2
and
10 matches
Mail list logo
