## Collect wrt to fraction...

Dear All

How can we collect coefficient wrt certain differential ration in an expression?

See for detail:

 >
 >
 (1)
 >
 (2)
 >
 (3)

How ca extract coefficients of fraction  in (3) ????

Regards

## Characteristic polynomial from solution...

solution of system of differential equations

is {a(t) = 1, b(t) = exp(t), c(t) = 0}

1 = exp(t), since t = 0, is characteristic equation = (y-0)?
0 = exp(t), how to find this characteristic equation when division is zero?
exp(t) = exp(t), is characteristic equation = (y-t)?

## Collecting coefficients in differential expression...

Dear all

I have problem related to collection of coefficient of differtials in differential expression containing multiple dependent variables and we want to collect coefficient wrt to selected dependent variables. Please see attached Maple file for details.

 >
 >
 (1)
 >
 (2)

Suppose we differential expression like:

 >
 (3)

We can collect coefficients of differential like in following manner:

 The Procedure

Now how can we collect coefficents with respect to so that differential expression (3) appear as

 >

Regards

## Coefficients in differential expression...

Dear Friends

In differential expressions(See Maple file) how to find coefficiets of dependent variable "u(x,t)" and "v(x,t)" and of their differentials ? There is command "dcoeffs(function)", but that work for single dependent variable only but in our case there are two dependent variables in consideration. There are other options like "indets", "specindex" but those do not work.

 >
 >
 (1)
 >
 (2)
 >
 (3)

In above differential expressions how to find coefficiets of dependent variable "u(x,t)" and "v(x,t)" and of their differentials ? There is command "dcoeffs(expr,u(x,t))", but that work for single dependent variable only but in our case there are two dependent variables in consideration. There are other options like "indets", "specindex" but those do not work.

 >

## What is short form of derivative? ...

mydoubt-to-mapleprime.mw
an any one help me understand and correct the following errors in maple document. Thanks.

Ramakrishnan V

rukmini_ramki@hotmail.com

## Coefficients of differential polynomial...

Hello everyone,

I am trying to extract the coefficients from a differential poynomial. In general, this poynomial is in two variables, say u and v along with their differentials, i.e D(u) or D@@2(u) or so on.

Coefficients of this polynomials are rational funcitons.

For instance- consider the following example:

a(x)*v*u+v*D(u)-D(v)*u

then output should be [a(x), 1, -1].

## Solve and plot differential equations...

I want to solve these two differential equations. I have the initial conditions:
x(0)=0
y(0)=0
x'(0)=5.7
t'(0)=8.1
What am I doing wrong?

## New: Factorizing with non-commutative variables

by:

One of the interesting things about the Physics package is that it was designed from scratch to extend the domain of operations of the Maple system from commutative variables to one that includes commutative, anticommutative and nonocommutative variables, as well as abstract vectors and related (nabla) differential operators. In this line we have, among others, the following Physics commands working with this extended domain: `*` , `.` , `^` , diff , Expand , Normal , Simplify , Gtaylor , and Coefficients .

More recently, Pascal Szriftgiser (from Laboratoire PhLAM, Université Lille 1, France), suggested a similar approach to factorize expressions involving noncommutative variables. This is a pretty complicated problem though. Pascal's suggestion, however, spinned around an idea beautiful for its simplicity, similar to what is done in the experimental Physics command, PerformOnAnticommutativeSystem , that is, to transform the problem into one that can be treated with the command that works only with commutative variables and from there extract the result for noncommutative ones.The approach has limitations but it is surprising how far one can go using imaginative algebraic manipulations to extend these commands that otherwise only work with commutative variables.

In brief, we now have a new command, Physics:-Factor, with already powerful performance for factorizing algebraic expressions that involve commutative, noncommutative and anticommutative variables, making Maple's mathematical capabilities more advanced in very interesting directions. This command is in fact useful not just in advanced theoretical physics, but for instance also when working with noncommutative symbols representing abstract matrices (that can have dependency, and so they can be differentiated before saying anything about their components, multiplied, and be present int  expressions that in turn can be expanded, simplified and now also factorized), and also useful with expressions that include differential operators, now that within Physics you can compute with the the covariant and noncovariant derivatives D_  and d_ algebraically. For instance, how about solving differential equations using Physics:-Factor (reducing their order by means of factoring the involved differential operators) ? :)

