The best and simplest way to compute a quaternionic derivative is as follows. We represent a quaternion argument p = x+y⋅i+z⋅j+u⋅k (i,j,k are basic quaternion units; "⋅" is the quaternion multiplication) and a quaternion-differentiable (holomorphic) function of that argument ψ(p) = ψ
1(x,y,z,u)+ψ
2(x,y,z,u)⋅i+ψ
3(x,y,z,u)⋅j+ψ
4(x,y,z,u)⋅k in the Cayley–Dickson doubling form: p = a+b⋅j, where a´= x+y⋅i ; b = z+u⋅i and ψ(p)=ψ(a,b)=Φ
1(a,b)+Φ
2(a,b)⋅j, where Φ
1(a,b)=ψ
1(a,b)+ψ
2(a,b)⋅i and Φ
2(a,b)=ψ
3(a,b)+ψ
4(a,b)⋅i. Each expression for ψ(p) is initially to be obtained from a complex function of the same kind by means of the direct replacement of a complex variable with a quaternion variable in the expression for the complex function. For example, ψ(p)=p
-1. Just as a complex- holomorphic function satisfies Cauchy-Riemann's equations in complex analysis, a quaternion- holomorphic function satisfies the following quaternionic generalization of Cauchy-Riemann's equations:
(1) ∂
aΦ
1 = ∂
b×Φ
2×, (2) ∂
aΦ
2 = - ∂
b×Φ
1×,
(3) ∂
aΦ
1 = ∂
bΦ
2, (4) ∂
a×Φ
2 = - ∂
b×Φ
1
after doing a = a× = x,
where the complex conjugation is denoted by × and the
partial differentiation with respect to some variable s is denoted by ∂
s. For example, by ∂
b×Φ
2× is denoted the partial derivative of the complex conjugate of a function Φ
2 with respect to the complex conjugate of a complex variable b. Firstly, we compute the partial derivatives of functions Φ
1, Φ
2, Φ
1×, Φ
2× (with respect to variables a, b, a
×, b
×); secondly, we put a = a
× =x in the computed expressions of partial derivatives; and thirdly, we check whether equations (1) - (4) hold. One of the formulae to compute the first quaternionic derivative of the quaternion-holomorphic function is the following:
ψ(p)(1) = (∂aΦ1 + ∂a×Φ1) + (∂aΦ2 + ∂a×Φ2)⋅j .
Higher derivatives of quaternion-holomorphic functions can be computed analogically and they are holomorphic like the first derivative.
For details and examples I refer to
http://vixra.org/abs/1609.0006