Items tagged with expression expression Tagged Items Feed


Suppose I have a triplet a[5], a[6], a[7]and an expression given as below:

V[5]*a[5]+V[6]*a[6]+V[7]*a[7]+V[2]+V[9]

How I can list all possible forms of above expression for triplet "(a[5], a[6], a[7])->(a[5], a[6], a[7]), (a[5], a[6], 0), (a[5], 0, a[7]), (0, a[6], a[7]), (a[5], 0, 0), (0, a[6], 0 ), (0, 0, a[7]), (0, a[6], 0), (0, 0, 0) "etc


Download triplet.mw

Regards

Hi,

can you please help me with the usage of applyrule? I have the following problem that I cannot isolate. I have a rule that I want to apply, but instead of applying the rule to the expression, the rule seems to be applied to itself. In an isolated worksheet everything is fine:

restart;

rule:=abs(''a''::algebraic)^2=''a''^2;

abs('a'::algebraic)^2 = 'a'^2

(1)

 

myexp:=abs(548.477146186283171377723+radius_motor*q_mot_vec_2(t)-l_wire_0[2])^2

abs(548.477146186283171377723+radius_motor*q_mot_vec_2(t)-l_wire_0[2])^2

(2)

applyrule(rule,myexp);

(548.477146186283171377723+radius_motor*q_mot_vec_2(t)-l_wire_0[2])^2

(3)

rule

abs(a::algebraic)^2 = a^2

(4)

 


Download applyrule.mw

What happens when I try to use this rule in my script is this:

restart;

#read "some_long_script.mpl":

rule:=abs(''a''::algebraic)^2=''a''^2;

abs('a'::algebraic)^2 = 'a'^2

(1)

applyrule(rule,abs(x)^2)

abs(x)^2

(2)

rule

a::algebraic^2 = a^2

(3)

 

>
 

I want to clarify, that the script that is read before applying those rules does not apply a similar rule. Also the variables rule and x were free. Can you please help me to locate the problem?

Thanks!

Honigmelone

Hi all,

I am trying to do some calculation using Maple. I have different expressions containing some power of y. I would like that Maple automatically substitutes y^4 = f(y), in each one of these. So for example, if my expression contains y^5, I want it substituted by y * f(y), and so on. How can I do?

Dear All

I want to use subs command to replace products in expression with entries from matrix. When I use "subs" command it returns same result as original without any substitution. But when I use "algsubs", it works fine, but is complicated to apply if there are large number of substitutions.

Please see followings:

 

M := Matrix([[U[1], U[2], -epsilon*U[1]+U[3], U[4], U[5], U[6]], [U[1], U[2], -epsilon*U[2]+U[3], epsilon*U[1]+U[4], epsilon*U[2]+U[5], U[6]], [e^epsilon*U[1], e^epsilon*U[2], U[3], U[4], U[5], U[6]], [U[1], -epsilon*U[1]+U[2], U[3], U[4], -epsilon*U[4]+U[5], U[6]-epsilon*(U[3]+2*U[5])+epsilon^2*U[4], U[1], e^(-epsilon)*U[2], U[3], e^epsilon*U[4], U[5], e^(-epsilon)*U[6]], [U[1], e^(-epsilon)*U[2], U[3], e^epsilon*U[4], U[5], e^(-epsilon)*U[6]], [epsilon*U[2]+U[1], U[2], U[3], U[4]+epsilon*(U[3]+2*U[5])+epsilon^2*U[6], epsilon*U[6]+U[5], U[6]]])