What follows are some basic algebraic examples illustrating the novelty, and as usual to reproduce the results in this worksheet you need to update your Physics library with the one available in the Maplesoft R&D Physics webpage.

 >
 (1)
 >
 (2)

First example, because of using mathematical notation, noncommutative variables are displayed in different color (olive)

 >
 (3)
 >
 (4)

A more involved example from a physics problem, illustrating that the factorization is also happening within function's arguments, as well as that we can also correctly expand mathematical expressions involving noncommutative variables

 >
 (5)
 >
 (6)

So first expand ...

 >
 (7)

Now retrieve the original expression by recursing over the arguments and so factoring the integrand

 >
 (8)

This following one looks simpler but it is actually more complicated:

 >
 (9)
 >
 (10)

The complication consists of the fact that the standard factor  command, that assumes products are commutative, can never deal with factors like  because if products were commutative these factors are equal to 0. Of course we not just us factor but include a number of algebraic manipulations before using it, so that the approach handles these cases nicely anyway

 >
 (11)

This other one is more complicated:

 >
 (12)

When you expand,

 >
 (13)

you see that there are various terms involving the same noncommutative operands, just multiplied in different order. Generally speaking the limitation (n this moment) of the approach is: "there cannot be more than 2 terms in the expanded form containing the same operands" . For instance in (13) the 1st and 4th terms have the same operands, that are actually also present in the 5th term but there you also have  and for that reason (involving some additional manipulations) it can be handled:

 >
 (14)

Recalling, in all these examples, the task is actually accomplished by the standard factor  command, and the manipulations consist of ingeniously rewriting the given problem as one that involves only commutative variables, and from extract the correct result for non commutative variables.

To conclude, here is an example where the approach implemented does not work (yet) because of the limitation mentioned in the previous paragraph:

 >
 (15)
 >
 (16)

In this expression, the 1st, 2nd, 4th and 5th terms have the same operands  and then there are four terms containing the operands . We do have an idea of how this could be done too ... :) To be there in one of the next Physics updates.

 >

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

## How to solve this BVP problem?...

Hi everyone,

I'm kinda new here, and I really hope you guys can help me through this.

I'm doing a term project and I need to solve this differential equation. I tried to make it work, but since I'm new with maple I couldn't get to any result. So I really appreciate if you can help me solve the problem, since I am eally stuck here and all other results sort of depend on this equation.

here is the equation:

diff(y(x), x, x) = -(8*omega*(1-exp(-8*x))*exp(-8*x)/(2*x^2)-pi^2/(32*x))*y(x)

with following boundary conditions:

y(0)=0,

y(1/2)=0 ,

y'(0)=1.

thanks again.

## PDEs and Boundary Conditions - new developments

by: Maple

The PDE & BC project , a very nice and challenging one, also one where Maple is pioneer in all computer algebra systems, has restarted, including now also the collaboration of Katherina von Bülow.

Recapping, the PDE & BC project started 5 years ago implementing some of the basic methods found in textbooks to match arbitrary functions and constants to given PDE boundary conditions of different kinds. At this point we aim to fill gaps, and the first one we tackled is the case of 1st order PDE that can be solved without boundary conditions in terms of an arbitrary function, and where a single boundary condition (BC) is given for the PDE unknown function, and this BC does not depend on the independent variables of the problem. It looks simple ... It can be rather tricky though. The method we implemented is a simple however ingenious use of differential invariants  to match the boundary condition.

The resulting new code, the portion already tested, is available for download in the Maplesoft R&D webpage for Differential Equations and Mathematical Functions (the development itself is bundled within the library that contains the new developments for the Physics package, in turn within the zip linked in the webpage).

The examples that can now be handled, although restricted in generality to "only one 1st order linear or nonlinear PDE and only one boundary condition for the unknown function itself", illustrate well how powerful it can be to use more advanced methods to tackle these tricky situations where we need to match an arbitrary function to a boundary condition.

