In[]:=
τθ=0.01;τq=0.015;ϵ=.001;
In[]:=
m=1+τqs+
1
2
2
(τqs)
;n=(1+τθs);θ0=
1
s
;
In[]:=
k1=
m
n
+
2
s
+
mϵ
n
-
-4mn
2
s
+
2
-m-n
2
s
-mϵ
n
2
;​​k2=
m
n
+
2
s
+
mϵ
n
+
-4mn
2
s
+
2
-m-n
2
s
-mϵ
n
2
;
In[]:=
h11=
-m+
2
k1
n
k1mϵ
;​​h12=
-m+
2
k2
n
k2mϵ
;​​h21=
-m+
2
k1
n-mϵ
mϵ
;​​h22=
-m+
2
k2
n-mϵ
mϵ
;
In[]:=
q1=
-m-mϵ+
1
2
n
m
n
+
2
s
+
mϵ
n
+
-4mn
2
s
+
2
-m-n
2
s
-mϵ
n
θ0
-4mn
2
s
+
2
(-m-n
2
s
-mϵ)
;​​q2=
m-n
2
s
+mϵ+
-4mn
2
s
+
2
(-m-n
2
s
-mϵ)
θ0
2
-4mn
2
s
+
2
(-m-n
2
s
-mϵ)
;