Hellow ,

any body can help me to remove the error from the plot, here is my codes .

Note. 1. *Title {Effect of alpha when 'Gamma2'=Gamma2}

2. *on the y axis i want  bar(deltap)

restart;

h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
lambda:=unapply(Int(K1,z=0..1), Gamma2):
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plots:-display(Vector[row]([seq(plot([seq(eval(lambda(Gamma2), Nb=j), j=[0.1,0.2,0.3])], delta2=0.02..0.1, legend=[Nb=0.1,Nb=0.2,Nb=0.3],labels = ["d1 ", "Dp"],labelfont = ["SYMBOL", 12], title=(Effect of alpha when 'Gamma2'=Gamma2)), Gamma2=[10,20,30])]),axes=boxed);
 

How to proof classical and quantum both p summation equal to one in maple?

with(linalg):
A := matrix([[1,0],[0,1]]);
AA := exponential(t*A);

A := matrix([[-i/h,0],[0,-i/h]]);
AA := exponential(t*A);
Error, (in evalf/matrixexp) cannot determine if this expression is true or false: 0.9482523555e34*abs(t*i) < 1

arctan(sqrt(3)/2/(1/2));
cos(Pi/3) + i*sin(Pi/3);

exp(i*Pi/3);

cos(Pi/3) + i*sin(Pi/3) - exp(i*Pi/3);

with(ScientificConstants):
Constant(hbar);
h := evalf(%);
dsolve(diff(s(t),t) + i/h*s(t)=0,s(t));
dsolve(diff(s(t),t) + i*s(t)=0,s(t));
 

after solve differential equation, there is no Pi/3, 

can s(t) be a matrix ?

where do Pi/3 come from?

then i try Ket

with(Physics);
K1 := Ket(&psi;t);
dsolve(diff(K1,t) + i/h*K1=0,K1);

Error, invalid neutral operator;
Ket is matrix or complex number ?

and how to solve?

I am still buffled by a Record in Maple.

I simply want to make one variable be a Record. Then make separate variable, by making a COPY of the content of the first variable, then change the second variable. Sounds easy, right?

When I do this, I find that the changes in the second variable are being made also to the first variable.

I used eval() to make a copy of all field of first variable to the second.  May be this is not the correct way, but do not know how else (other by explicit copying one field one at a time).

I looked at help, and see nothing there. It talks about making the Record packed. I tried that, but it is still not working. Not a single example in help of how to copy the "content" of one record to another.

I simply want to make a new instance of the first variable, a copy. and change the second instance without the first changing at same time. Why is this so hard in Maple?

restart;
c:=Record('a','b'):
c:-a:=0:
c:-b:=-1:

c1:=eval(c): #make a copy of c??
print("c1=",eval(c1)); #yes, made a copy

c1:-b:=99:  #change one field  in the second variable

print("c1=",eval(c1));
print("c=",eval(c)); #why did c change?

I tried using Record[packed]('a' , 'b'): since help says something about same Record sharing same memory but that make no difference.

How to make a copy of a record to another in Maple without them being the same? I know I could do this

restart;
c:=Record('a','b'):
c:-a:=0:
c:-b:=-1:

c1:=Record('a','b'):
c1:-a:=c:-a;
c1:-b:=c:-b;

c1:-b:=99:  #change one field  in the second variable

print("c1=",eval(c1));
print("c=",eval(c)); #Now it did not change.

But my actual record has many fields. and I do not want to do the above each time.

aab := int(exp(-a*x)*cos(b*x), x=0..infinity);
simplify(eval(aab));
 

why above equation is not equal to a/(a^2+b^2) ?

and how to solve below equation

z assuming complex;
solve(z^n - z = 0, z);

it only return 1 or 0, but not complex equations in terms of n

I am trying to find formula of the sequence f(n)

f(n+1) = (1+f(n))/(2*f(n)), f(1) = -7/13

I tried
restart;
rsolve({f(1) = -7/13, f(n+1) = (1+f(n))/(2*f(n))}, {f})

But I couldn't get the result. With Mathematica, I got the result.