To illustrate the idea, consider first a linear example, among the simplest one could imagine:

 >
 (1)
 >
 (2)

Input now a boundary condition (bc) for the unknown such that this bc does not depend on the independent variables ; this bc can however depend on arbitrary symbolic parameters, for instance

 >
 (3)

With the recent development, this kind of problem can now be solved in one go:

 >
 (4)

Nice! And how do you verify this result for correctness? With pdetest , which actually also tests the solution against the boundary conditions:

 >
 (5)

And what has been done to obtain the solution (4)? First the PDE was solved regardless of the boundary condition, so in general, obtaining:

 >
 (6)

In a second step, the arbitrary function  got determined such that the boundary condition  is matched. Concretely, the mapping _F1 is what got determined. You can see this mapping reversing the solving process in two steps. Start taking the difference between the general solution (6) and the solution (4) that matches the boundary condition

 >
 (7)

and isolate here

 >
 (8)

So this is the value  that got determined. To see now the actual solving mapping _F1, that takes for arguments  and  and returns the right-hand side of (8), one can perform a change of variables introducing the two parameters  and  of the _F1 mapping:

 >
 (9)
 >
 (10)
 >
 (11)

So the solving mapping _F1 is

 >
 (12)

Wow! Although this pde & bc problem really look very simple, this solution (12) is highly non-obvious, as is the way to get it just from the boundary condition  and the solution (6) too. Let's first verify that this mapping is correct (even when we know, by construction, that it is correct). For that, apply (12) to the arguments of the arbitrary function and we should obtain (8)

 >
 (13)

Indeed this is equal to (8)

 >
 (14)

Skipping the technical details, the key observation to compute a solving mapping is that, given a 1st order PDE where the unknown depends on  independent variables, if the boundary condition depends on  arbitrary symbolic parameters , one can always seek a "relationship between these parameters and the differential invariants that enter as arguments in the arbitrary function _F1 of the solution", and get the form of the mapping _F1 from this relationship and the bc. The method works in general. Change for instance the bc (3) making its right-hand side be a sum instead of a product

 >
 (15)
 >
 (16)
 >
 (17)

An interesting case happens when the boundary condition depends on less than  parameters, for instance:

 >
 (18)
 >
 (19)

As we see in this result, the additional difficulty represented by having few parameters got tackled by introducing an arbitrary constant _C1 (this is likely to evolve into something more general...)

 >
 (20)

Finally, consider a nonlinear example

 >
 (21)
 >
 (22)

Here we have 2 independent variables, so for illustration purposes use a boundary condition that depends on only one arbitrary parameter

 >
 (23)

All looks OK, but we still have another problem: check the arbitrary function _F1 entering the general solution of pde when tackled without any boundary condition:

 >
 (24)

Remove this RootOf to see the underlying algebraic expression

 >
 (25)

So this is a pde where the general solution is implicit, actually depending on an arbitrary function of the unknown  The code handles this problem in the same way, just that in cases like this there may be more than one solution. For this very particular bc (23) there are actually three solutions:

 >
 (26)

Verify these three solutions against the pde and the boundary condition

 >
 (27)

:)

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

## Skew product does not work properly...

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?

## Derivative from pdsolve/numeric solution?...

Is it possible to somehow extract a derivative from numeric solution of partial differential equation?

I know there is a command that does it for dsolve but i couldn't find the same thing for pdsolve.

The actual problem i have is that i have to take a numeric solution, calculate a derivative from it and later use it somewhere else, but the solution that i have is just a set of numbers an array of some sort and i can't really do that because obviously i will get a zero each time.

Perhaps there is a way to interpolate this numeric solution somehow?

I found that someone asked a similar question earlier but i couldn't find an answer for it.

## How to solve delay differential equations with Map...

How to solve delay differential equations with Maple?

Example:

diff(x(t),t) = 3*x(t)^2 + 0.3*x(x-0.03)

## DELAY DIFFERENTIAL EQUATIONS...

Hello everyoene, please i have a problem solving this delay differential equation:

y'(x)=cos(x)+y(y(x)-2)    0<x<=10

y(x)=1      x<=0

tau=x-y(x)+2

please i need the solution urgently