MaplePrimes Questions

For a month I did do not get email notfications.

I checked spam folders, settings and email filters.

Is it only me?

The attached sheet contains equations H1 to H6 and K1 to K3. What data do I need to modify to ensure that the values of H1 to H6 fall between 0 and 1, and K1 to K3 are negative? The parameter ranges are also given in the sheet. Is there a method to achieve this?

rouhg.mw

i did two case of this equation and odetest is worked good but in this case the odetest is not worked well anyone can determine what is mistake ?

F_P_Correct_case_three.mw

Dear Maple user I am facing error while running the codes  to plot the graph for two data sets .

I am attaching the files.

Error_in_Display_1.mw

Why does dsolve not call odetest by default before a solution is returned?

I mean, why do I have test each result separately. dsolve could have an odetest option (default=true).

In case of discrepancies dsolve could inform the user and suggest to call dsolve with odetest=false and run odetest separately to analyse the problem.

Set up this way, dsolve would never return potentially incorrect solutions that do not pass odetest.

HTR.mw

In above problem, Additionally How to  plot  heat transfer rate  Q versus L^2  for distinct porosity parmeters(Sh) , using  heat transfer rate formula, Q = (q*L)/(k*A*T[b])=theta'(1).

using  [Sh = 0.1, L^2 = 0.1, Nr =0 .1, Ha =0 .1, Pe = 0.1],  [Sh = 0.3, L^2 = 0.3, Nr = 0.1, Ha = 0.1, Pe =0 .1],   [Sh = 0.5, L^2 =0 .5, Nr =0 .1, Ha = 0.1, Pe =0 .1].

How can I find Mean, Median, Quartiles, Variance, StandardDeviation of data in this table

I use Mathamatica and get the result

Clear["Global`*"]
boundaries = Range[0, 10, 5/2];
frequencies = {18, 11, 13, 6};
binMeans = Mean /@ Partition[boundaries, 2, 1];
weighted = WeightedData[binMeans, frequencies];
weightedHist = HistogramDistribution[weighted, {5/2}];
Through[{Mean, Median, Quartiles, Variance, StandardDeviation}[
  weightedHist]]

I have number of worksheets open. I set each to use its own engine.

If I have one worksheet running a long loop that takes hrs to complete, or even if I have cmaple.exe running long script in command line, and then at same time I open new worksheet and type

             Physics:-Version(latest)

Then I always get this error

This is described here

It seems because Maple was busy running another worksheet. I never see this error if I am not running any other thing in Maple at the time. ( I actually try not to update Physics while running something else in Maple, but sometimes I do not notice that I have something else running).

When I wait for my other worksheet to finish, then close all of Maple and reopen Maple, I see the latest version is installed, at least this is what Physics:-Version() now says.

My question is:  Is it safe to do Physics:-Version(latest) even when Maple is busy running  computation in other separate worksheet or even running cmaple.exe from windows command line? Is there any concurrency/locking issue on the Physics library being updated while it is being open by another worksheet or another server.exe running in background?

All this is on windows 10 pro.

I trying to simplify expressions for lines so no higher order terms. factor and op seperate out what I need but how do I select the one with the variables in this case x,y. I cant depent on this always been the last one returned from the op command.

restart

 

vars=[x,y]

vars = [x, y]

(1)

k:=17

17

(2)

expr:=-2*sqrt(118)*(((-4*x + y + 51/32)*sqrt(k) + (k*x)/4 - (51*y)/4 + 153/32)*sqrt(-4012 + 1003*sqrt(k)) + ((x + 4*y)*sqrt(k) - (85*x)/4 - (17*y)/4)*sqrt(4012 + 1003*sqrt(k)))*k^(1/4)/(17051*(-1 + sqrt(k)))

-2*118^(1/2)*(((-4*x+y+51/32)*17^(1/2)+(17/4)*x-(51/4)*y+153/32)*(-4012+1003*17^(1/2))^(1/2)+((x+4*y)*17^(1/2)-(85/4)*x-(17/4)*y)*(4012+1003*17^(1/2))^(1/2))*17^(1/4)/(-17051+17051*17^(1/2))

