```python
import sympy as sym
import matplotlib.pyplot as plt
import stanpy as stp

x = sym.Symbol("x")
l = 4  # m
E = 2.1e8
b, ha, hb = 0.2, 0.3, 0.4  # m
hx = ha + (hb - ha) / l * x  # m
cs = stp.cs(b=b, h=hx)

s = {"E": E, "cs": cs, "q": 10, "l": l, "bc_i": {"w": 0, "M": 0, "H": 0}, "bc_k": {"w": 0, "phi": 0}}

print(cs)
print("I_y(0) = ", cs["I_y"](0))
print("I_y(l) = ", cs["I_y"](l))

offset = 0.2
fig, ax = plt.subplots()
stp.plot_cs(ax, b, hx.subs(x, 0))
stp.plot_cs(ax, b, hx.subs(x, l), dy=0.3)

ax.set_xlim(-offset, b + offset)
ax.set_ylim(-offset, float(hx.subs(x, l)) + offset)
ax.grid(linestyle=":")
ax.axis('equal')

plt.show()

fig, ax = plt.subplots()
stp.plot_system(ax, s, render=True, facecolor="gray", alpha=0.5, render_scale=0.4)
ax.grid(linestyle=":")
ax.axis('equal')

plt.show()
```


```python

```


```python

```


```python

```


```python

```


```python

```

1470034f05ac745ac832272e1be110ea2fccaf80
ipynb
Jupyter Notebook
tests/docs/user_guide/jupyter/non_constant/ex01.ipynb
DavidNaizheZhou/stanpy
257072bd52154e9e4d68be957fd12eee1ad3dc56
MIT
tests/docs/user_guide/jupyter/non_constant/ex01.ipynb
DavidNaizheZhou/stanpy
257072bd52154e9e4d68be957fd12eee1ad3dc56
MIT
tests/docs/user_guide/jupyter/non_constant/ex01.ipynb
DavidNaizheZhou/stanpy
257072bd52154e9e4d68be957fd12eee1ad3dc56
MIT
2022-01-30T17:25:38.000Z
2022-01-30T17:25:38.000Z
Qwen/Qwen-72B
1. YES
2. YES