@Carl Love 

Don't you think an interpretation of the question could be the one given in the attached file ?

The idea is : psi1 and psi2 are functions of x and t. Provided Fubini's conditions hold, differentiating eq1 (eq2) with respect to t and eq3 (eq4) with respect to x and substracting the results pairwise gives a couple of PDE in psi1 and psi2.

psi.mw

@weidade37211 
Your problem remains obscure for me. 
Could you give a concrete example ?


PS : I'm not even sure I could be of any help to you ... me request is simple curiosity

@vv 
Great !

@Kitonum 

Hi, 
A very interesting answer indeed !
I tried to apply it to randomly generated polynomials and I ontained "strange" results (attached file).
For instance m[0,0] sometimes appear, or monomials in only one indet, such as "y", is not represented as m[0,1].

I guess your answer fits well MarcoFerro's question , but could you generalize it to handle the test cases in the attached file ?
 

Great thanks in advance

(PS : I use Maple 2015.2)

CanYouExplainThisToMe.mw

@shaif 

PS : the problems MAPLE faces in manipulating the GBD come from the exponents (alpha-1), alpha and (b-1).
For rational exponents some tests suggest MAPLE could be capable to return the closed form (while extremely complicated, see next example) of PDF(prod, t)

pp := int( subs({alpha=1, b=1}, subs(t=u/t, p1) * p4 / t), t=0..+infinity) assuming u > 0:
pp := int( subs({alpha=3, b=1}, subs(t=u/t, p1) * p4 / t), t=0..+infinity) assuming u > 0:
pp := int( subs({alpha=3/2, b=2}, subs(t=u/t, p1) * p4 / t), t=0..+infinity) assuming u > 0:
pp := int( subs({alpha=3/8, b=4/7}, subs(t=u/t, p1) * p4 / t), t=0..+infinity) assuming u > 0:

@shaif 

First point, my previous reply contained some mistakes, with this major one
beta(1+((x-gamma)/beta)^alpha)^(k+1) 
instead of 
beta*(1+((x-gamma)/beta)^alpha)^(k+1)

apologies ...


Second point: your Burr's distribution is a generalizae one (GBD for short in the sequel).
Let me denote SBD the standard Burr's distribution 

StandardBurrPdf          := piecewise(x<0, 0, alpha*k*(x^(alpha-1) / (1+x^alpha)^(k+1)) );
StandardBurrDistribution := Distribution(PDF = unapply(StandardBurrPdf, x), Conditions=[alpha > 0, k > 0]);

Please note the piecewised definition of StandardBurrPdf.
Some trails show that MAPLE is able to compute some statistics of a SBD random variable (Mean, CDF, ...)

What about your GBD ?
Here (x-gamma)/beta has to be positive which, because of the positiveness of beta and gamma, reduces to x > gamma.
I think the correct definition of the GBD then is :
GeneralizedBurrPdf := piecewise(x<gamma, 0, alpha*k*((x-gamma)/beta)^(alpha-1) / (beta*(1+((x-gamma)/beta)^alpha)^(k+1)) );
GeneralizedBurrDistribution := Distribution(PDF = unapply(burrpdf, x), Conditions=[alpha > 0, beta > 0, gamma >= 0, k > 0]);

The problem is that MAPLE can do very little with this distribution (I hope someone here will read this and contradict me) unless we help it (see a simple example in the attached file).


Finally, with this "negative" observation, it is very unlikely for MAPLE to be able to return a closed form of PDF(prod, t)

I remain at your disposal
 

burr2.mw

@Preben Alsholm 
Thanks fot having corrected me

@Kitonum 
I use to use  sort(L, key=(x=x[1]))  instead.
Is your syntax docummented elsewhere ?

PS : see Mariusz Iwaniuk's answer to your previous post.
This will explain why CDF(burr_distribution, x), or Mean(burr_distribution), or any other statistics has no closed form for general values of the parameters

Maybe this ?

burr.mw

@Les 

The fit you propose seems odd ...

Maybe a typo error in your f:= t->439e^-2.9t function ?

 

Download FIT_or_fit.mw

Here is your complete worksheet in WORKSHEET-MODE style

EulerLagrange.mw


Don't be afraid if the copy-paste operation transforms the first line into the second one ... the results are the same

If you want to build the solution "sol2" from t=1000 to t=3000, the simplest way is to set the IC equal to sol(1000) and solve in the range 1000..3000.
If you do this (brown code in the attached file), you will see that sol2 is equal to sol in the range 1000..3000

aa_problem_MWE2.mw

I did not analyze your code for solution "sol2".

But, if you try to solve the same problem from t=1000 to t=3000, the only thing tou have to do is to set the new IC at t=1000 equal to sol(1000) and solve in the range 1000..3000
If you do this, sol2 returns the solution sol gave in the range 1000..3000 (see the code in brown in the attached file).

By the way, I added legend=... in the graphs for a smarter plot and change method=rkf45 (default choice) ny method=rosenbrock (to hande potential stiff problems).

So I guess your code for "sol2" could contain some mistake(s) ?


aa_problem_MWE2.mw
 

I did not analyze your code for solution "sol2".

But, if you try to solve the same problem from t=1000 to t=3000, the only thing tou have to do is to set the new IC at t=1000 equal to sol(1000) and solve in the range 1000..3000
If you do this, sol2 returns the solution sol gave in the range 1000..3000 (see the code in brown in the attached file).

By the way, I added legend=... in the graphs for a smarter plot and change method=rkf45 (default choice) ny method=rosenbrock (to hande potential stiff problems).

So I guess your code for "sol2" could contain some mistake(s) ?


aa_problem_MWE2.mw