Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

I have a physics question I need to program it by Maple , Can you help me to solve ? Its Problem 9.34 from Griffiths

I am working on an Isoperimetric problem.

I have two différential dL and dA. And I am seeking the function y for wich the two différentials are collinear.

So I want "dL/dA=constant"


I tried many expression, but I don't know how to express the "= constant"
I tried also (dL/dA)'=0 but no workable answer...


Thank's for your help !


I have the following problem

> restart:with(plots): Digits := 60:
> A:=42403:
> B:=269.95:
> C:=269.95:
> J1:=0.55:
> J2:=1.02:
> a:=196200:
> b:=14342220:
> c:=1589220:
> k:=157.08/0.0875:
> sys := {A-B*diff(x(t),t)=J1*diff(x(t),t,t)+C*(x(t)-y(t)),k*(J2*diff(y(t),t,t)-C*(x(t)-y(t)))=-(a+b*(diff(y(t),t)^2)/k^2+c*diff(y(t),t,t)),  D(x)(0) = 157.08, D(y)(0)=0, D(y)(10)=157.08,  D(x)(10) = 154.98}:
> dsn1 := dsolve(sys,numeric):
Error, (in dsolve/numeric/BVPSolve) matrix is singular

I need a help from someone who knows the GRTensor commands.

How do I write the last two terms of this equation of motion?How do I define Lagrangian matter, too?

My master said that I will do these calculations in the classic worksheet  maple18.

KGmn= (1/2) Tmn + (C0/2) Rabgl Tag Tbl gmn+ 3 C0 Rbrs(mTn)b Tsr – C0 Lm Rab Tab gmn + C0 Rab Tab Tmn+ 2 C0  Lm Rmrns Trs + C0ar(T(mr Tn)s)- C0sr(Tsr Tmn)

Lm =-(1/4) Fmn Fmn 

I'd like to view the Maple source code in .maple file. I am not able to find how to do this.

I tried commands in 136697-Can-One-Get-The-Source-Code-Back-From

which applies to .mla using the march command. But it did not work for .maple. It gives an error when trying it on .maple file.

I can open the .maple file in the Maple user interface OK using File->Open

but other than documentation, I see no Maple code.

I am looking now at Physics updates file, which can be downloaded as .maple file from 


Since most Maple code can be viewed, I assume one can also view the source code in the above.

Any idea how to do this?

Hey there,

I'm trying to numerically intergrate a function s(K,i,j) dK using runge kutta over a 2D grid of i,j values. Essentialy, performing the same sort of integral many different times for slightly different combinations of i and j. The function is more or less gaussian, and so the bulk of the result will come from the values of K around the peak of said gaussian. For some combinations of i and j, the function seems to have a singularity on the right edge of the gaussian peak, which causes my script to spit out an error, telling me the calculation cannot be performed further to the right past the singularity. Now, like i said before, the singularity is on the very edge of the gaussian and therefore I am perfectly happy to stop the integration before the singularity, because anything past it wont contribute very much to the result.

How can I use dsolve events to halt my integration just before hitting the singularity?

RK := (i, j) -> dsolve({diff(n(K), K) = K*s(K, a[i], b[j]), n(0) = 0}, numeric, method = rkf45)





I am trying to (again) extract the equations of the systems which have been done by me, but the equation extraction gave me wrong results.


For example;

The DC motor and the inverted pendulum equation must be 3. (third) order and only a differential equation.

The pneumatic motor and the inverted pendulum must be a differential equation (any order).

Also the hydraulic motor and the inverted pendulum must be a differential equation (any order).


I tried to use "isolate" command but it did not work for me in this case. How can I get the differential equations of the systems ?


I have added the MapleSim simulations and the differential equations that I had.


Please help me for this problem chain.

trying some examples from 209970-Exact-Solutions-For-PDE-And-Boundary--Initial-Conditions-2018

I noticed some answer change when adding assumptions which should not really cause the answer to change. Numerically, both answers seem to be the same. But one is series solution and the other with assumptions added uses integrals.

In wave pde on string where left end at x=0 and adding assumptions x>0,t>0 to pdsolve should produce the same form and it does on some problems but no on others.