(3)

indets(expr)

{x, y}

(4)

factor(expr)

(1/64192)*(-4012+1003*17^(1/2))^(1/2)*118^(1/2)*17^(1/4)*(17^(1/2)+5)*(40*x-24*y-3)

(5)

op(factor(expr))

1/64192, (-4012+1003*17^(1/2))^(1/2), 118^(1/2), 17^(1/4), 17^(1/2)+5, 40*x-24*y-3

(6)

select(has,vars, factor(expr))

select(has,vars, op(factor(expr)))

Error, invalid input: has expects 2 arguments, but received 7

 
 

 

Download 2024-09-11_Has_Select_Question.mw

I get my on results but the results are not the same please help me if i did any mistake in my code

 

symmetry_PDESYS_3_time_fraction[1].mw

This code is working for function f1 but not for f2
f2 := (x,y)->9*x^2-24*x*y+16*y^2+10*x-70*y + 175;
Why this code is not working for f2 ?
unprotect(D);
f1:= (x, y) -> 3*x^2 - 3*y*x + 6*y^2 - 6*x + 7*y - 9;
coeffs(f(x, y));
A, B, C, D, E, F := %;
theta := 1/2*arctan(B/(A - C));
solve({-2*A*xc - B*yc = D, -B*xc - 2*C*yc = E});
assign(%);
x := xcan*cos(theta) - ycan*sin(theta) + xc;
y := xcan*sin(theta) + ycan*cos(theta) + yc;
Eq := simplify(expand(f1(x, y)));
xcan^2/simplify(sqrt(-tcoeff(Eq)/coeff(Eq, xcan^2)))^`2` + ycan^2/simplify(sqrt(-tcoeff(Eq)/coeff(Eq, ycan^2)))^`2` = 1;

Thank you

In Maple 2024, using commands like 
with(LinearAlgebra):
M := Matrix([[1, 2, 3, 10], [4, 15, 6, 7], [0, 0, 0, 0]]);
ReducedRowEchelonForm(M);

we see output of the 3x4 matrices fine at first, but when I and my students try to print or save to a PDF, ellipsis (three dots) sometimes show up in the bottom row. I checked the settings for elision and they are as usual. This is only when trying to print---the ellipsis show up in both the Maple file and the printed copy. We've tested on both Mac and PC. This did not happen in Maple 2023.

I have side-by-side comparisons in an image below: the Maple looking fine before we try to print or save as a PDF, what happens in the Maple file after with the bottom row, and the elision settings. Any help would be greatly appreciated as we would like the full matrices to show up in printing or a PDF we save to.

I'm trying to solve a system of coupled differential equations numerically, but I'm getting the following error

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

The error occurs at the dsolve step, despite trying to ensure that all equations and conditions are in the correct form (sets/lists).

Could someone help me identify what I'm missing here?

Thanks in advance!

L := 200; K := 99; kappa := 1; omegaD := 1; beta := 1; delta := 0.5e-1; j := 2; tmax := 3000; h := L/(K+1); nsp := [`$`(-(1/2)*L+h*i, i = 0 .. L/h)]; km := nops(nsp); omegaD2 := h^2*omegaD^2; deltaHat := h*delta; a := 2; var := seq(x[i](t), i = 1 .. km); initialPositions := seq(x = a*sin(j*h*Pi*nsp[i]/L), i = 1 .. km); initialVelocities := seq((D(x[i]))(0) = 0, i = 2 .. km-1)

boundaryConditions := [x[1](t) = 0, x[km](t) = 0]

equations := seq(diff(x[n](t), t, t)-kappa*(x[n+1](t)-2*x[n](t)+x[n-1](t))+deltaHat*(diff(x[n](t), t))-omegaD2^2*(x[n](t)-beta*x[n](t)^3) = 0, n = 2 .. km-1)

sol := dsolve({equations, boundaryConditions, initialPositions, initialVelocities}, var, numeric, method = rkf45, range = 0 .. tmax)

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

 

NULL

Download dsolve_error.mw

