Carl Love

Carl Love

21383 Reputation

24 Badges

8 years, 287 days
Natick, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

Could you give an example? It is not clear to me why you presented the lists three times.

If they are lists, you should enclose them in []: [a,a,a,b,c]. Otherwise, they are called sequences, which are more difficult to work with than true lists.

Do you want the rules to apply to sublists also? That is, should [b,a,a,a,b,c,e] become [b,b,b,b,c,e,e]?

@Markiyan Hirnyk That is not an example, and you know it! I ask for a solution set returned by solve which contains a trivial equation. And that is what the Asker is referring to in his followup to your Answer.

@Markiyan Hirnyk That is not an example, and you know it! I ask for a solution set returned by solve which contains a trivial equation. And that is what the Asker is referring to in his followup to your Answer.

@Markiyan Hirnyk Please kindly show me an example of a solution set returned by solve, or a similar command, which contains a trivial equation and yet represents only a finite number of solutions.

@Markiyan Hirnyk Please kindly show me an example of a solution set returned by solve, or a similar command, which contains a trivial equation and yet represents only a finite number of solutions.

Here's how you can quickly check if any solution set (or list) S contains a trivial equation of the form z = z. Do

evalb(`*`(map(lhs-rhs, S)[]) = 0);

The result will be true if and only if S contains a trivial equation.

Here's how you can quickly check if any solution set (or list) S contains a trivial equation of the form z = z. Do

evalb(`*`(map(lhs-rhs, S)[]) = 0);

The result will be true if and only if S contains a trivial equation.

@fixedeggs But dsolve itself is evaluating the derivatives by using the right sides of dsys8. There is no "built-in interpolant" for this particular step. 

But if you really want to do it your way, you just need to add one equation and one initial condition to the system for each derivative that you want to evaluate. For example,

Sys:= {diff(y(t),t) = y(t)}: IC:= {y(0)=1}:
Sys1:= {diff(y(t), t) = w(t)}: IC1:= {w(0)=1}:
Sol:= dsolve(`union`(Sys,Sys1,IC,IC1), {y(t), w(t)}, numeric);
          proc(x_rkf45_dae)  ...  end;
Sol(1);
          [t = 1., w(t) = 2.7182818283577044,  y(t) = 2.7182818283577044]

@fixedeggs But dsolve itself is evaluating the derivatives by using the right sides of dsys8. There is no "built-in interpolant" for this particular step. 

But if you really want to do it your way, you just need to add one equation and one initial condition to the system for each derivative that you want to evaluate. For example,

Sys:= {diff(y(t),t) = y(t)}: IC:= {y(0)=1}:
Sys1:= {diff(y(t), t) = w(t)}: IC1:= {w(0)=1}:
Sol:= dsolve(`union`(Sys,Sys1,IC,IC1), {y(t), w(t)}, numeric);
          proc(x_rkf45_dae)  ...  end;
Sol(1);
          [t = 1., w(t) = 2.7182818283577044,  y(t) = 2.7182818283577044]

@Preben Alsholm This is a new usage of "so" as an introduction rather than as a continuation, which seems quite prevalent among academics. There is even a weekly radio show Le Show, hosted by Harry Shearer, that has a segment each week (about 2-3 minutes) devoted to playing clips of interviewees introducing their comments with "so". (Recorded episodes available at harryshearer.com/le-show/)

@LEETZ 

The boundary conditions from your original question will work just fine.

S := pdsolve({pde, bc});

@LEETZ 

The boundary conditions from your original question will work just fine.

S := pdsolve({pde, bc});

@digerdiga 

Sorry, I misunderstood what you meant by "parametrically", but I understand now. My fault: You didn't use the word incorrectly. You're right: The display method will not work to create the parametric plot of one dependent variable versus another.

@digerdiga 

Sorry, I misunderstood what you meant by "parametrically", but I understand now. My fault: You didn't use the word incorrectly. You're right: The display method will not work to create the parametric plot of one dependent variable versus another.

@digerdiga 

Any plots can be merged with plots:-display as long as they have the same dimensionality, i.e., all 2D or all 3D.

First 536 537 538 539 540 541 542 Last Page 538 of 593