Thomas Richard

Mr. Thomas Richard

3556 Reputation

13 Badges

14 years, 164 days
Technical professional in industry or government
Aachen, North Rhine-Westphalia, Germany

MaplePrimes Activity


These are answers submitted by Thomas Richard

Simply add this line:

plot(rhs(solution), t = 0 .. 5);

Your worksheet was saved with Maple 18 which is quite old.

The problem was fixed in a later version - not sure when exactly, but Maple 2023 gets it right. So the recommendation is to upgrade.

The following will return two formal solutions for nu, but I doubt they will be of much use:

Note that I have cleaned up parentheses a bit, and replaced U[w] by U__w (literal subscript, as opposed to indexing).

eta := (x,y,t) -> y/((nu*t*cos(alpha)+nu*x/U__w*sin(alpha))^(1/2));
psi := (x,y,t) -> U__w*(nu*t*cos(alpha)+nu*x/U__w*sin(alpha))^(1/2)*f(eta(x,y,t));

transf := [u(x,y,t) = diff(psi(x,y,t),y), v(x,y,t) = -diff(psi(x,y,t),x)];

eq1 := diff(u(x,y,t),x) + diff(v(x,y,t),y) = 0;
eq2 := diff(u(x,y,t),t) + u(x,y,t)*diff(u(x,y,t),x) + v(x,y,t)*diff(u(x,y,t),y) - nu*diff(u(x,y,t),y,y) = 0;
sys := [eq1,eq2]:

newsys := eval(sys,transf);
# newsys[1] is trivial: 0=0
neweq := simplify(newsys[2]);
solve(neweq,nu);

Does this help?

Perhaps it's more useful to run pdsolve on the original system sys.

There are several ways to obtain that; I find this one quite easy:

p := x^2-a*x-b*x+a*b;
                                           2
                   p := a b - a x - b x + x 

eval(p,[a=2,b=3]);
                           2          
                          x  - 5 x + 6

L := [[2,3],[-3,7],[9,10]]:

seq(eval(p, [a,b] =~ para), para in L);
           2             2              2            
          x  - 5 x + 6, x  - 4 x - 21, x  - 19 x + 90

See the help pages on eval and the tilde (elementwise operator) for more background information.

You can query an event and optionally install your own event handler like this:

NumericStatus(real_to_complex);
sqrt(-5);
NumericStatus(real_to_complex);

MyHandler:=proc(operator,operands,default_value)
   WARNING("You must assume the sqrt argument is positive."):
   return(default_value):
end proc:

NumericEventHandler(real_to_complex=MyHandler):

sqrt(-5);

For more details, please enter ?NumericEventHandler to open the help page.

If you additionally want sqrt(-5) to return the symbol undefined (as opposed to I*sqrt(5)), just load the RealDomain package. It is not perfect, however: it may not catch every situation where complex numbers arise.

That's right, they are not supported in Maple Flow.

For this particular example, you can get the result by entering 4*~L, however. The tilde operator will take care of mapping.

Hi,

do you really want to print the tests?

If working interactively in Maple is also an option, take a look at the PracticeSheet command. It's covering several topics (mainly from calculus), and will take care of the randomization.

This is covered by the expand command:

restart

eq := (cosh(zeta+d)-2*sinh(zeta+d))/sinh(zeta+d)

(cosh(zeta+d)-2*sinh(zeta+d))/sinh(zeta+d)

(1)

eeq := expand(eq)

cosh(zeta)*cosh(d)/(sinh(zeta)*cosh(d)+cosh(zeta)*sinh(d))+sinh(zeta)*sinh(d)/(sinh(zeta)*cosh(d)+cosh(zeta)*sinh(d))-2*sinh(zeta)*cosh(d)/(sinh(zeta)*cosh(d)+cosh(zeta)*sinh(d))-2*cosh(zeta)*sinh(d)/(sinh(zeta)*cosh(d)+cosh(zeta)*sinh(d))

(2)

normal(eeq)

(cosh(zeta)*cosh(d)+sinh(zeta)*sinh(d)-2*sinh(zeta)*cosh(d)-2*cosh(zeta)*sinh(d))/(sinh(zeta)*cosh(d)+cosh(zeta)*sinh(d))

(3)

simplify(eeq)

((cosh(d)-2*sinh(d))*cosh(zeta)-2*sinh(zeta)*(cosh(d)-(1/2)*sinh(d)))/(sinh(zeta)*cosh(d)+cosh(zeta)*sinh(d))

(4)

``

Download trigh-expansion.mw

One way to obtain these coefficients is as follows:
 

restart:

n:=3:

ls:=expand((3+2*x+3*x^2)^n)

27*x^6+54*x^5+117*x^4+116*x^3+117*x^2+54*x+27

(1)

rs:=sum(b[i]*x^i,i=0..2*n);

x^6*b[6]+x^5*b[5]+x^4*b[4]+x^3*b[3]+x^2*b[2]+x*b[1]+b[0]

(2)

sol:=solve(identity(ls=rs,x));

{b[0] = 27, b[1] = 54, b[2] = 117, b[3] = 116, b[4] = 117, b[5] = 54, b[6] = 27}

(3)

assign(sol); # optional step

b[4];

117

(4)

 


Please see ?solve,identity for more details.

Download solve-identity.mw

 

Yes, the free Maple Player is covering this. It is actually a bit more than a static viewer; please see the description. All mainstream platforms (Windows, Linux, macOS) are supported.

If I understand that correctly, you can prescribe some k<>0, and then search for the remaining variables:

k := 17:
eqn := (x+a)*(x+b)*(x+c)*(x+d)+k = 0;
sol := isolve(eqn);
eval(sol, _Z1=-4); # plug in any integer for _Z1
eval(eqn, %); # check eqn

Is that what you meant?

Maple is trying to tell you that the input framed in red does not make sense syntactically.

You need to insert a statement separator (colon or semicolon) at the end of the first assignment in the loop.

Next time, please upload a worksheet, as opposed to a screenshot.

[It seems that you had edited your question while I wrote my reply.]

I can only comment on the technical aspects. For purchasing the upgrade, please contact the sales department (or the reseller for India, I guess).

Like all toolboxes, the QCT requires the same Maple version it was made for. In particular, QCT 2023 will only work with Maple 2023 - and I vaguely guess that's what you are asking about...

In Maple 2023, you can tackle some of them with new (or newly documented) int option method=Pseudoelliptic (see ?updates,Maple2023,AdvancedMath).

I have not tried to simplify the results, though.

No, it cannot be suppressed. Don't you find it useful?

4 5 6 7 8 9 10 Last Page 6 of 45