here is a matrix

A := Matrix(2, 2, {(1, 1) = ((1/2)*I)*E[0]*d[ba]/`ℏ`, (1, 2) = -((1/2)*I)*delta, (2, 1) = ((1/2)*I)*delta, (2, 2) = ((1/2)*I)*E[0]*d[ab]/`ℏ`})

it is clear that there is a common coeff with each element by ((1/2)*I,i just want to show the matrix with traditional way.

#A:="((1/2)*I*Matrix(2, 2, {(1, 1) = -E[0]*d[ba]/`ℏ`, (1, 2) = delta, (2, 1) = -delta, (2, 2) = -E[0]*d[ab]/`ℏ`})"

A1:=((1/2)*I;

I have a complicated function (combination of exponential and sin), and I want to find the infinite integral of that from 0 to infinity which I know that the function itself converges to zero at infinity. I used evalf to force maple to do that but it fails to give the numerical answer. How can I resolve the problem?

Thanks,

Hamid

Maple uses the LGPL licensed GMP as a support library for doing large integer and high precision float computations. There exists is a later version of GMP than was bundled and shipped in Maple 14.x. A natural question is whether it would improve Maple's performance. Another question is whether it can be used within Maple 14.

I  have to solve this initial value problem:

pde := (-x+y)*(diff(u(x, y), x))+(-y*(diff(u(x, y), y)))+c*u(x, y)=0

subject to u(0,y)=g(y)

I have had a go and this is what I have so far:

a := -x(s)+y(s)

b := -y(s)

c := 1

pde := a*(diff(u(x, y), x))+b*(diff(u(x, y), y))+c*u(x, y) = 0

chareq1 := diff(x(s), s) = a

chareq2 := diff(y(s), s) = b

here is an equation

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

req := Int(conjugate(f[m](r))*rhs(eq), r)

leq := Int(conjugate(f[m](r))*lhs(eq), r)

sup:=Int(conjugate(f[m](r))*f[n](r), r) = delta[m, n]

ans:=applyrule(sup,leq=req)

how can get the result as

ans := I*`ℏ`*(diff(c[m...

Hello!

I have solved a set of ODEs. For example I have obtained L(t). Now, I have a relation between "t" and for example "z", such as "t = exp(z)". How can I get for example L(z=1)?

Thanks!

Download for_mapleprimes.mw

I used dolve for dsystem,

but I got this error:cannot evaluate right than .....

is there anyway to overcome this problem?

thanks

GCD := proc (a, b) local sol, x, y, eqs, vars; vars := [x, y]; eqs := [Typesetting:-delayDotProduct(cos(x), sin(y)) = 0, y+sin(x) = b]; sol := solve(eqs, vars); sol end proc

I tried to convert this code into C, using C(GCD) and got an error message that target language may not recognize commands like solve etc., which is obvious. Please let me know if there is any way to convert procedures involving such commands into C code.

I am using CodeGeneration[Fortran] to generate a Fortran subroutine.  The subroutine simply takes a vector state as input and generates a large matrix as output.  However, upon calling "Fortran" I have been getting the following warning that I have not been able to figure out:

"Warning, character limit for Fortran77 statements exceeded; please filter input through codegen[split] first."

If anyone can provide some tips, I would appreciate it.

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...

Hi, I have recently started working with Maple. I need some help with solve command.

Now if I run solve on above two equations for values...

I keep getting this error,

not matter what command. 

How can I solve this please?

We have been asked to write a general code for the revised simplex method in Maple.  We have been trying for several days but haven't got anywhere.  We think we have written (a very basic) code for the first iteration but we don't know how to do a loop so that it keeps going until the solution has been found. 

PS its a maximisation problem.

Using SVD to underst.mw

Hola ev1;
In the procedure PLOTSVD, I get the error
"Error, (in SVDPLOT) invalid input: linearalgebra" SingularValues Expects output value for keyword parameter to Be of type {list (identical (U, S, Vt, list)), identical (U, S, Vt, list)}