Hyperbolic iterated functions computations
Hyperbolic iterated functions computations
PrintDynamic[f_Symbol,z_,n_,Ω_,terms_]:= Module{}, Clear[DD];DD[f,n,Ω,1]=Derivative[1][f][Ω]^n;DD[f,m_,Ω,1]:=DD[f,n,Ω,1]/.nm; Forderv=2,derv≤terms,derv++, Printf[[z]]; Printf[[z]]//.Derivative[derv][][z]; Printf[[z]]//.Derivative[derv][][z]//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]; Printf[[z]]//.Derivative[derv][][z]//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]/.zΩ; Printf[[z]]//.Derivative[derv][][z]//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]/.zΩ/.0/.Derivative[1][f][Ω]ℒ; Printf[[z]]//.Derivative[derv][][z]//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]/.zΩ/.0/.Derivative[1][f][Ω]ℒ/.-1(-); Printf[[z]]//.Derivative[derv][][z]//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]/.zΩ/.0 /.Derivative[1][f][Ω]ℒ/.-1(-)/.ℒDerivative[1][f][Ω]; DD[f,n,Ω,derv]=f[[z]]//.Derivative[derv][][z]//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]/.zΩ/.0/.Derivative[1][f][Ω]ℒ/.-1(-)/.ℒDerivative[1][f][Ω]; sderv=ToString[IntegerPart[derv]]; ToExpression["DD[f,m_, Ω, derv]:= DD[f,n, Ω,"<>sderv<>" ] /. n m"]; ;;
∂
{z,derv}
n-1
f
∂
{z,derv}
n-1
f
-1+n
f
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
b_n+c_
ℒ
b+c-1
ℒ
b-1
ℒ
bn
ℒ
n
ℒ
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
b_n+c_
ℒ
b+c-1
ℒ
b-1
ℒ
bn
ℒ
n
ℒ
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
b_n+c_
ℒ
b+c-1
ℒ
b-1
ℒ
bn
ℒ
n
ℒ
General::spell1:Possible spelling error: new symbol name "sderv is similar to existing symbol derv.
Clear[DD,Dyn,Dynamic];
terms=2;DD[f,n,Ω,0]=Ω;DD[f,n,Ω,1]=Derivative[1][f][Ω]^n;DD[f,m_,Ω,1]:=DD[f,n,Ω,1]/.nm;If[Abs[Derivative[1][f][Ω]]==1, cyclotomic[f]:=True, cyclotomic[f]:=False];
General::spell1:Possible spelling error: new symbol name "cyclotomic is similar to existing symbol Cyclotomic.
Forderv=2,derv≤terms,derv++, term=f[[z]]; term=term//.Derivative[derv][][z]; term=term//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]; term=term//.zΩ; term=term//.[Ω]Ω; term=term//.Derivative[1][f][Ω]ℒ; Print[term]; term=term/.0; DD[f,n,Ω,derv]=term//.; DD[f,n,Ω,derv]=Collect[DD[f,n,Ω,derv],]; sderv=ToString[IntegerPart[derv]]; ToExpression["DD["<>sf<>",m_,"<>sΩ<>", derv]:= DD[" <>sf<>","<>sn<>", "<>sΩ<>","<>sderv <>" ] /. "<>sn<>" m " ];; Dyn[f,n,Ω]:=Sum[DD[f,n,Ω,cnt]cnt!,{cnt,0,terms}]/.ℒDerivative[1][f][Ω];
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
m_
f
b_n+c_
ℒ
b*(n-1)
∑
i=(b-1)*(n-1)
i
ℒ
n
ℒ
cnt
(z-Ω)
Forderv=2,derv≤terms,derv++,term=f[[z]]; PrintHoldFormf[[z]],"=",term; term=term//.Derivative[derv][][z]; Print[term]; term=term//.Derivative[vcnt_][][z]DD[f,-1+n,Ω,vcnt]; Print[term]; term=term//.zΩ; Print[term]; term=term//.[Ω]Ω; Print[term]; term=term//.Derivative[1][f][Ω]ℒ; Print[term]; term=term/.0; Print[term]; DD[f,n,Ω,derv]=term//.; DD[f,n,Ω,derv]=Collect[DD[f,n,Ω,derv],];;
∂
{z,derv}
n-1
f
∂
{z,derv}
n-1
f
-1+n
f
-1+n
f
m_
f
b_n+c_
ℒ
b*(n-1)
∑
i=(b-1)*(n-1)
i
ℒ
n
ℒ
∂
{z,derv}
n-1
f
2
′
()
-1+n
f
-1+n
f
′′
()
-1+n
f
2[z]+2[z]
-1+n
f
2
′
()
-1+n
f
-1+2n
2
-2+2n
Ω
-1+n
f
-1+2n
2
-2+2n
Ω
-1+n
f
2Ω+
-1+2n
2
-2+2n
Ω
ℒ+
-1+2n
2
-2+2n
Ω
-1+2n
2
-2+2n
Ω
cyclotomic[f]
cyclotomic[f]
f
f
PrintDynamic[f,z,n,Ω,3];
2
′
()
-1+n
f
′′
f
-1+n
f
′
f
-1+n
f
′′
()
-1+n
f
[[z]]+[z][[z]]
′
f
-1+n
f
2
′
()
-1+n
f
′′
f
-1+n
f
[[z]]+[Ω][[z]]
′
f
-1+n
f
-2+2n
′
f
′′
f
-1+n
f
[[Ω]]+[Ω][[Ω]]
′
f
-1+n
f
-2+2n
′
f
′′
f
-1+n
f
-2+2n
ℒ
′′
f
-1+n
f
(-+)[[Ω]]
n
ℒ
2n
ℒ
′′
f
-1+n
f
(-1+ℒ)ℒ
(-[Ω]+[Ω])[[Ω]]
n
′
f
2n
′
f
′′
f
-1+n
f
(-1+[Ω])[Ω]
′
f
′
f
3[z][[z]][z]+[z][[z]]+[[z]][z]
′
()
-1+n
f
′′
f
-1+n
f
′′
()
-1+n
f
3
′
()
-1+n
f
(3)
f
-1+n
f
′
f
-1+n
f
(3)
()
-1+n
f
[[z]]+3[z][[z]][z]+[z][[z]]
′
f
-1+n
f
′
()
-1+n
f
′′
f
-1+n
f
′′
()
-1+n
f
3
′
()
-1+n
f
(3)
f
-1+n
f
[[z]]++[Ω][[z]]
′
f
-1+n
f
3[Ω]-[Ω]+[Ω][[Ω]][[z]]
-2+n
′
f
-1+n
′
f
2(-1+n)
′
f
′′
f
-2+n
f
′′
f
-1+n
f
-1+[Ω]
′
f
-3+3n
′
f
(3)
f
-1+n
f
[[Ω]]++[Ω][[Ω]]
′
f
-1+n
f
3[Ω]-[Ω]+[Ω][[Ω]][[Ω]]
-2+n
′
f
-1+n
′
f
2(-1+n)
′
f
′′
f
-2+n
f
′′
f
-1+n
f
-1+[Ω]
′
f
-3+3n
′
f
(3)
f
-1+n
f
The Stuff
The Stuff