Why Maple returns y(x)=0 as solution to this ode when even odetest do not validate the solution? Solution satisfies the ode but not the given BC.

Is this a known bug?

I've updated the worksheet now to include few more examples, all using same ode but with different BC's. In all 4 examples, Maple gives solutions when there should not be solution and in one case gives solution which can not even be used for odetest. These are random BC's generated just for testing.
 

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1804. The version installed in this computer is 1802 created 2024, September 3, 11:35 hours Pacific Time, found in the directory C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib\`

libname;

"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"

restart;

#EXAMPLE 1

ode:=diff(diff(y(x), x), x) + 2*diff(y(x), x)/x - y(x) = 0;
IC:=y(-1)=0,D(y)(-infinity)=-1/exp(1);

diff(diff(y(x), x), x)+2*(diff(y(x), x))/x-y(x) = 0

y(-1) = 0, (D(y))(-infinity) = -1/exp(1)

maple_sol:=dsolve([ode,IC]);

y(x) = 0

odetest(maple_sol,[ode,IC])

[0, 0, -exp(-1)]

restart;

#EXAMPLE 2

ode:=diff(diff(y(x), x), x) + 2*diff(y(x), x)/x - y(x) = 0;
IC:=y(infinity)=0,D(y)(-infinity)=-1/exp(1);

diff(diff(y(x), x), x)+2*(diff(y(x), x))/x-y(x) = 0

y(infinity) = 0, (D(y))(-infinity) = -1/exp(1)

maple_sol:=dsolve([ode,IC]);

y(x) = 0

odetest(maple_sol,[ode,IC])

[0, 0, -exp(-1)]

restart

#EXAMPLE 3

ode:=diff(diff(y(x), x), x) + 2*diff(y(x), x)/x - y(x) = 0;
IC:=y(-infinity)=0,D(y)(-infinity)=-1/exp(1);

diff(diff(y(x), x), x)+2*(diff(y(x), x))/x-y(x) = 0

y(-infinity) = 0, (D(y))(-infinity) = -1/exp(1)

maple_sol:=dsolve([ode,IC]);

y(x) = signum((sinh(x)+cosh(x))/x)*infinity

odetest(maple_sol,[ode,IC])

Error, (in signum) signum is not differentiable at 0

#EXAMPLE  4

restart;

ode:=diff(diff(y(x), x), x) + 2*diff(y(x), x)/x - y(x) = 0;
IC:=y(-infinity)=0,D(y)(infinity)=-1/exp(1);

diff(diff(y(x), x), x)+2*(diff(y(x), x))/x-y(x) = 0

y(-infinity) = 0, (D(y))(infinity) = -1/exp(1)

maple_sol:=dsolve([ode,IC]);

y(x) = 0

odetest(maple_sol,[ode,IC])

[0, 0, -exp(-1)]

 


 

Download wrong_sol_dsolve_sept_10_2024.mw

Update sept 13, 2024

Here is one more example I found. In this ode, the IC given leads to division by zero in the solution. Yet Maple for some reason removes the offending part of the solution with that constant, and returns the rest. 

Is one allowed to do this? odetest also do not validate the solution. There should not be solution returned in this case, since it is not possible to find values for the constants of integration given these initial conditions. Right?

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

restart;

ode:=x^2*diff(y(x),x$2)+x*diff(y(x),x)-9*y(x)=0;
IC:=y(0)=1,D(y)(0)=0;
maple_sol_no_IC:=dsolve(ode);

x^2*(diff(diff(y(x), x), x))+x*(diff(y(x), x))-9*y(x) = 0

y(0) = 1, (D(y))(0) = 0

y(x) = c__1*x^3+c__2/x^3

maple_sol_with_IC:=dsolve([ode,IC]); #maple just removed the c2/x^3 part. Since at x=0 gives problem?

y(x) = c__1*x^3

odetest(maple_sol_with_IC,[ode,IC]); #shows solution does not verify

[0, 1, 0]

 


 

Download another_example_solution_given_when_none_exist.mw

 

 

1 2 3 4 5 6 7 Last Page 2 of 2356