I am currently working on a calculation in Mathematica that involves various parameters, including Nb, k, and Ns. I want to obtain an approximate form of the calculation by considering the following approximations: k << 1, Ns << 1 and Nb >> 1.
ClearAll["Global`*"]
S = 2 Ns + 1; B = 2 Nb + 1; A = 2*k*Ns + B; Cq = 2*Sqrt[Ns*(Ns + 1)];
vH01 = B;vH02 = S;
vH1[k_] = ((-1)^k (S - A) + Sqrt[(A + S)^2 - 4 k*(Cq)^2])/2
G[s_, x_] = (2^s)/((x + 1)^s - (x - 1)^s)
G1 = G[1/2, vH01];G2 = G[1/2, vH02];G3 = G[1/2, vH1[1]];G4 = G[1/2, vH1[2]];
L[s_, x_] = ((x + 1)^s + (x - 1)^s)/((x + 1)^s - (x - 1)^s)
L1 = L[1/2, vH01]; L2 = L[1/2, vH02]; L3 = L[1/2, vH1[1]];L4 = L[1/2, vH1[2]];
xp = Sqrt[0.5 ((A + S) +Sqrt[(A + S)^2 - 4*k*(Cq)^2])/(Sqrt[(A + S)^2 - 4*k*
(Cq)^2])]//FullSimplify
xm = Sqrt[0.5 ((A + S) -Sqrt[(A + S)^2 - 4*k*(Cq)^2])/(Sqrt[(A + S)^2 - 4*k*
(Cq)^2])]//FullSimplify
SA ={{xp, 0, xm, 0}, {0, xp, 0, -xm}, {xm, 0, xp,0}, {0, -xm, 0, xp}}
VAMid = {{L1, 0, 0, 0}, {0, L1, 0, 0}, {0, 0, L2, 0}, {0, 0, 0, L2}}
VA = (SA).(VAMid).(SA); Vdet = Det[VA]
These are the equations that I am trying to use:
This is the determinant equation that I am trying to simplify:
Could you please guide me on how to implement these approximations in Mathematica? Any suggestions, explanations, or code examples would be highly appreciated.