I have got an equation of Sum.i want to split the eqution three parts of Sum,how can do?

ans:=sum(I*`ℏ`*(diff(c[n](t), t))*phi(r)*exp(-I*omega[n]*t)-(1/2)*E[0](e_.r_)*e*c[n](t)*phi(r)*omega[n]*exp(-I*omega[n]*t+I*Omega*t)-(1/2)*E[0](e_.r_)*e*c[n](t)*phi(r)*omega[n]*exp(-I*omega[n]*t-I*Omega*t), n = l .. k)

i want to  result like this

ans1:=sum(I*'ℏ'*(diff(c[n](t), t))*phi(r)*exp(-I*omega[n]*t), n = l .. k) = sum((1/2)*E[0...

Hey folks, I'll give you the actual problem I'm trying to solve before I show you my code in case anyone can think of a better idea...

Long read but you can ignore the examples and additional info if you don't need it.

Take the doubling map,

f(x) = 2x if 0 <= x <= 1/2

f(x) = 2x - 1 if 1/2 < x <= 1.

Or in maple code, f := x -> piecewise(x<=1/2, 2*x, 2*x - 1):

This looks...

I have a problem changing a maple matrix to a matlab code without getting the error

Error, (in PrintTarget) assigning to a long list, please use Arrays ( I don't know how to use arrays ) thanks Big in advance

Here is my code:

> J[p1] := Matrix(3, 4, {(1, 1) = -sin(theta[1](t))*`&Delta;x`[1]-cos(theta[1](t))*`&Delta;z`[1], (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (2, 1) = cos(theta[1](t))*`&Delta;x`[1]-sin(theta[1](t))*`&Delta;z`[1], (2, 2...

Integrals.mw

I tried solving the following equation but when i do i get

> solve((500*.87)*(1-(1/x)^.2385)-(288*(x^.3174-1))/(.85), x);

Warning, solutions may have been lost

Is it not possible to solve for highly nonlinear equations ?

i want to use numerical technique to solve non linear ode.

are there some built in numerical techniques in maple13.?

kindly name few and there syntax.

i have seen that we can solve linear ode numerically. how can i check that what technique does maple use?

I have 2 queestions:

1- concerned about the title font format:

fot the title

tit:=sprintf("W=%g,D=%g,f=%g,M=%g,N=%g,q=%a",Omega,Delta,phi,M,N,q);
                        "W=0.2,D=0,f=0,M=0,N=0,q=-10"

the little problem in the last value q=-10 it appear theta=10 even the sign (-) invisible

2- the sign(-) at the tickmars...

Error_redirect_pa.mw  i got the problem to get repeated value of unknown parameter ''ALPHA"" using pade approximation in maple. i uploted my worksheet kind watch out my problem. and give some suggessions...Error_redirect_pa.mw 

hello ,

I have figured out how to solve for corresponding values but when i try to solve the set of equations by using for example,

a :=[10 values]

b:=[10 values]

w:= [10 values]

##here i want to solve the set of equations for  [a1 b1 w1] ,[a2 b2 w1] , [a3 b3 w1]....like i want to solve the equations for all correponding values of a,b @ each w ##

set of equations (function of a,b,c)

i cant make it work...

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/ad_plot_for_derivati.mw .

Download ad_plot_for_derivati.mw

restart; with(Physics); Setup(mathematicalnotation = true):

am := Annihilation(N); ap := Creation(N)

Setup(op = {b[1], b[2], b[3]})

alg := seq(%Commutator(ap, b[i]) = 0, i = 1 .. 3), seq(%Commutator(am, b[i]) = 0, i = 1 .. 3), %Commutator(b[1], b[2]) = b[3], %Commutator(b[1], b[3]) = -Physics[`*`](2, b[1]), %Commutator(b[2], b[3]) = Physics[`*`](2, b[2])

q := %Commutator(Physics[`.`](ap+am, b[2]), b[1])

p := expand(q)

ans := sort(p, [a, b])

Hey folks I have a question about how to use solutions from fsolve in future calculations without copying and pasting.

Here is the code...

Order := 13:
f := x -> -x + x^2 - x^3 + add(a[i]*x^i,i=4..11):
A := n -> coeff(series(f(f(x)), x), x^n):
fsolve({seq(A(k) = 0, k=5..12)}, {seq(a[j], j=4..11)});

Now this gives values of a[k] from a[4] to a[11] but then if I were to type a[5] in the...

I got the range result as limit

ki:= t -> (sum(alpha1(S[j])*(exp(2*I*k*delta*t)-exp(S[j]*t))/((2*I*k*delta-S[j])*Q1(S[j])),j=1..3));limit(ki(t),t=infinity);
where : k,delta are constant and sum(alpha1([sj]),j=1..3) is number

the output is:

-.1294653014455612741589137389024249282843455188351695736221586639528345838460087707530271384645825938-0.8233085340292186435280681029574142108148445032573172243750416705060582056155292536631305980808392462e-1*I...

p := (g+Delta*g*Zeta/omega-g*Zeta)*a[1]*b[2]+(g+Delta*g*Zeta/omega-g*Zeta)*a[2]*b[1]+(-Delta*g*Zeta/omega+g-g*Zeta)*a[2]*b[2]+(-Delta*g*Zeta/omega+g-g*Zeta)*a[1]*b[1]

q := factor(add(op(i, p), i = 1 .. 2))

q := g*(omega+Delta*Zeta-Zeta*omega)*(a[1]*b[2]+a[2]*b[1])/omega

here,throuth factor ,simplify the expression .yet,if i want the expression like this:

(g+Delta*g*Zeta/omega-g*Zeta))*(a[1]*b[2]+a[2]*b[1])

how do i?

thanks

what can I do I stuck in it,

f:=c1*a+c2*b:

where c1 and c2 is constant of three terms

a,b are variables in x