```python
import sympy as sp
sp.init_printing(use_latex = True)
%matplotlib inline
```


```python
M_s, x, w, y, h, a, beta, t, nu = sp.symbols('M_s, x, w, y, h, a, beta, t, nu')
```


```python
H_x = M_s/(4*sp.pi) * (sp.log(((x+w)**2 + (y-h)**2)/((x+w)**2 + (y+h)**2)) - sp.log(((x-w)**2 + (y-h)**2)/((x-w)**2 + (y+h)**2)))
H_x
```


```python
H_y = M_s/(2*sp.pi) * (sp.atan((2*h*(x+w))/((x+w)**2 + y**2 - h**2)) - sp.atan((2*h*(x-w))/((x-w)**2 + y**2 - h**2)))
H_y
```


```python
H = sp.sqrt(H_x**2 + H_y**2)
H
```


```python
Hx = sp.diff(H, x)
HX = Hx.subs(y, x)
print(HX)
```

    (M_s**2*(2*(2*h*(-2*w - 2*x)*(w + x)/(-h**2 + x**2 + (w + x)**2)**2 + 2*h/(-h**2 + x**2 + (w + x)**2))/(4*h**2*(w + x)**2/(-h**2 + x**2 + (w + x)**2)**2 + 1) - 2*(2*h*(-w + x)*(2*w - 2*x)/(-h**2 + x**2 + (-w + x)**2)**2 + 2*h/(-h**2 + x**2 + (-w + x)**2))/(4*h**2*(-w + x)**2/(-h**2 + x**2 + (-w + x)**2)**2 + 1))*(-atan(2*h*(-w + x)/(-h**2 + x**2 + (-w + x)**2)) + atan(2*h*(w + x)/(-h**2 + x**2 + (w + x)**2)))/(8*pi**2) + M_s**2*(-2*((-2*w + 2*x)/((h + x)**2 + (-w + x)**2) + (2*w - 2*x)*((-h + x)**2 + (-w + x)**2)/((h + x)**2 + (-w + x)**2)**2)*((h + x)**2 + (-w + x)**2)/((-h + x)**2 + (-w + x)**2) + 2*((-2*w - 2*x)*((-h + x)**2 + (w + x)**2)/((h + x)**2 + (w + x)**2)**2 + (2*w + 2*x)/((h + x)**2 + (w + x)**2))*((h + x)**2 + (w + x)**2)/((-h + x)**2 + (w + x)**2))*(-log(((-h + x)**2 + (-w + x)**2)/((h + x)**2 + (-w + x)**2)) + log(((-h + x)**2 + (w + x)**2)/((h + x)**2 + (w + x)**2)))/(32*pi**2))/sqrt(M_s**2*(-log(((-h + x)**2 + (-w + x)**2)/((h + x)**2 + (-w + x)**2)) + log(((-h + x)**2 + (w + x)**2)/((h + x)**2 + (w + x)**2)))**2/(16*pi**2) + M_s**2*(-atan(2*h*(-w + x)/(-h**2 + x**2 + (-w + x)**2)) + atan(2*h*(w + x)/(-h**2 + x**2 + (w + x)**2)))**2/(4*pi**2))



```python
H1 = H.subs(x, 0)
H1
```


```python
H2 = H1.subs(y, 0)
H2
```


```python
H3 = H2.subs(h, 12.5e-6)
H3
```


```python
H4 = H3.subs(w, 25e-6)
H4
```


```python
H5 = H4.subs(M_s, 8.6e5)
H5
```


```python
H5.evalf()   #H5 = H0
```


```python
print(H)
```

    sqrt(M_s**2*(-log(((-h + y)**2 + (-w + x)**2)/((h + y)**2 + (-w + x)**2)) + log(((-h + y)**2 + (w + x)**2)/((h + y)**2 + (w + x)**2)))**2/(16*pi**2) + M_s**2*(-atan(2*h*(-w + x)/(-h**2 + y**2 + (-w + x)**2)) + atan(2*h*(w + x)/(-h**2 + y**2 + (w + x)**2)))**2/(4*pi**2))


15f43c405e2aaba353f6054116c995688f3c4a8e
ipynb
Jupyter Notebook
FHD/FHD H Calc.ipynb
ali-kin4/JupyterNotes-1
aeb944b43da4f731b46f25758f804110dd421fec
MIT
FHD/FHD H Calc.ipynb
ali-kin4/JupyterNotes-1
aeb944b43da4f731b46f25758f804110dd421fec
MIT
2021-05-17T07:45:07.000Z
2021-05-17T09:25:10.000Z
FHD/FHD H Calc.ipynb
ali-kin4/JupyterNotes-part1
aeb944b43da4f731b46f25758f804110dd421fec
MIT
Qwen/Qwen-72B
1. YES
2. YES