Unanswered Questions

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

I am considering the following PDE and I am getting an error, please suggest a better numerical method than the default one used in maple:

 

the PDE is:

u_{xx}u^3 - sin(xt)u_{tt} = u(x,t)

u(x, 0) = sin(x), (D[2](u))(x, 0) = cos(x), u(0, t) = cos(t), (D[1](u))(0, t) = sin(t)

Please suggest me a method that will also work for the following PDEs:

u^m* u_{xx} - sin(xt)u_{tt} = u^n

for m,n =0,1,2,3,... for the cases m=n and m not equal n

Here's the code:

 

pde := u(x, t)^3*(diff(u(x, t), x, x))-sin(x*t)*(diff(u(x, t), t, t)) = u(x, t);

u(x, t)^3*(diff(diff(u(x, t), x), x))-sin(x*t)*(diff(diff(u(x, t), t), t)) = u(x, t)

(1)

ibc := u(x, 0) = sin(x), (D[2](u))(x, 0) = cos(x), u(0, t) = cos(t), (D[1](u))(0, t) = sin(t);

u(x, 0) = sin(x), (D[2](u))(x, 0) = cos(x), u(0, t) = cos(t), (D[1](u))(0, t) = sin(t)

(2)

pds := pdsolve(pde, [ibc], numeric, time = t, range = 0 .. 1, spacestep = 0.1e-1)

module () local INFO; export plot, plot3d, animate, value, settings; option `Copyright (c) 2001 by Waterloo Maple Inc. All rights reserved.`; end module

(3)

pds:-plot3d(u(x, t), t = 0 .. 1, x = 0 .. 1, labels = [t, x, u(x, t)], labelfont = [times, bold, 20], axesfont = [times, bold, 16])

Error, (in pdsolve/numeric/plot3d) unable to compute solution for t>HFloat(0.0):
Newton iteration is not converging

 

``

 

Download nonlinear_hyperbolic_PDE.mw

Hello

anyone plz help me to find adomian decomposition method (ADM) maple code, also how to find the adomian polynomials separately using maple, i search it for a long time but did'nt succeeded,

I made reference to the show the enemy below and included part of the script an interesting part in the show. 

I did a mapleprimes search - only two returns none of which were what I was looking for.  I also did a google search which turned up nothing. 

Was my post deleted?

Good day everyone,

 

I want to construct groebner bases over rings of differential Operators.

Thus I used the following:

 

with(Ore_algebra);

with(Groebner); N := 3;

A := skew_algebra(diff = [D[1], x[1]], diff = [D[2], x[2]], diff = [D[3], x[3]], comm = i, alg_relations = i^2+1); T := MonomialOrder(A, tdeg(D[1], D[2], D[3]));

A["polynomial_indets"];  [returns    {D[1], D[2], D[3]}]

A["rational_indets"];       [returns    {i, x[1], x[2], x[3]}]

 

So far everything seems as it should be: The ring i wanted to define here is the third Ring of  Differential operators over the field of complex 'rational' functions and maple returns that indeed it will handle the D[i] as monomials and the rest as coefficients for them.

 

Hover, when i use the Skew product, the following happens:

skew_product(x[1],D[1],A)                    [returns D[1]x[1]+1}]

skew_product(x[1],D[1],A)                    [returns D[1]x[1]]

Both is wrong, but maple seems to 'know' this. It used the correct relation to in the first product, the only thing it didn't do was switch D[1] and x[1]. I think maple handles the result of the skew product as if it was a commutative product and always places the D[i] at the left but still 'knows' what the actual result is.

The actual results should have been x[1]D[1]+1 for the first and x[1]*D[1] for the second product.

 

In the second product, it seems like maple treated x[1]*D[1] as if the Elements were switched already.

 

What i want though is for maple to correctly display the skew products and return the products so that in every summand the D[i] are at the right side and their coefficients are at the left side (and if possible (i do not know how to do that yet) sort the result of a skew product in a way that displays every different Power product of the D[i] with their coefficients, so that i get (x[1]+x[2])*D[2] and not x[1]*D[2]+x[2]*D[2]).

 

Can anyone help me here?

 


                              
                     
                   
                 

I have been unable to successfully open a text file for writing in MAPLEPLAYER, whether it is requested in the startup code section or via a document tools button.  Is there a way to write to a file in MAPLEPLAYER?  It returns a "no write access" statement.

