@yangtheary There is no aspect of the problem statement or diagram that distinguishes the up-down dimension from the left-right dimension. Therefore, the outer lattice and the inner lattice are each squares. Therefore, the answer is a sum of two squares.

@casperyc I haven't looked at it yet, but when you say "failed" do you mean that the answer is wrong or that it takes too long to get?

The reverse order trick only worked becuase of a special form of the relations: Each one was the previous one essentially to a higher power (except that first one was independent of the rest). In general, if the relations are applied sequentially, the results will be wrong, i.e., not as simplified as if they were applied all together.

Please make this a separate Question thread. It is some trouble to dig up old threads.

@casperyc I haven't looked at it yet, but when you say "failed" do you mean that the answer is wrong or that it takes too long to get?

The reverse order trick only worked becuase of a special form of the relations: Each one was the previous one essentially to a higher power (except that first one was independent of the rest). In general, if the relations are applied sequentially, the results will be wrong, i.e., not as simplified as if they were applied all together.

Please make this a separate Question thread. It is some trouble to dig up old threads.

@Alejandro Jakubi I guess that it is too much to hope for at this point. But solve or allvalues will sometimes show a general pattern for a trancendental equation (with, _Zn~, _Bn~, _NNn~) that allows you to compute as many of an infinite sequence of roots as you want.

Something for future development: For the equation tan(x)=x, it could tell you that there is exactly one root in each interval (2*_Z1~-1)*Pi/2..(2*_Z1~+1)*Pi/2. Note that the notation for this already exists in RootOf: The second argument can be a range.

@Alejandro Jakubi I guess that it is too much to hope for at this point. But solve or allvalues will sometimes show a general pattern for a trancendental equation (with, _Zn~, _Bn~, _NNn~) that allows you to compute as many of an infinite sequence of roots as you want.

Something for future development: For the equation tan(x)=x, it could tell you that there is exactly one root in each interval (2*_Z1~-1)*Pi/2..(2*_Z1~+1)*Pi/2. Note that the notation for this already exists in RootOf: The second argument can be a range.

@Andriy It means the same thing as

nn:= proc(i,j,sigma)
     if type(i,integer) and type(j,integer) and type(sigma,integer) then
          return ap(i,j,sigma) . am(i,j,sigma)
     else
          return 'nn'(i,j,sigma) #return unevaluated
     end if
end proc;

@Andriy It means the same thing as

nn:= proc(i,j,sigma)
     if type(i,integer) and type(j,integer) and type(sigma,integer) then
          return ap(i,j,sigma) . am(i,j,sigma)
     else
          return 'nn'(i,j,sigma) #return unevaluated
     end if
end proc;

@Mathematix The answer will be a function for i(t). Are you sure that there is not some other equation in the problem? possibly an equation with a derivative of q? Now, I can use the above as the value of q(t) and get a solution for i(t), but I have a feeling that that is not what was intended.

Also, please clarify if the term is (1/C)*q(t) or 1/(C*q(t)). By default it would be the former.

With the information provided so far, the answer will be

ode:= L*D(i)(t) + (R[1]+R[2])*i(t) + 1/C*q(t) = U;

dsolve(
     eval(
          {ode, i(0)=0},
          [L= 2, R[1]= 2, R[2]= 8, C= 1/12,
           U= 4*exp(-t),
           q(t) = -3*exp(-3) + 4*exp(-2*t) + exp(-t)
          ]
     )
);

@Andriy Could you show an example of it not working?

@Andriy Could you show an example of it not working?

@Jimmy It's easy. Just say

d:= 0.28

before calling NonlinearFit. If you have specified d in the parameterranges, then you'll need to remove that also.

@Jimmy It's easy. Just say

d:= 0.28

before calling NonlinearFit. If you have specified d in the parameterranges, then you'll need to remove that also.

It is easy to solve for arbitrary q(t), but not possible to plot when q(t) is arbitrary. So, do you have more information about q(t)? We also need an initial condition for a plot, something like i(0) = ???.

The above solution is missing some roots. Another sublety of NextZero is that it helps to have a continuous function. The function in this case can easily be made continuous: The roots of tan(sqrt(z)) - sqrt(z) are equivalent to the roots of sin(sqrt(z)) - cos(sqrt(z))*sqrt(z).

Root:= table():  Root[0]:= 0:
for k to 20 do
     Root[k]:= RootFinding:-NextZero(
          z-> sin(sqrt(z))-cos(sqrt(z))*sqrt(z),
          Root[k-1],
          maxdistance= evalf(((k+1)*Pi)^2)
     )
od:
Root:= convert(Root, list);

The above solution is missing some roots. Another sublety of NextZero is that it helps to have a continuous function. The function in this case can easily be made continuous: The roots of tan(sqrt(z)) - sqrt(z) are equivalent to the roots of sin(sqrt(z)) - cos(sqrt(z))*sqrt(z).

Root:= table():  Root[0]:= 0:
for k to 20 do
     Root[k]:= RootFinding:-NextZero(
          z-> sin(sqrt(z))-cos(sqrt(z))*sqrt(z),
          Root[k-1],
          maxdistance= evalf(((k+1)*Pi)^2)
     )
od:
Root:= convert(Root, list);