how to convert (a+b)^3 to a list [a+b, a+b, a+b]?

for example

(a+b)^2 to [a+b, a+b]

I would like work in Riemann Normal Coordinates, and derive expansion in number of derivatives of the metric.

For example, starting with the expansion of the metric in Riemann Normal Coordinates (assuming this needs not be derived)

I would like to express the lapliacian, or square-root of the determinant ... etc., in terms of this metric and derive an expansion for them.

However, I cannot even define the metric as such in Maple-Physics, because the coefficients of x^k depend themselfs upon the metric to be defined.

Is there a way to do such calculations in Maple ?

how to count how many terms or items are equal when compare two lists of polynomial terms when length of two list may not be equal

cannot find the error in loop

> restart; u[0] := (4/3)*c^2*cos((1/4)*x)^2; alpha := 2;
 
> iteration := 3;
> for k from 0 while k <= Iteration do u[s] := eval(u[k], t = xi); u[k+1] := simplify(u[k]-(int(diff(u[s], [`$`(xi, alpha)])+diff(u[s]*u[s], x)+diff(u[s]*u[s], x, x, x), xi = 0 .. t))) end do;

Here's a short tensor manipulation which goes totally bananas. Basically I have a metric, and I define a vector k, and right at the end I calcualte the covariant derivative of it. In the metric and elsewhere I have a constant epsilon, a function of time a(t), and 2 functions of all coords, Phi and psi. Then at the end it gives the covariant derivative of k, but now with epsilon and a as functions of all coordinates. 

Any idea what's going on??

code:

restart;
with(Physics);with(PDEtools):
Setup(coordinatesystems=spherical);

ds2:=expand(a(t)^2*(-(1+2*epsilon*Phi(r,theta,phi,t))*dt^2+(1-2*epsilon*psi(r,theta,phi,t))*(dr^2+r^2*dtheta^2+r^2*sin(theta)^2*dphi^2)));
declare(%,(H)(t),(k0,k1,k2,k3)(X));
Setup(metric = ds2);

Define(k[~mu]=[-1/a(t)^2+epsilon*k1(X),epsilon*k2(X),epsilon*k3(X),1/a(t)^2+epsilon*k0(X)]);

Define(NullGDE[nu]=(k[~mu](X)*D_[mu](k[nu](X))),OpticalDeformationMatrix[mu,nu]=(D_[mu](k[nu](X))));

(OpticalDeformationMatrix[1,1,nonzero]);

the final output is Matrix(1, 1, [[(-(2*((diff(epsilon(X), r))*(psi(X))(X))+2*epsilon(X)*(diff((psi(X))(X), r)))*(epsilon(X)*(k1(X))(X)*((a(t))(X)^2)-1)*a(t)-(2*epsilon(X)*(psi(X))(X)-1)*((diff(epsilon(X), r))*(k1(X))(X)*((a(t))(X)^2)+epsilon(X)*((diff((k1(X))(X), r))*((a(t))(X)^2))+2*epsilon(X)*(k1(X))(X)*(a(t))(X)*(diff((a(t))(X), r)))*a(t)+epsilon*psi[r]*(epsilon*k1(X)*(a(t)^2)-1)*a(t)-epsilon^2*psi[theta]*(a(t)^3)*k2(X)-epsilon^2*psi[phi]*(a(t)^3)*k3(X)-(2*(diff(a(t), t))*psi(X)*epsilon+psi[t]*a(t)*epsilon-(diff(a(t), t)))*(epsilon*k0(X)*(a(t)^2)+1))*(1/a(t))]])

where epsilon is now a function!

This may be a silly question, but does there exist some simple way of (Taylor) expanding an expression of 'small' functions in terms of these functions.

A simple example: Assume that diff(f(x),x) and g(x) are two functions both with range, say, in [-a,+a], where a << 1, and consider the following expression:

sqrt(1 + diff(f(x),x)) * (2 + g(x));

Its expansion to first order in terms of diff(f(x),x) and g(x) should be 2 + diff(f(x),x) + g(x). My problem is that mtaylor does not accept functions as variables to expand on, and I would prefer not to have to substitute back and forth with some 'placeholders'.

The help text for dsolve,numeric,events describes many kinds of triggers and actions and other event specifications, but shows only a very limited selection of examples of these possibilities.

Please tell me of any sources containing extensive examples of the wide variety of available specifications and the problems they are useful in solving.

