I'm trying to solve a Quantum Mechanics problem that requires me to normalize a wavefunction before I proceed. I have the book solutions to these problems (I'm just trying to practice solving using Maple) so I know what the answers are SUPPOSED to be. But I'm not getting the answers I expected when dumping it into Maple.

What I put into Maple:

What I was SUPPOSED to get:

I assume this is a radical simplication I need to make. But in highlighting the expression, right clicking, and trying the different simplification options, nothing worked.

Thanks!

http://www.maplesoft.com/applications/view.aspx?SID=4229

from book example, it seems assumed that input size of data such as list size or matrix size is the same as

trained data set size, but this need to hard code infinite number of types of size

What is the method to programming neural network when input size is smaller or changing and not equal to size of trained data set?

Hello,

I have a list of equations. I would like to display this list in column.
Problem : as each term of my list is an equation, i can not transform my list in a vector.
How can I do to display a list of equations in column ?

Thank you for your help.

Hi guys,

I run my code and after a point all the numerical values are followed by a dot. Why is that? Also the symbolic variables have a 1. before them.

Thanks a lot!

Probably not Maple's territory, but maybe a challenge?  Can we can get maple to do this?  Done by Matlab found here..

http://stackoverflow.com/questions/28279337/matlab-cuda-ocean-wave-simulation

i got 2 curves
a := abs(x);
b := (3/4)*x^2+1/4;

how can i get the max distance between them from x = -1 until x =1?

 heeeeeeeeeelp
 

Hello! Hope everyone would be fine. I want to solve the following system of ODEs please help to find the numerical solution

N := .6; alpha := .4; beta := .1; Nt := .2; Pr := .5; Nb := .1; s := .2; lambda[1] := 1; delta := .5; gm := 1; Sc := 1:L:=1:

Eq1 := (alpha*s+1)*(diff(F(eta), eta, eta, eta))-(F(eta)+(1/2)*s*eta)*(diff(F(eta), eta, eta))+((1/2)*(diff(F(eta), eta))-s)*(diff(F(eta), eta))-2*(G(eta)^2-(1-gm)^2)-2*lambda[1]*(H(eta)+N*Y(eta))-(alpha+beta-(1/4)*delta*(diff(F(eta), eta, eta, eta)))*(diff(F(eta), eta, eta))^2-(alpha-2*beta)*(diff(F(eta), eta))*(diff(F(eta), eta, eta, eta))-(2*(alpha-beta-(1/4)*delta*(diff(F(eta), eta, eta, eta))))*(diff(G(eta), eta))^2-(2*(alpha-(1/4)*delta*(diff(F(eta), eta, eta))))*G(eta)*(diff(G(eta), eta, eta)) = 0; Eq2 := (alpha*s+1)*(diff(G(eta), eta, eta))-F(eta)*(diff(G(eta), eta))+G(eta)*(diff(F(eta), eta))+s*(1-gm-G(eta)-(1/2)*eta*(diff(G(eta), eta)))-(1/2)*alpha*s*eta*(diff(G(eta), eta, eta, eta))+((3/2)*alpha+beta)*G(eta)*(diff(F(eta), eta, eta, eta))-((1/2)*alpha+beta)*(diff(F(eta), eta))*(diff(G(eta), eta, eta))-delta*((diff(F(eta), eta, eta))^2+6*(diff(G(eta), eta))^2)*(diff(G(eta), eta, eta)) = 0; Eq3 := (diff(H(eta), eta, eta))/Pr-F(eta)*(diff(H(eta), eta))+(1/2)*H(eta)*(diff(F(eta), eta))-s*(2*H(eta)+(1/2)*eta*(diff(H(eta), eta)))+Nb*(diff(H(eta), eta))*(diff(Y(eta), eta))+Nt*(diff(H(eta), eta))^2 = 0; Eq4 := (diff(Y(eta), eta, eta))/Sc-F(eta)*(diff(Y(eta), eta))+(1/2)*Y(eta)*(diff(F(eta), eta))-s*(2*Y(eta)+(1/2)*eta*(diff(Y(eta), eta)))+Nt*(diff(H(eta), eta, eta))/Nb = 0;

IC1 := F(0) = 0, (D(F))(0) = 0, G(0) = gm, H(0) = 1, Y(0) = 1; IC2 := (D(F))(L) = 0, G(L) = 1-gm, (D(G))(L) = 0, H(L) = 0, Y(L) = 0; dsys1 := {Eq1, Eq2, Eq3, Eq4, IC1, IC2}; dsol1 := dsolve(dsys1, numeric, output = listprocedure, range = 0 .. L);

dsol1f := subs(dsol1, F(eta));

dsol1g := subs(dsol1, G(eta)); dsol1h := subs(dsol1, H(eta)); dsol1y := subs(dsol1, Y(eta));

With my best regards and sincerely.

Let:

f:=x->1/sqrt(2*Pi)*exp(-x^2/2);

I.e. f is a standard Gaussian PDF.

Then (in Maple 2016.1):

Int(convert(f(x)*f(y)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity):
evalf(%);

Returns:

1.692568751

However (again in Maple 2016.1):

int(convert(f(x)*f(y)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity):
evalf(%);

Returns:

-0.5641895835

This is clearly incorrect, as the integral of a positive function must be positive.

This also seems to be a problem in which ever version of Maple is used behind the scenes on this forum.

gives:

Hi, there

How can I find the recurrence relation  for second derivative of sequence of functions  f-{n}(x)=\frac{(1-x^2)^n}{n!} in  maple 15?

please specify the commands.

we know the solution f"_{n}(x)=2(1-2n)f_{n-1}(x)+4f_{n-2}(x)

Regards

M.R. Yegan

Dear all

I have created a script code in maple. I also have contructed a power circuit in matlab simulink. How I use my code in matlab?

Hi,

Seem to be a bit stuck. Here's my code:
 

Thanks in advance :-) 

It is suggested  

hypergeom([1/3, 2/3], [3/2], (27/4)*z^2*(1-z)) = 1/z

if z > 1. Here is my try to prove that with Maple:

Download simplification.mw

Maple provides efficient vectorization and automatic parallelization for many common operators. For example

x -> 2*~x*~cos~(x*~x)

But in my application it is common to want to create rather long vectorized operators starting from some complicated symbolic computations. Doing conversions by hand from symbolic expressions to element-wise operations is laborious and error prone.

As a very simple example consider that it is possible to obtain (almost) the same result as above by writing the following as a vectorized operation

D(x->sin(x^2))~

But there are at least two problems with this. First of all it is not nearly as efficient as the first operator and second, perhaps not unrelated, is that the datatype returned when applying this operator to a Vector/rtable of hardware floats (e.g. datatype=float[8]) becomes something  more general.

My question is how can I convert a complicated symbolic expression into an efficient numeric element-wise vector operation?

I have tried several different approaches but so far without success. In the case above for example it seemed natural to expect that the following derivative

D(x->sin~(x^~2))

would produce a vectorized result, but this is not the case. In another attempt I was unable to see how to perform substitions into an expression, e.g. like this

unapply(subs(`*`=`*`~, cos=cos~, diff(sin(x),x)), x)

I would be glad to receive suggestions and/or references to relevant documentation. 

I use the example procedure when search. Clock in help 

but elapsed function can only run one time

because it return clock is not running

need to run clock start again and calculate from beginning again

how elapsed function can run more times

Hi!

Everyone,

I want to draw  phase plane of system of three fractional order equations. 

Note that 

Also want the  phase portrait when the values of alpha are not same....

Also

Thanks