M := Vector(4, {(1) = ` 6 x 12 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

(1)

M[5, 1]; 1; M[5, 3]

U[1]

 

U[3]

(2)

expand((w*T[5]+z*T[4]+T[6])*(U[1]*a[1]+U[2]*a[2]+U[3]*a[3]+U[4]*a[4]+U[5]*a[5]+U[6]*a[6]))

w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]

(3)

subs(T[5]*U[1] = M[5, 1], T[5]*U[2] = M[5, 2], T[5]*U[3] = M[5, 3], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6])

w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]

(4)

Why this command subs in not replacing  products  T[5]*U[1], T[5]*U[2], T[5]*U[3] from matrix (1)NULL?????

algsubs(T[5]*U[1] = M[5, 1], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6])

z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]+(w*U[2]*a[2]+w*U[3]*a[3]+w*U[4]*a[4]+w*U[5]*a[5]+w*U[6]*a[6])*T[5]+(z*T[4]*a[1]+w*a[1]+T[6]*a[1])*U[1]

(5)

It worked !!!! But If I try "algsubs" like :

algsubs(T[5]*U[1] = M[5, 1], T[5]*U[2] = M[5, 2], T[5]*U[3] = M[5, 3], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6])

Error, invalid input: algsubs expects its 3rd argument, x, to be of type {name, function, list({name, function}), set({name, function})}, but received T[5]*U[3] = U[3]

 

See it is giving error !!! but still I can use "algsubs" as follow:

 

expand(algsubs(T[6]*U[6] = M[6, 6], expand(algsubs(T[6]*U[5] = M[6, 5], expand(algsubs(T[6]*U[4] = M[6, 4], expand(algsubs(T[6]*U[3] = M[6, 3], expand(algsubs(T[6]*U[2] = M[6, 2], expand(algsubs(T[6]*U[1] = M[6, 1], expand(algsubs(T[5]*U[6] = M[5, 6], expand(algsubs(T[5]*U[5] = M[5, 5], expand(algsubs(T[5]*U[4] = M[5, 4], expand(algsubs(T[5]*U[3] = M[5, 3], expand(algsubs(T[5]*U[2] = M[5, 2], expand(algsubs(T[5]*U[1] = M[5, 1], expand(algsubs(T[4]*U[6] = M[4, 6], expand(algsubs(T[4]*U[5] = M[4, 5], expand(algsubs(T[4]*U[4] = M[4, 4], expand(algsubs(T[4]*U[3] = M[4, 3], expand(algsubs(T[4]*U[2] = M[4, 2], expand(algsubs(T[4]*U[1] = M[4, 1], expand(algsubs(T[3]*U[6] = M[3, 6], expand(algsubs(T[3]*U[5] = M[3, 5], expand(algsubs(T[3]*U[4] = M[3, 4], expand(algsubs(T[3]*U[3] = M[3, 3], expand(algsubs(T[3]*U[2] = M[3, 2], expand(algsubs(T[3]*U[1] = M[3, 1], expand(algsubs(T[2]*U[6] = M[2, 6], expand(algsubs(T[2]*U[5] = M[2, 5], expand(algsubs(T[2]*U[4] = M[2, 4], expand(algsubs(T[2]*U[3] = M[2, 3], expand(algsubs(T[2]*U[3] = M[2, 3], expand(algsubs(T[2]*U[2] = M[2, 2], expand(algsubs(T[2]*U[1] = M[2, 1], expand(algsubs(T[1]*U[6] = M[1, 6], expand(algsubs(T[1]*U[5] = M[1, 5], expand(algsubs(T[1]*U[4] = M[1, 4], expand(algsubs(T[1]*U[3] = M[1, 3], expand(algsubs(T[1]*U[2] = M[1, 2], expand(algsubs(T[1]*U[1] = M[1, 1], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))

a[5]*epsilon*U[6]+a[1]*epsilon*U[2]+U[6]*a[6]+U[2]*a[2]+U[3]*a[3]+U[4]*a[4]+U[5]*a[5]+z*a[1]*U[1]+epsilon^2*U[6]*a[4]+epsilon*U[3]*a[4]+2*epsilon*U[5]*a[4]+U[6]*w*a[6]/e^epsilon+U[4]*e^epsilon*w*a[4]+w*a[1]*U[1]+U[2]*w*a[2]/e^epsilon+epsilon^2*z*U[4]*a[6]-epsilon*z*U[3]*a[6]-2*epsilon*z*U[5]*a[6]+w*U[3]*a[3]+w*U[5]*a[5]-z*a[5]*epsilon*U[4]-z*a[2]*epsilon*U[1]+z*U[2]*a[2]+z*U[3]*a[3]+z*U[4]*a[4]+z*U[5]*a[5]+z*U[6]*a[6]+a[1]*U[1]

(6)

But this procedure is complex, I want to use "subs" only. I know it complexity of terms in (3) might me creating problem.

collect(a[5]*epsilon*U[6]+a[1]*epsilon*U[2]+U[6]*a[6]+U[2]*a[2]+U[3]*a[3]+U[4]*a[4]+U[5]*a[5]+z*a[1]*U[1]+epsilon^2*U[6]*a[4]+epsilon*U[3]*a[4]+2*epsilon*U[5]*a[4]+U[6]*w*a[6]/e^epsilon+U[4]*e^epsilon*w*a[4]+w*a[1]*U[1]+U[2]*w*a[2]/e^epsilon+epsilon^2*z*U[4]*a[6]-epsilon*z*U[3]*a[6]-2*epsilon*z*U[5]*a[6]+w*U[3]*a[3]+w*U[5]*a[5]-z*a[5]*epsilon*U[4]-z*a[2]*epsilon*U[1]+z*U[2]*a[2]+z*U[3]*a[3]+z*U[4]*a[4]+z*U[5]*a[5]+z*U[6]*a[6]+a[1]*U[1], [U[1], U[2], U[3], U[4], U[5], U[6]])

(-epsilon*z*a[2]+w*a[1]+z*a[1]+a[1])*U[1]+(a[1]*epsilon+a[2]+w*a[2]/e^epsilon+z*a[2])*U[2]+(-epsilon*z*a[6]+epsilon*a[4]+w*a[3]+z*a[3]+a[3])*U[3]+(epsilon^2*z*a[6]+e^epsilon*w*a[4]-z*a[5]*epsilon+z*a[4]+a[4])*U[4]+(-2*epsilon*z*a[6]+2*epsilon*a[4]+w*a[5]+z*a[5]+a[5])*U[5]+(a[5]*epsilon+a[6]+epsilon^2*a[4]+w*a[6]/e^epsilon+z*a[6])*U[6]

(7)

``

 

Download Subs_not_giving_result.mwSubs_not_giving_result.mw

Regards

 

2*t*exp(t) + 2*t

would like to apply op recurively 

and output list

[2,t,exp(t),2,t]

Hi all!

as shown below, how can get a result without 'R':

p_com(z,t):=Re(exp(I*omega*t-I*k*(lambda[r]+I*lambda[i])*z)) assuming omega::real,t::real,k::real,lambda[r]::real,lambda[i]::real,z::real

Thanks very much!

 

Hi Guys, 

I'm trying to find stationary points of a numeric function. Any help would be appreciated. 

Assume I have a numeric function g(x). I'm attempting the following:

deriv:=(x)->fdiff(g(t),t=x); (I can't use the D operator as it doesn't like g(x))

Now deriv(x) <- Can be computed and computes the value quickly and easily. However plot and other functions require algebraic functions. I can plot this (with difficulty) via the following

plot(deriv,-10..10) <- This avoids converting to algebraic function and runs as a procedure. 

I can't use convert(deriv,algebraic) as it fails. 

I want to use the Student[NumericalAnalysis]Roots command, but it requires an algebraic function and I can't use the procedure trick which I did in the plot command. 

Does anyone know a better way of doing this? or is there a way I can write the numerical differivative as an algebraic function. (I've tried fsolve <- But this guy doesn't give me the correct answers generally). 

Thanks guys. 

 

 

 

Dear All

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

See for detail:


with(PDEtools):

u := a[0]+a[1]*(diff(phi(xi), xi))/phi(xi)

a[0]+a[1]*(diff(phi(xi), xi))/phi(xi)

(1)

-k*(diff(u, `$`(xi, 2)))+alpha*(diff(u, `$`(xi, 3)))

-k*(a[1]*(diff(diff(diff(phi(xi), xi), xi), xi))/phi(xi)-3*a[1]*(diff(diff(phi(xi), xi), xi))*(diff(phi(xi), xi))/phi(xi)^2+2*a[1]*(diff(phi(xi), xi))^3/phi(xi)^3)+alpha*(a[1]*(diff(diff(diff(diff(phi(xi), xi), xi), xi), xi))/phi(xi)-4*a[1]*(diff(diff(diff(phi(xi), xi), xi), xi))*(diff(phi(xi), xi))/phi(xi)^2+12*a[1]*(diff(diff(phi(xi), xi), xi))*(diff(phi(xi), xi))^2/phi(xi)^3-3*a[1]*(diff(diff(phi(xi), xi), xi))^2/phi(xi)^2-6*a[1]*(diff(phi(xi), xi))^4/phi(xi)^4)

(2)

expand(dsubs(diff(phi(xi), `$`(xi, 2)) = -lambda*(diff(phi(xi), xi))-mu*phi(xi), -k*(a[1]*(diff(diff(diff(phi(xi), xi), xi), xi))/phi(xi)-3*a[1]*(diff(diff(phi(xi), xi), xi))*(diff(phi(xi), xi))/phi(xi)^2+2*a[1]*(diff(phi(xi), xi))^3/phi(xi)^3)+alpha*(a[1]*(diff(diff(diff(diff(phi(xi), xi), xi), xi), xi))/phi(xi)-4*a[1]*(diff(diff(diff(phi(xi), xi), xi), xi))*(diff(phi(xi), xi))/phi(xi)^2+12*a[1]*(diff(diff(phi(xi), xi), xi))*(diff(phi(xi), xi))^2/phi(xi)^3-3*a[1]*(diff(diff(phi(xi), xi), xi))^2/phi(xi)^2-6*a[1]*(diff(phi(xi), xi))^4/phi(xi)^4)))

-a[1]*(diff(phi(xi), xi))*alpha*lambda^3/phi(xi)-a[1]*alpha*lambda^2*mu-7*a[1]*(diff(phi(xi), xi))^2*alpha*lambda^2/phi(xi)^2-8*a[1]*(diff(phi(xi), xi))*alpha*lambda*mu/phi(xi)-a[1]*(diff(phi(xi), xi))*k*lambda^2/phi(xi)-2*a[1]*alpha*mu^2-a[1]*k*lambda*mu-12*a[1]*(diff(phi(xi), xi))^3*alpha*lambda/phi(xi)^3-8*a[1]*(diff(phi(xi), xi))^2*alpha*mu/phi(xi)^2-3*a[1]*(diff(phi(xi), xi))^2*k*lambda/phi(xi)^2-2*a[1]*(diff(phi(xi), xi))*k*mu/phi(xi)-6*a[1]*(diff(phi(xi), xi))^4*alpha/phi(xi)^4-2*a[1]*(diff(phi(xi), xi))^3*k/phi(xi)^3

(3)

How ca extract coefficients of fraction (diff(phi(xi), xi))/phi(xi) in (3) ????


Download Coefficients_of_Fractions.mw

Regards

I am wondering why Maple does this.

> f:=x^2
> f(3)
output: x(3)^2

I understand the difference between an expression and a function. If f is an expression, shouldn't it ouput  x^2(3). Why is the output x(3)^2?

When I enter f*3 or f*(3) then I get the correct expression.

Here is a screenshot. http://prntscr.com/a7u9hm

Here is image inserted with a slightly different function

 

Also while I am here, what exactly does g(x):= x^2 do? when i enter g(3) I get g(3) back.

g(x) is neither a function nor an expression.

screenshot http://prntscr.com/a7ua75



 I am haunted by an equation which I could simplify it with hand and it equals to −i (complexe number). But, I don't know how to use maplesoft to simplify it. The equation takes the form of:

 

-Omega*a*sqrt(2)*sqrt(-Omega^2*a^2-2*k*m+sqrt(Omega^2*a^2*(Omega^2*a^2+4*k*m)))/(-Omega^2*a^2+sqrt(Omega^2*a^2*(Omega^2*a^2+4*k*m)))



I could simplify this equation with hand calculation and it equals to -i (the complexe number). I'am sorry for not clarifying that the a, k, Ω and m are positive variables.

Thank you in advance for taking a look.

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.

 


with(PDEtools):

DepVars := [u(x, t), v(x, t), a[1](t), a[2](t), a[3](t), b[1](t), b[2](t), b[3](t), r(x, t), s[1](x, t), p[1](x, t), s[2](x, t), p[2](x, t)]

[u(x, t), v(x, t), a[1](t), a[2](t), a[3](t), b[1](t), b[2](t), b[3](t), r(x, t), s[1](x, t), p[1](x, t), s[2](x, t), p[2](x, t)]

(1)

alias(u = u(x, t), v = v(x, t), a[1] = a[1](t), a[2] = a[2](t), a[3] = a[3](t), b[1] = b[1](t), b[2] = b[2](t), b[3] = b[3](t), r = r(x, t), s[1] = s[1](x, t), p[1] = p[1](x, t), s[2] = s[2](x, t), p[2] = p[2](x, t))

u, v, a[1], a[2], a[3], b[1], b[2], b[3], r, s[1], p[1], s[2], p[2]

(2)

Suppose we differential expression like:

a[1]*(diff(u, x))*s[1]*u-2*a[1]*u*(diff(r, x))*(diff(u, x))+2*a[2]*(diff(v, x))*s[2]*v-2*a[2]*v*(diff(r, x))*(diff(v, x))-(diff(a[3], t))*r*(diff(u, x))/a[3]+diff(p[1], t)+a[3]*(diff(p[1], x, x, x))+(diff(r, t))*(diff(u, x))+(diff(s[1], t))*u-(diff(a[3], t))*s[1]*u/a[3]-s[1]*a[2]*v*(diff(v, x))-(diff(a[3], t))*a[1]*u*(diff(u, x))/a[3]-(diff(a[3], t))*a[2]*v*(diff(v, x))/a[3]-3*(diff(r, x))*p[1]+(diff(a[1], t))*u*(diff(u, x))+(diff(a[2], t))*v*(diff(v, x))+a[2]*(diff(v, x))*p[2]+a[2]*v^2*(diff(s[2], x))+a[2]*v*(diff(p[2], x))+a[1]*u*(diff(p[1], x))+a[1]*(diff(u, x))*p[1]+a[1]*u^2*(diff(s[1], x))+3*a[3]*(diff(s[1], x))*(diff(u, x, x))+3*a[3]*(diff(s[1], x, x))*(diff(u, x))+a[3]*(diff(r, x, x, x))*(diff(u, x))-(diff(a[3], t))*p[1]/a[3]-3*r*(diff(r, x))*(diff(u, x))-3*(diff(r, x))*s[1]*u+a[3]*(diff(s[1], x, x, x))*u+3*a[3]*(diff(r, x, x))*(diff(u, x, x)) = 0

3*a[3]*(diff(diff(r, x), x))*(diff(diff(u, x), x))+3*a[3]*(diff(s[1], x))*(diff(diff(u, x), x))+3*a[3]*(diff(diff(s[1], x), x))*(diff(u, x))+a[3]*(diff(diff(diff(r, x), x), x))*(diff(u, x))+a[3]*(diff(diff(diff(s[1], x), x), x))*u+diff(p[1], t)+(diff(r, t))*(diff(u, x))+(diff(s[1], t))*u-3*(diff(r, x))*p[1]+a[3]*(diff(diff(diff(p[1], x), x), x))-(diff(a[3], t))*a[1]*u*(diff(u, x))/a[3]-(diff(a[3], t))*a[2]*v*(diff(v, x))/a[3]+a[1]*(diff(u, x))*s[1]*u-2*a[1]*u*(diff(r, x))*(diff(u, x))+2*a[2]*(diff(v, x))*s[2]*v-2*a[2]*v*(diff(r, x))*(diff(v, x))-(diff(a[3], t))*r*(diff(u, x))/a[3]-(diff(a[3], t))*s[1]*u/a[3]-s[1]*a[2]*v*(diff(v, x))+(diff(a[1], t))*u*(diff(u, x))+a[1]*u*(diff(p[1], x))+a[2]*v*(diff(p[2], x))+a[2]*v^2*(diff(s[2], x))+a[2]*(diff(v, x))*p[2]+a[1]*(diff(u, x))*p[1]+a[1]*u^2*(diff(s[1], x))-(diff(a[3], t))*p[1]/a[3]-3*r*(diff(r, x))*(diff(u, x))-3*(diff(r, x))*s[1]*u+(diff(a[2], t))*v*(diff(v, x)) = 0

(3)

We can collect coefficients of differential like u[x], u[x, x], v[x], u, vin following manner:

The Procedure

   

 

 

Now how can we collect coefficents with respect to u[x], u[x, x], v[x], u, vso that differential expression (3) appear as
"(......)*u+(.......)*v+(......)*u[x]+(........)*uu[x]+(.........)vv[x]+(........)u[xx]  =0....................."????????""

``


Download Collecting_Coefficients_in_Differential_Expression.mw

Regards

guys,i computed a tensorial expression by maple but i think i made mistake.

 

 

tensorial.mw

sqrt(5) gives sqrt(5)

sqrt(1+sqrt(5)) gives "You have entered an invalid Maple expression"

sqrt(u) gives sqrt(u)

sqrt(1+u); gives "You have entered an invalid Maple expression"

 

when using the Maple Math icon. How can I get the correct input for the two expressions?

I have the following function

where A,B,Ψ, K1,K2,K3,α,β are all constants.

How to find the value of m for which the above expression is 0 or approximate to 0 for different values fo the constants.

e.g., Fixing all the parameters except A, I want to find the values of m for different values of A. How to do that in maple?

 

Hello! Hope every is fine. I want to expand the following expression

exp(2*c*t+2*d*n-d)*alpha*c*a[0]*b[1]^2-exp(2*c*t+2*d*n-d)*alpha*c*a[1]*b[0]*b[1]-exp(2*c*t+2*d*n-d)*alpha*a[0]*a[1]*b[1]+exp(2*c*t+2*d*n-d)*alpha*a[1]^2*b[0]+exp(2*c*t+2*d*n)*alpha*a[0]*a[1]*b[1]-exp(2*c*t+2*d*n)*alpha*a[1]^2*b[0]+exp(c*t+d*n)*alpha*c*a[0]*b[0]*b[1]-exp(c*t+d*n)*alpha*c*a[1]*b[0]^2-exp(2*c*t+2*d*n-d)*a[0]*b[1]^2+exp(2*c*t+2*d*n-d)*a[1]*b[0]*b[1]-exp(c*t+d*n-d)*alpha*a[0]^2*b[1]+exp(c*t+d*n-d)*alpha*a[0]*a[1]*b[0]+exp(2*c*t+2*d*n)*a[0]*b[1]^2-exp(2*c*t+2*d*n)*a[1]*b[0]*b[1]+exp(c*t+d*n)*alpha*a[0]^2*b[1]-exp(c*t+d*n)*alpha*a[0]*a[1]*b[0]-exp(c*t+d*n-d)*a[0]*b[0]*b[1]+exp(c*t+d*n-d)*a[1]*b[0]^2+exp(c*t+d*n)*a[0]*b[0]*b[1]-exp(c*t+d*n)*a[1]*b[0]^2

 

like this 

exp(2*c*t+2*d*n)*exp(-d)*alpha*c*a[0]*b[1]^2-exp(2*c*t+2*d*n)*exp(-d)*alpha*c*a[1]*b[0]*b[1]-exp(2*c*t+2*d*n)*exp(-d)*alpha*a[0]*a[1]*b[1]+exp(2*c*t+2*d*n)*exp(-d)*alpha*a[1]^2*b[0]+exp(2*c*t+2*d*n)*alpha*a[0]*a[1]*b[1]-exp(2*c*t+2*d*n)*alpha*a[1]^2*b[0]+exp(c*t+d*n)*alpha*c*a[0]*b[0]*b[1]-exp(c*t+d*n)*alpha*c*a[1]*b[0]^2-exp(2*c*t+2*d*n)*exp(-d)*a[0]*b[1]^2+exp(2*c*t+2*d*n)*exp(-d)*a[1]*b[0]*b[1]-exp(c*t+d*n)*exp(-d)*alpha*a[0]^2*b[1]+exp(c*t+d*n)*exp(-d)*alpha*a[0]*a[1]*b[0]+exp(2*c*t+2*d*n)*a[0]*b[1]^2-exp(2*c*t+2*d*n)*a[1]*b[0]*b[1]+exp(c*t+d*n)*alpha*a[0]^2*b[1]-exp(c*t+d*n)*alpha*a[0]*a[1]*b[0]-exp(c*t+d*n)*exp(-d)*a[0]*b[0]*b[1]+exp(c*t+d*n)*exp(-d)*a[1]*b[0]^2+exp(c*t+d*n)*a[0]*b[0]*b[1]-exp(c*t+d*n)*a[1]*b[0]^2

i.e., expand exp(2*c*t+2*d*n-d) into exp(2*c*t+2*d*n)*exp(-d) 

waiting your kind response 

1 2 3 4 5 6 7 Last Page 1 of 10