#example 19, wave pde, on string left end fixed, right end moving, with source, no initial velocity
pde := diff(u(x, t), t$2) = 4*diff(u(x, t), x$2)+(1+t)*x;
bc := u(0,t)=0,u(Pi,t)=sin(t);
ic := u(x,0)=0,eval(diff(u(x,t),t),t=0)=0;

sol:=pdsolve([pde, bc,ic],u(x,t));  #1
sol:=pdsolve([pde, bc,ic],u(x,t)) assuming t>0; #2
sol:=pdsolve([pde, bc,ic],u(x,t)) assuming x>0; #3
sol:=pdsolve([pde, bc,ic],u(x,t)) assuming t>0,x>0;  #4

(1,2,3) above give same answer which is 

but (4) gives

Tried random points and numerically they are the same (set the sum to n=100 to evaluate).

my question is why would adding no assuming, or adding t>0 on its own or adding x>0 on its own give same answer, but adding  t>0,x>0 at same time give such a different looking soution.

I typically adding such assumptions to be safe, but here it seems to had some side effect I have not noticed before.

 "C:\Maple_updates\Physics+Updates.maple", 2018, December 22, 

    10:41 hours, version in the MapleCloud: 265, version installed in this computer: not installed

  Maple 2018.2, X86 64 WINDOWS, Nov 16 2018, Build ID 1362973


Hello. Is there a built-in function to determine the required sample size for different distributions in a population?

I try to solve ODE with conditions, but it give answer only without conditions:


Thank you.

Is there any facility to apply Finite Volume Method to Partial idifferential equation on MAPLE?
Any comand?

Any Code?

Modify the procedure that implements the secant and tangent routines in such a way that instead of the number of iterations in the beginning give the given accuracy E, with which the approximation is to be determined.

As a result, the procedure should return the last approximation along with the iteration number.




I need help.


I try to install MapleToolbox2018.0LinuxX64Installer.run I run it as root or with sudo command, but everytime I finish with below message:

Choose Maple 2018 Folder

Please specify the path to your existing Maple 2018 Installation.

 [/opt/maple2018]:  /opt/exp_soft/local/generic/maple/2018

Choose MATLAB Folder

Select your MATLAB installation.

 []: /opt/exp_soft/local/generic/matlab/R2016b

Warning: The installer has detected a previous version of the Maple Toolbox in
your MATLAB installation.
Press [Enter] to continue:
Warning: The directory
is not writable by the current user

I'm not sure what to do with it.




I have a vector A. I want to re-index vector A.

After re-indexing it, I will use the elements of the vector in new calculations.

For example:




But the above code doesn' t work.

I have two list in which I need to substitute variables from on elist with variables in other list in an expression.

rootElements := [[a, Integer, c, Integer, -k*a, b], [Integer, c, a, Integer, b, -k*a], [Integer, a, b, -k*a, c, Integer], [Integer, a, Integer, c, -k*a, b]];

list1 := [[a[maximum], 0, c[maximum], 0, (-k*a)[minimum], b[minimum]], [0, c[minimum], a[maximum], 0, b[maximum], (-k*a)[maximum]], [0, a[minimum], b[maximum], (-k*a)[maximum], c[maximum], 0], [0, a[minimum], 0, c[minimum], (-k*a)[minimum], b[minimum]]]

let the term be

term :=a*c/(k*a+b);

Terms := seq(applyrule([seq(rootElements[i][j] = list1[i][j], j = 1 .. nops(list1[i]))], term), i = 1 .. nops(list1));

my outputs  are

a[maximum]*c[maximum]/(k*a[maximum]+b[minimum]),# k*a should be minimum, becuae I substituted a earlier it is failing
a[minimum]*c[maximum]/(k*a[minimum]+b[maximum]),# k*a should be maximum, becuae I substituted a earlier it is failing

is there a way to substitute larger terms first[like -k*a first] and goes to smaller variable??

 I am using applyrule since algsubs cannot do for more than one substitution.Also, I would be glad if any other alternative is provided.
thank you

First 12 13 14 15 16 17 18 Last Page 14 of 248