```python
import sympy as sp
from sympy.physics import vector as vec
sp.init_printing()
vec.init_vprinting()
```


```python
R = vec.ReferenceFrame('R')
R
```




    R




```python
R??
```


```python
x, y, z = R[0], R[1], R[2]
# x, y, z, t = sp.symbols('x y z t', real=True)
t = sp.symbols('t', real=True)
c, eps, mu = sp.symbols('c varepsilon_0 mu_0')
c, eps, mu
```


```python
Ex= sp.Function('Ex')(x, t)
Ey= sp.Function('Ey')(x, t)
Ez= sp.Function('Ez')(x, t)
Bx= 0
By= sp.Function('By')(x, t)
Bz= sp.Function('Bz')(x, t)
```


```python
E = Ex*R.x + Ey*R.y + Ez*R.z
B = Bx*R.x + By*R.y + Bz*R.z
vec.curl(E, R)
```


```python
jx= sp.Function('jx')(x, t)
jy= sp.Function('jy')(x, t)
jz= sp.Function('jz')(x, t)
j = jx*R.x + jy*R.y +jz*R.z
```


```python
# FARADAY = sp.Eq(vec.curl(E, R), -B.diff(t,R))
FARADAY = B.diff(t, R)- vec.curl(E, R)
FARADAY
```


```python
AMPERE = vec.curl(B, R) - mu * j - eps * mu * E.diff(t, R)
AMPERE
```


```python
(FARADAY + c * AMPERE).simplify()
```

30d0f7fc19a655bcf3e52e656035a44093b56b0f
ipynb
Jupyter Notebook
docs/ipynb/fields.ipynb
StanczakDominik/PIC3
583262cff0edfaee48b9540505bcd68983ec53ec
BSD-3-Clause
2016-03-29T09:07:07.000Z
2021-09-27T07:59:17.000Z
docs/ipynb/fields.ipynb
StanczakDominik/PIC3
583262cff0edfaee48b9540505bcd68983ec53ec
BSD-3-Clause
2017-02-14T13:27:24.000Z
2017-03-10T19:53:03.000Z
docs/ipynb/fields.ipynb
StanczakDominik/PythonPIC
583262cff0edfaee48b9540505bcd68983ec53ec
BSD-3-Clause
2016-09-11T19:31:20.000Z
2021-01-11T03:26:02.000Z
Qwen/Qwen-72B
1. YES
2. YES