jakubi

1369 Reputation

12 Badges

19 years, 338 days

MaplePrimes Activity


These are replies submitted by jakubi

There is a problem with the condition for psi(alpha,.3361035810) as it involves Omega[1] and Omega[2], whose values contain the variable 't'. This is the origin of the error message arising at this point:

...
pds := pdsolve(PDE, IBC, numeric);

Error, (in pdsolve/numeric/process_IBCs) improper specification 
of initial/boundary condition: additional unknown(s) {t} were 
found in psi(alpha,.3361035810) = .7905694150e-2*2^(1/2)*Pi^(3/4)*
exp(-25*Pi)*(exp(-.5000000000e-6*(-t+1/2*arccosh(100*exp(-2*alpha)))^2+
I*(-100*t+50*arccosh(100*exp(-2*alpha))-1/2*(10000-exp(4*alpha))^(1/2)-
1/4*Pi))+exp(-.5000000000e-6*(t+1/2*arccosh(100*exp(-2*alpha)))^2-
I*(100*t+50*arccosh(100*exp(-2*alpha))-1/2*(1000...

There is a problem with the condition for psi(alpha,.3361035810) as it involves Omega[1] and Omega[2], whose values contain the variable 't'. This is the origin of the error message arising at this point:

...
pds := pdsolve(PDE, IBC, numeric);

Error, (in pdsolve/numeric/process_IBCs) improper specification 
of initial/boundary condition: additional unknown(s) {t} were 
found in psi(alpha,.3361035810) = .7905694150e-2*2^(1/2)*Pi^(3/4)*
exp(-25*Pi)*(exp(-.5000000000e-6*(-t+1/2*arccosh(100*exp(-2*alpha)))^2+
I*(-100*t+50*arccosh(100*exp(-2*alpha))-1/2*(10000-exp(4*alpha))^(1/2)-
1/4*Pi))+exp(-.5000000000e-6*(t+1/2*arccosh(100*exp(-2*alpha)))^2-
I*(100*t+50*arccosh(100*exp(-2*alpha))-1/2*(1000...

It is quite inconvenient trying to extract the input from your last post. Could you just write down your input lines?

It is quite inconvenient trying to extract the input from your last post. Could you just write down your input lines?

I mean, after IBC, and before the line:

module () local INFO; export plot, plot3d, animate, value, settings; option

some input is missing in your post above.

I mean, after IBC, and before the line:

module () local INFO; export plot, plot3d, animate, value, settings; option

some input is missing in your post above.

Inserting a 'tracelast' after

P1:=DEplot3d([X,Y,Z],[x(t),y(t),z(t)],t=0..300,[[x(0)=3.675,y(0)=0,z(0
)=-0.67]],numpoints=30000, thickness=1, linecolor=black,maxfun=0):

shows:

...
 #(`dsolve/numeric/SC/firststep`,24): f(neq,x,y,yp)
 f called with arguments: 3, 0., Array(1..3, [...], datatype = float[8], order
= C_order), Array(1..3, [...], datatype = float[8], order = C_order)
 #(f,9): YP[3] := .523598775598298*(-7*piecewise(0 <= Y[2],1/2*(2+1.41421356237
...
Error, (in f) unable to store '-1.06028752058655+0.231185221888386e-6*I' when datatype=float[8]

and

showstat(`dsolve/numeric/SC/firststep`,23..24);

`dsolve/numeric/SC/firststep` := proc(misc::array)
local W, V, neq, ndiff, npar, ncon, x, y, yp, f, Df, fx, fy, eca, eev, rh, Jm, Jp, rp, rv, i, j;
       ...
  23       if has(eval(f),'call_external') then
  24         f(neq,x,y,yp)
           else
             ...
           end if;
       ...
end proc

Now, the last library update for Maple 13.01, i.e. update.mla, is largely about 'dsolve/numeric'. So, my guess is that these code changes are related to this observed change in behavior from Maple 12 to 13.

Which is this complicated differential equation?

What about k and w?

What about k and w?

These (sub)threads provide some more comments about software companies, customers, EULAs and lawsuits:

old bugs

known bugs

When an upgrade isn't

Is it I(0,t) or Y(0,t)?

Is it I(0,t) or Y(0,t)?

Yes, the code that maps to:

sum((-1)^n, n = 0 .. infinity)

is

Typesetting:-mrow(`Typesetting:-munderover`(`Typesetting:-mo`
("&Sum;", mathvariant = "normal", fence = "false", separator = 
"false", stretchy = "true", symmetric = "false", largeop = 
"true", movablelimits = "true", accent = "false", lspace = 
"0.0em", rspace = "0.1666667em"), `Typesetting:-mrow`
(`Typesetting:-mi`("n", italic = "true", mathvariant = "italic"),
 `Typesetting:-mo`("=", mathvariant = "normal", fence = "false", 
separator = "false", stretchy = "false", symmetric = "false", 
largeop = "false", movablelimits = "false", accent = "false", 
lspace = "0.2777778em", rspace = "0.2777778em"), 
`Typesetting:-mn`("0", mathvariant = "normal")), 
`Typesetting:-mi`("∞", italic = "true", mathvariant = "italic"), 
accent = "false", accentunder = "false"), `Typesetting:-mspace`
(height = "0.0ex", width = "5.0", depth = "0.0ex", 
linebreak = "auto"), `Typesetting:-msup`(`Typesetting:-mfenced`
(`Typesetting:-mrow`(`Typesetting:-mo`("&uminus0;", 
mathvariant = "normal", fence = "false", separator = "false", 
stretchy = "false", symmetric = "false", largeop = "false", 
movablelimits = "false", accent = "false", lspace = 
"0.2222222em", rspace = "0.2222222em"), `Typesetting:-mn`("1", 
mathvariant = "normal")), mathvariant = "normal"), 
`Typesetting:-mi`("n", italic = "true", mathvariant = "italic"),
 superscriptshift = "0"));

which is rather different from the other code.

The underlying Typesetting code for this -1^n is:

Typesetting:-msup(`Typesetting:-mrow`(`Typesetting:-mi`(""), 
`Typesetting:-mrow`(`Typesetting:-mo`("&uminus0;", executable = 
"true", mathvariant = "normal", fence = "false", separator = 
"false", stretchy = "false", symmetric = "false", largeop = 
"false", movablelimits = "false", accent = "false", lspace = 
"0.2222222em", rspace = "0.2222222em"), `Typesetting:-mn`("1", 
executable = "true", mathvariant = "normal")), 
`Typesetting:-mi`("")), `Typesetting:-mi`("n", italic = "true",
 executable = "true", mathvariant = "italic"), 
superscriptshift = "0");

I.e. no brackets. You paste it in Standard 1D input, execute, and the output is like -1^n. Select the output, paste into 1D input and get again -1^n. On the other hand, including the code for the sum sign before:


Typesetting:-mrow(`Typesetting:-munderover`(`Typesetting:-mo`
("&Sum;", executable = "true", mathvariant = "normal", 
fence = "false", separator = "false", stretchy = "true", 
symmetric = "false", largeop = "true", movablelimits = "true", 
accent = "false", lspace = "0.0em", rspace = "0.1666667em"), 
`Typesetting:-mrow`(`Typesetting:-mi`("n", italic = "true", 
executable = "true", mathvariant = "italic"), 
`Typesetting:-mo`("=", executable = "true", mathvariant = 
"normal", fence = "false", separator = "false", stretchy = 
"false", symmetric = "false", largeop = "false", 
movablelimits = "false", accent = "false", 
lspace = "0.2777778em", rspace = "0.2777778em"), 
`Typesetting:-mn`("0", executable = "true", mathvariant = 
"normal")), `Typesetting:-mn`("∞", executable = "true", 
mathvariant = "normal"), accent = "false", accentunder = 
"false"), `Typesetting:-mi`(""), `Typesetting:-mrow`
(`Typesetting:-mi`(""), `Typesetting:-mrow`
(`Typesetting:-msup`(`Typesetting:-mrow`(`Typesetting:-mi`(""),
 `Typesetting:-mrow`(`Typesetting:-mo`("&uminus0;", executable 
= "true", mathvariant = "normal", fence = "false", separator 
= "false", stretchy = "false", symmetric = "false", largeop = 
"false", movablelimits = "false", accent = "false", lspace = 
"0.2222222em", rspace = "0.2222222em"), `Typesetting:-mn`("1", 
executable = "true", mathvariant = "normal")), 
`Typesetting:-mi`("")), `Typesetting:-mi`("n", italic = 
"true", executable = "true", mathvariant = "italic"), 
superscriptshift = "0"))));

It shows again no brackets on the output. But now, copy & pasted this output into 1D input produces brackets:

(&sum;)(-1)^n
First 38 39 40 41 42 43 44 Last Page 40 of 123