kappa := Vector(7, [1,w[1]*(1-phi+phi*(1-1/(1+exp(-mu[p]-tau[p3]))))+(1-w[1])*
(1-phi+phi*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2])))),w[1]*phi/(1+exp(-mu[p]-tau[
p3]))+(1-w[1])*phi/(1+exp(-mu[p]-tau[p3]-eta[p2])),w[1]*(1-phi+phi*(1-1/(1+exp
(-mu[p])))*(1-phi)+phi^2*(1-1/(1+exp(-mu[p])))*(1-1/(1+exp(-mu[p]-tau[p3]))))+
(1-w[1])*(1-phi+phi*(1-1/(1+exp(-mu[p]-eta[p2])))*(1-phi)+phi^2*(1-1/(1+exp(-
mu[p]-eta[p2])))*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2])))),w[1]*phi^2*(1-1/(1+exp
(-mu[p])))/(1+exp(-mu[p]-tau[p3]))+(1-w[1])*phi^2*(1-1/(1+exp(-mu[p]-eta[p2]))
)/(1+exp(-mu[p]-tau[p3]-eta[p2])),w[1]*(phi/(1+exp(-mu[p]))*(1-phi)+phi^2/(1+
exp(-mu[p]))*(1-1/(1+exp(-mu[p]-tau[p3]))))+(1-w[1])*(phi/(1+exp(-mu[p]-eta[p2
]))*(1-phi)+phi^2/(1+exp(-mu[p]-eta[p2]))*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2]))
)),w[1]*phi^2/(1+exp(-mu[p]))/(1+exp(-mu[p]-tau[p3]))+(1-w[1])*phi^2/(1+exp(-
mu[p]-eta[p2]))/(1+exp(-mu[p]-tau[p3]-eta[p2]))]);

Download kappa.txt

Here is the expression, I am trying to simplify, given a set of rules. NEW_Cole.mw

I have tried different substitutions, using simplify with side rules, applyrule, eval, subs, algsubs.

But none seem to be working as the way I want them to be.

Is there a better way?

Thanks!

I am trying the following:

restart; 
with(Physics); 
with(Library); 
Setup(mathematicalnotation = true); 
Setup(op = A);
assume(n::integer);

type(A^n, Physics:-`^`(PhysicsType:-ExtendedQuantumOperator, integer));
type(A^n, PhysicsType:-ExtendedQuantumOperator^integer);

the result is

false
false

However, I expected to get true in both cases. How to construct correct type for A^n?
Thank you.

1) Consider the followin Integral:

int(cosh(a*x)*cos(b*x),x)

How can I tell maple to give the results in terms of hyperbolic and harmonic funtions?
The maple returns the answer in term of exponentials with imaginary powers.
Also, I could not use the "simplify" command to turn the answer returned by maple into an answer in terms of what I want! Please see the code below.


2) Is there a way that I can update the maple integration database manaully?

Example.mw

Hi,

I am new to the Maple Forum and have the following question.I am calculating some Gaussian error propagations and other equations and did this in this way:

[eval(diff(x(y,z),y),[y=...,z=...])²*(dy)²+eval(diff(x(y,z),z),[y=...,z=...])²*(dz)²]^1/2

Now I am having several dozens of different y and z for which I have to calculate the above equation. How can I define these variables y and z only once at the beginning before all equations and press the button !!! that everything afterwards is just calculated with the defined variables? This way I would only have to change the variables once instead of several times in each term and for each equation.

Thank you very much. I appreciate your help!

L

Is it possible to evaluate a function at multiple points described by an array or something of that sort and have Maple return the evaluations as an array. I need approximations of a function at various values of its argument so it would be nice to do it with a single command.

Thanks

How do I make "simplify" function to recognize complex exponentials. Maple created some exponentiol functions with a square root of negative expression in the denominator of the exponent. It refuses to simplify it. If I manually edit it, by taking negative sign under the root and adding immaginary number, the simplification works. How do I simplify it without manual edit?

The Eq 33 is manually edited line. Simplify function correctly uses cos function.

The Eq 34 is original expression. I tried to tell maple to assume real numbers hoping that it will wotk, but it did not.

I have question about how to write the expression to get the order in the sequence?

For exmple, if I have this sequence,

E:=[0,2,7,15,26,40]

I want to write a procedure that tells me the location of any input,x in the sequence E.

For example, if I choose x:=3, then the number 3 is actually between 2 and 7, so it will give me the answer C(x)=3 (because 7 is the third elemet).

I try to write the procedure but there always wrong! I do not know what expression should I write!

C:=proc(x)
local
for i from 1 to nops(E) do
if x<=E[i] then
c:=E[i]
fi;od;
end;

Please help! 

hi

I need your help,it is very necessary.I want to write two dimensional arrays to a file ,and then I should be able to read it again.I appericiate your help.

How do i use d'Alembert formula to solve,plot and animate with Maple software

Hi!

I seem to run into problems with (quantum) perturbation theory. In the following minimal working example, an energy denominator, as occuring in perturbation theory is not evaluated if I previously assume that the quantum number is a positive integer. It's supposed to return an error, as the energy denominator is 0.

Download energy_denominator.mw

Best regards,

Sören

how to calculate the cobordism between two ideals

Hi,

I was helping for a computing session today.

Maple 17 and 18 were not plotting the following:

restart:

with(plots):
spacecurve([x,cos(x), sin(x)], x = 0..10);

On some of the machines, with campus site license, windows 7, x64.

The plot came up with a "frame" only. And nothing was plotted.

We really can't figure out why?

Something to do with the settings(options)?

Had anyone had a similar experience?

Thanks,

casper

Find a series solution of the differental equation

y'=2*x*y

do graph have invariant for invariant of cayley graph

for example, how to convert below word relation into graph

a^2 = 1

a*b = b*a