The write works perfectly in MAPLE 2015.1

MRB

Hi All,

 

I have o problem with simplify. A variable cp1r has been assumed to be positive. Why simplify still has csgn(cp1r) for it? Here is my code:

tmp := subs(cp1t(t)=cp1r, cp2t(t)=cp2r, Ca[2]);
1 / 2 2
----------- |-cp2r sin(x[1]) sin(x[7]) cp1r
2 2 |
cp1r cp2r |
\

2
+ 2 cp2r sin(x[1]) cos(x[1]) cos(x[7]) sin(x[7]) cp1r +

1 / 2 2 /
-------------- \cp2r cos(x[1]) cos(x[7]) sin(x[7]) \
(1/2)
/ 2\
2 \cp1r /
2 \\\
-2 cos(x[1]) cos(x[7]) sin(x[1]) + 2 sin(x[1]) cos(x[1])//|
|
|
/
assume(cp1r > 0, cp2r > 0);
simplify(tmp);
1 / / 3 3
---------- \sin(x[1]) sin(x[7]) \-cos(x[1]) cos(x[7])
2
cp1r cp1r

+ 2 cos(x[1]) cos(x[7]) cp1r csgn(cp1r) cp1r

2 3 \ \

- cp1r csgn(cp1r) cp1r + cos(x[1]) cos(x[7])/ csgn(cp1r)/

 

should csgn(cp1r) be simplified to 1 already? What is wrong with my script?

 

Thanks 

Everett

Maple does not cope with the following simple example:

with(geom3d):

point(A,0,0,0), point(B,1,0,0), point(C,2,0,0), point(E,2,1,0):

AreCoplanar(A,B,C,E);

           Error, (in geom3d:-plane) the points may not be AreCollinear

 

Should we interpret this behavior as a bug? I think I met with this yet 10-12 years ago, but unfortunately since then nothing has changed.

Just purchased Maple 2015 and playing with it for the first time.  I'm running the 64-bit version on Win 8.1.

Anytime the program generates a pop-up dialog box, the pop-up seems to get stuck behind the main program window.  I can't alt+tab to get to the pop-up window and I can't click on anything in the main program window because it's frozen while waiting for me to aknowledge the pop-up dialog box.  So I have to kill everything from the task manager and lose anything that's unsaved.

I'm talking about pop-ups for things like "Error, (in @@) invalid arguments"... I click on the link and see a brief flash while that pop-up quickly flies behind the main window.  Same thing happens when I click on a link that asks what web browser I want to use.

 

I can't be the only one that has this problem, so is there a fix or workaround for this somewhere?  It's not really useable like this.  Thanks!

Why am I getting different results in these two cases ?

 

I have the following in a text file.

a:="方程式";

When I read it in command-line maple, I am getting the following error.

on line 1, syntax error, non-ASCII character in unicode file 1:
Error, while reading `123456.txt`


Is it possible to input foreign characters in maple? (I dont want to enter directly in maple worksheet)

Currently I am running a maple engine using a java script and sending equations to the engine.

Thanks in advance.

i was to design a GUI using embedded components, i started as shown in the screenshot i uploaded but i had alittle problem trying to code the `evaluate at params` bottton i wanted the botton to return the evaluation of the mathematical expression at the value of the various parameters [delta,gamma,alpha,G,C] into the mathcontainer infront of it i.e evaluating at these points which can be change as new vualues are type in the textfeilds. such as

this is the code i used and it poping out an error message

use DocumentTools in

module()
local f,q,delta,gamma,alpha,G,C;

Do(f=%MathContainer0);
Do(delta=value(%MathContainer1));
Do(gamma=value(%MathContainer2));
Do(alpha=value(%MathContainer3));
Do(G=value(%MathContainer4));
Do(C=value(%MathContainer5));



q:=eval(f,value,delta,gamma,alpha,G,C);
Do(%MathContainer6=q);




end module;
pls i need some one to correct me

Hello everybody.

In the attached file, you find 6 equations. All of these parameters are known except "pd" and "qd". How can I find these two unknowns from 6 equations??? It should be pointed out that, "pd" and "qd" must contain "ud", "vd","wd" and "rd".

Thanks in advance.

D.mw

VIM.3.mwHi,How do I make a graph and error table for exact and by using this method i.e VIM 

First 221 222 223 224 225 226 227 Last Page 223 of 365