Commit 17fc905e authored by Blais, Chris's avatar Blais, Chris
Browse files

from cleanup, moved files to appropriate folders

parent c6cc1688
Loading
Loading
Loading
Loading
+450 −0

File added.

Preview size limit exceeded, changes collapsed.

+54 −0
Original line number Diff line number Diff line
%% Cell type:code id:ccffd49a tags:

``` python
### minimum working example
import sympy as sy
m2, m3, m4, m5, m6, m7, m8 = sy.symbols("m2 m3 m4 m5 m6 m7 m8")
x2, x3, x4, x5, x6, x7, x8 = sy.symbols("x2 x3 x4 x5 x6 x7 x8")
b1, b1p, c1, c1p, e2, h2, g2, j2 = sy.symbols("b1 b1p c1 c1p e2 h2 g2 j2")
dp35, dp78 = sy.symbols("dp35 dp78")

eqs = []
eqs.append(m2 - m3 + m8)
eqs.append(m2*x2 - m3*x3 + m8*x8)
eqs.append(m3 - m4 - m5)
eqs.append(m3*x3 - m4*x4 - m5*x5)
eqs.append(b1*m4*x4 - b1p*(x3 - x4))
eqs.append(m5 - m6 - m7)
eqs.append(-m7 + m8)
eqs.append(c1*dp35 - c1p*(m3+m5)**2)
eqs.append(j2*dp78 - e2*m7**2 - g2 - h2*m7)
eqs.append(dp78-dp35)

x0_measured = {
    'm1': 50, 'm2': 50,'m4': 10, 'm6': 40,
    'x5': 50, 'x6': 50,'x7': 50, 'x8': 50,
    'b1': 1, 'b1p': 5,
    'c1': 338,'c1p': 1,
    'e2': 50, 'g2': 50, 'h2': 50, 'j2': 421
    }

eq_sys = []
for eq in eqs:
    eq_sys.append(sy.Eq(eq.subs(x0_measured), 0))

sy.solve(eq_sys, domain=sy.Reals)
```

%% Output

    [{dp35: 98109025/7235208,
      dp78: 98109025/7235208,
      m3: 147785/3804,
      m5: 109745/3804,
      m7: -42415/3804,
      m8: -42415/3804,
      x2: 3462010/81063,
      x3: 1097450/27021,
      x4: 1097450/81063},
     {dp35: 50, dp78: 50, m3: 70, m5: 60, m7: 20, m8: 20, x2: 43, x3: 45, x4: 15}]

%% Cell type:code id:d5b56247 tags:

``` python
```
+264 −0

File added.

Preview size limit exceeded, changes collapsed.