> restart:with(plots):blt:=7:

> lambda:=2:m:=3:s:=1:

> Eq1:=(diff(f(eta),eta,eta,eta))+(f(eta)*(diff(f(eta),eta,eta)))-((diff(f(eta),eta))^2)+lambda*(((f(eta)*(diff(f(eta),eta,eta,eta))))-2*(diff(f(eta),eta))*(diff(f(eta),eta,eta,eta))^2)-(M/(1+m^2))*((diff(f(eta),eta)+ms))=0;

>

> Eq2:=(diff(h(eta),eta,eta))+(f(eta)*(diff(h(eta),eta)))-((diff(f(eta),eta))*(h(eta)))+lambda*(((f(eta)*(diff(h(eta),eta,eta,eta))))+(h(eta)*(diff(f(eta),eta,eta,eta)))+(diff(f(eta),eta,eta))*(diff(h(eta),eta))-2*(diff(f(eta),eta))*(diff(h(eta),eta,eta))+(M/(1+m^2)))*(m*(diff(f(eta),eta)-h))=0;

>

> Eq3:=((f(eta))*(diff(theta(eta),eta)))+Pr*((diff(theta(eta),eta,eta)))=0;

>

> bcs1 := f(0) = 0, (D(f))(0) = 1, (h)(0) = 0, (theta)(0) = 1, (D(f))(blt) = 0,  h(blt) = 0, theta(blt) = 0, D(D(f))(blt) = 0, (D(h))(blt)=0;

>

> L := [1,2,3];

> for k to 3 do R := dsolve(eval({Eq1, Eq2, Eq3, bcs1}, M = L[k]), [f(eta),h(eta),theta(eta)], numeric, maxmesh=10000, output = listprocedure);;Y || k := rhs(R[3]); YA || k := rhs(R[6]);YB || k := rhs(R[5]);YC || k := rhs(R[4]);YD || k := rhs(R[7]);end do:

>

>

>

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

for example is there an existing package for reducing groups of matrices like the one below to only its unique elements, or do i basically need to use  linear algebra matrix operations ie finding the basis of the set via echeleon reduction blah blah

{Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = -673/2880}), Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = -5/96}), Matrix(2, 2, {(1, 1) = 0, (1, 2) = -(2521/17920)*Zeta(5), (2, 1) = 0, (2, 2) = -(2521/17920)*Zeta(3)-2087/1920}), Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = 0}), Matrix(2, 2, {(1, 1) = 0, (1, 2) = -(7/320)*Zeta(5), (2, 1) = 0, (2, 2) = -(7/320)*Zeta(3)-499/840})}

i want a scheme of fractional differential equation so that i solve my questions and make a code of it.

please provide me the scheme

I was learning about Maple, but I do not understand this code below. This does not look correct.

http://www.maplesoft.com/support/help/Maple/view.aspx?path=Task/ConvertExpressionToFunction

When I type it on my worksheet, I get these errors

Sometimes I wonder if any one at Maplesoft actually looks at their own help web pages and try to use them, or is just the poor users who do that.

I noticed Maple gives different order of term (but still correct) when variable y or z is used, vs. other letters in this computation.

In this example below, integral results are given, using one letter of the alphabet at a time in the integrand. 

Only when the variable is y or z, Maple shows expression in different form from all the others (negative sign pulled out). The result are correct ofcourse, but was wondering why this hapens only for these two variables? I would have expected that variable name (letter used) should not make any difference to the final form.

restart;
interface(rtablesize=infinity):
chars:=[seq(parse(StringTools:-Char(i)),i=97..122)]: #generate a..z variables, may be there is better way?
Matrix([seq([chars[i0],simplify(int((c0^2 - c1^2)/(chars[i0] - x0+I*c3), x0))],i0=1..nops(chars))]);

Gives (this below. Notice the very last 2 entries below)

Is this something common in Maple? Does it depend on the computation being performed? Can a user do something to tell Maple not to change order/form of expression depending on what variable letter is used?

Download integral.mw

Download integral.mwintegral.mw

hi

please help me for dsolve differential equations...after much time dont answer!!!

thanks..

dsolve.mw

Download dsolve.mw

Can you tell me the signification of the number 5 in output below.

Input: maximize(x^2+3,x=-1..2,location)

Output: 7,{[{x=2},5]}

Thanks