Carl Love

Carl Love

28050 Reputation

25 Badges

12 years, 335 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

@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);

Unfortunately NextZero does not seem very reliable. In this case, it seems to not always find what is truly the next zero. I will try again by converting this to a continuous problem---trivial to do in this case.

Unfortunately NextZero does not seem very reliable. In this case, it seems to not always find what is truly the next zero. I will try again by converting this to a continuous problem---trivial to do in this case.

@abbeykabir Please start a new Question thread with this question. Also, please try to put some details of the question into the text itself, rather than leaving it all in the worksheet. I will look at your worksheet in the meantime.

Update: I took a look at it. It is a very interesting question. Please post it as a separate thread. Also, please attach a worksheet that I can execute. The current one does not give me access to the b; it merely displays its value.

@abbeykabir Please start a new Question thread with this question. Also, please try to put some details of the question into the text itself, rather than leaving it all in the worksheet. I will look at your worksheet in the meantime.

Update: I took a look at it. It is a very interesting question. Please post it as a separate thread. Also, please attach a worksheet that I can execute. The current one does not give me access to the b; it merely displays its value.

First 619 620 621 622 623 624 625 Last Page 621 of 709