hi 

i cant comute the double integrate with infinity integral of my function and just evaluating

The calculation steps are too time consuming and will not be answered

Hi, I am very new to Maple and do therefore probably propose a very simple question, so I would kindly ask you to excuse me in advance.

Basically I want so solve a simple and discrete game by comparing expected utility payoffs. There are only three states of the world (0,1,2) and I must determine for what cost-price-ratio what quantitiy is chosen from the same set (0,1,2).

Mathematically it's easy and I solved it, but Maple does not let me present my findings in a way I would like, as presented in my picture (preview from LaTex):

To achieve this I would have to get all w over p and then define it as z (cost-price-ratio). But when I try solving for w/p or respectively z doesn’t work. Additionally, when I try to apply some rearrangements manually (like multiply by 1/p in the first inequality) it doesn't work either, I just receive a '*' in front of round brackets with my equation inside.

So, is there a possibilty to do these simple rearrangement of my inequalities in Maple? And is it possible to aggregate the solutions in the same way as I pictured it for my LaTeX code?

Thank you very much!

Well I also come to think of it, disagree with calling the double factorial function the double factorial function, i would have thought bi-factorial to be better. or singular parity factorial function. 

But my interface (2016)  queries me everytime i enter "n!!" as to whether I am referring to the double factorial function or the application of the factorial twice in iteration, ie :  (n!)!

I just think this was a bad decision to make common place and the notation cabal need to review and change this. 

Hi,

I m trying to work some Maple Formulas in the Excel 2016 environnement, but i m always getting Error Mesages ( See attachments)

Thanks

I need to move a term in the denominator, inside the radical sign in the numerator.  After some struggle, I can do it for a simple expression. But once the expression gets little more complicate, the method does not work.

Given expr:=sqrt(x*y)/x;  I want to force the "x" in the denominator inside the sqrt so I can simplify things. I can do this like this

restart;
expr1  := sqrt(x*y)/x;
expr1  := simplify(expr1,symbolic);
combine(expr1,symbolic);
algsubs(y/x=u,%);

But for the following similar expression, it does not work. Given expr:=sqrt(x^2+y^2)/x, I want to change to sqrt( (x^2+y^2)/x^2) or sqrt(1+ y^2/x^2) or sqrt(1 + u^2).  So the first step is as above, to force the x in denomiantor to go inside the radical. But the same steps do not work on this now

expr1:=sqrt(x^2+y^2)/x;
expr1:=simplify(expr1,symbolic);
combine(expr1,symbolic);

I tried radsimp(expr1); simplify(expr1,radical,symbolic); radnormal, I tried adding assuming positive everywhere. Nothing is working.

I am trying to convert an expression to be   (y/x) and then replace (y/x) to u.  BUt the first step is to get the x inside the radical and I am not sure how to do it for the second example and why it is failing. It is always hard for me in Maple to figure the right method to manipulate an expression.

Any suggestions?

I need to know if an expression contain a radical in it. For example, expr:=1/x*sqrt(x*y)  or expr:=3*(x*y)^(2/3)+5; or expr := sqrt(5), etc.. all these will return true, but for expr:=3+sin(x) it will return false.

I need to find if the expression has radical, becuase then I want to do special simplify if that is the case.

I do not want to use pattern matching, since the expression can be complicated, and thought to ask if there might be a simpler way, one of those hidden Maple commands which can do it.

I looked at  "has", "hasfun". with(PolynomialIdeals) has function IsRadical but it is only for polynomials.

Is there a way other than using patmatch?  

is possible to solve these partial differential equations in maple via pdsolve?

Thanks

Hello,

I tried to operate a transformation of some random variables but Maple doesn't compute the PDF of the outcome RV. Do you have a suggestion on how I should edit my code ? (the goal is to get the PDF of S)

Thanks in advance,

Antoine

Way i get this error

plot([x, cos(x), x = -Pi .. Pi, y = -.5 .. .5])

Error, (in plot) incorrect first argument [x, cos(x), x = -Pi .. Pi, y = -.5 .. .5]