i want to solve an ode , but maple return an integral in result, how can i have an answer?

Download dsolve.mw

 Dear Maple users,

In Maple 18, I want to  label y axis as $\hat{\sigma}_y$.

labels = [x, sigma[y]] works fine but I have no idea how to put a hat on sigma[y] .

regards,

Using Tools -> Check for Updates, I have just tried to update my Maple installation from 2017.0 to 2017.1. The download itself went without a hitch, but then nothing happened. In view of the fact that I clicked a button titled 'Download and Install', I would certainly expect the installation to have started automatically after download. Are anyone else experiencing the same odd behaviour?

PS: Related question of mine is Automatic update of help pages?

Limit of function f(x)=[arcsin(x)]/tan(πx/2).  At point x=1

I can't take the inverse laplace transform of one kind of functions.

Would you like to help me?

Thanks.

hello,

i went to plot a complex numerical solution , i used odeplot but did not work

please help

Thank you very much 

Download mapleprime.mwmapleprime.mw

I created a loop to calculate, using random numbers, the probability of a die roll. I got 0, which I originally attributed to the fact that the probability might have just been that low, but in comparing it to a binomial distribution, it shouldn't be 0.

I tried editing the loop to do a coin toss instead, where the probability is 0.5, and that worked fine. So I don't think my loop structure is necessarily wrong. But perhaps I wrongly input getting a 0.167 probability somehow.

Any advice would be appreciated!

dierollfail.mw

https://drive.google.com/file/d/0Bxs_ao6uuBDUNmd5ZVJtX29GT3c/view?usp=sharing
https://drive.google.com/file/d/0Bxs_ao6uuBDUSUdHSTcwSEQtS3M/view?usp=sharing

would like to return K map of P1

Summation expression for logic only consider 1 but how about wildcard x ?

if consider wildcard x as 1 too, then will use below

source = [[0,0,1,0],[0,0,1,1],[0,1,0,1],[0,1,1,0],[0,1,1,1],[1,0,0,0],[1,0,0,1],[1,0,1,0],[1,0,1,1],[1,1,0,0],[1,1,0,1],[1,1,1,0],[1,1,1,1]];
i use Quine Mccluskey algorithm

got result below

Hi, friends

I need to calculate the double integral over a non-rectangular domain.

Say, the domain is the triangle (in red)

When I enter

int(int((x^2+y^2)*`if`((y-x-1/2 <= 0) and (y+2*x-2<=0) and (y+x/2-1/2>=0), 1, 0), x=0..1, numeric = true), y=0..1, numeric = true);

or 

int(int((x^2+y^2)*eval(`if`((y-x-1/2 <= 0) and (y+2*x-2<=0) and (y+x/2-1/2>=0), 1, 0)), x=0..1, numeric = true), y=0..1, numeric = true);

an error occur (Error, (in int) cannot determine if this expression is true or false: y-x <= 1/2 and y+2*x <= 2 and 0 <= y+(1/2)*x-1/2)

For me, it is desirable to write boundary conditions in the int operator itself, not as a separate expession.

Regarding my recent question http://www.mapleprimes.com/questions/221909-How-To-Extract-Data-From-Implicit-Function I would like to share an interesting observation. Here the code of the program:

restart;
R0 := ln(y)+Re(Psi(1/2+(2*(p^2+(1/2)*sqrt(2*I+4*ksi_fs^2*p^2)*tanh(sqrt(2*I+4*ksi_fs^2*p^2)*x)/(tau+0.5e-2*a)))/y))+gamma+2*ln(2)
tau:= 10.000:ksi_fs:=10:p:=0.037:
R0p:= unapply(R0, [a,x]):
R0f:= proc(a,x)
local r:= fsolve(R0p(a,x), y= 0..1);
   `if`(r::float, r, Float(undefined))
end proc:
M:= Matrix(
   (100,100),
   (i,j)-> R0f(i, 1 + (j-1)*(0.5-0)/(100-1)),
   datatype= float[8]
);

After approximately 2 hours of calculations I get a message window

But I repeat this calculations on another computer with the same Windows 7 64 bit and Maple 17 I don't get such error and I obtain desired data.

So can Maple be sensitive to the hardware? 

I don't like the font times new roman in Math 2D input. Can't find where to set a font anywhere except in the head at a new worksheet. But that alters immediately after executing a command. There should be a place in preferences at least to do a persistent change of font.

Hi,

I want to calculate the Jacobian of my function in terms of x,y,z. It is a Gerstner function and I need to caculate the normal of the displaced point. This can be achieved by getting the Jacobian Matrix. However when using Maple it isn't caculating it. The way I'm calling it is below. I have a vector function for tx,ty,tz:

with(VectorCalculus);
Jacobian([tx, ty, tz], [x, y]);

I'm not well versed in Maple I was wondering if anyone can help me out on how to get the Jacobian Matrix for tx,ty,tz

Best Regards

Paul

Hello everyone.

Please I am trying to obtain series expansion of the expression below in u and v up to order 30 but encounter difficulties cum maple is slow to display solution. Can I get help on the code and what to do to optimize the displayed time of maple?

Thank you in anticipation of your quick and positive responses and suggestions.

convert(series(convert(series((y[n]+((-8 h u^2 v^2-4 u^3 sin(u) h+2 sin(2 u) h u^3+2 sin(2 v) h v^3-4 v^3 h sin(v)+2 v^3 h sin(2 u+v)+2 u^3 h sin(u-2 v)+2 u^3 h sin(u+2 v)-2 v^3 h sin(2 u-v)-u^3 h sin(2 u+2 v)-v^3 h sin(2 u+2 v)-u^3 h sin(2 u-2 v)+v^3 h sin(2 u-2 v)+4 h u^3 v^2 sin(2 u)+4 h u^2 v^3 sin(2 v)-4 h u^3 v^2 sin(u-v)+4 h u^2 v^3 sin(u-v)-4 h u^3 v^2 sin(u+v)-4 h u^2 v^3 sin(u+v)+4 h u^2 v^2 cos(u)+4 h u^2 v^2 cos(2 u)+4 h u^2 v^2 cos(2 v)+4 h u^2 v^2 cos(v)-4 h u^3 v cos(2 u-v)-2 h u^2 v^2 cos(2 u-v)+2 h u v^3 cos(2 u-v)+4 h u^3 v cos(2 u+v)-2 h u^2 v^2 cos(2 u+v)-2 h u v^3 cos(2 u+v)+2 h u^3 v cos(u-2 v)-2 h u^2 v^2 cos(u-2 v)-4 h u v^3 cos(u-2 v)-2 h u^3 v cos(u+2 v)-2 h u^2 v^2 cos(u+2 v)+4 h u v^3 cos(u+2 v)+4 h u^3 v cos(u-v)+4 h u v^3 cos(u-v)-4 h u^3 v cos(u+v)-4 h u v^3 cos(u+v)+4 u sin(u) v^2 h-2 sin(2 u) h u v^2-2 sin(2 v) h u^2 v+4 v h sin(v) u^2-2 v u^2 h sin(2 u+v)-2 u v^2 h sin(u-2 v)-2 u v^2 h sin(u+2 v)+2 v u^2 h sin(2 u-v)+v u^2 h sin(2 u+2 v)+u v^2 h sin(2 u+2 v)-v u^2 h sin(2 u-2 v)+u v^2 h sin(2 u-2 v)) f[n])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((-8 h u^2 v^2+8 u^3 sin(u) h-4 sin(2 u) h u^3-4 sin(2 v) h v^3+8 v^3 h sin(v)-4 v^3 h sin(2 u+v)-4 u^3 h sin(u-2 v)-4 u^3 h sin(u+2 v)+4 v^3 h sin(2 u-v)+2 u^3 h sin(2 u+2 v)+2 v^3 h sin(2 u+2 v)+2 u^3 h sin(2 u-2 v)-2 v^3 h sin(2 u-2 v)+8 h u^3 v^2 sin(u)+8 h u^2 v^3 sin(v)-4 h u^3 v^2 sin(2 u+v)+4 h u^2 v^3 sin(2 u+v)+4 h u^3 v^2 sin(u-2 v)+4 h u^2 v^3 sin(u-2 v)+4 h u^3 v^2 sin(u+2 v)-4 h u^2 v^3 sin(u+2 v)-4 h u^3 v^2 sin(2 u-v)-4 h u^2 v^3 sin(2 u-v)+8 h u^2 v^2 cos(u)+8 h u^2 v^2 cos(v)+4 h u^3 v cos(2 u-v)-4 h u^2 v^2 cos(2 u-v)-8 h u v^3 cos(2 u-v)-4 h u^3 v cos(2 u+v)-4 h u^2 v^2 cos(2 u+v)+8 h u v^3 cos(2 u+v)-8 h u^3 v cos(u-2 v)-4 h u^2 v^2 cos(u-2 v)+4 h u v^3 cos(u-2 v)+8 h u^3 v cos(u+2 v)-4 h u^2 v^2 cos(u+2 v)-4 h u v^3 cos(u+2 v)-2 h u^3 v cos(2 u+2 v)+4 h u^2 v^2 cos(2 u+2 v)-2 h u v^3 cos(2 u+2 v)+2 h u^3 v cos(2 u-2 v)+4 h u^2 v^2 cos(2 u-2 v)+2 h u v^3 cos(2 u-2 v)-8 u sin(u) v^2 h+4 sin(2 u) h u v^2+4 sin(2 v) h u^2 v-8 v h sin(v) u^2+4 v u^2 h sin(2 u+v)+4 u v^2 h sin(u-2 v)+4 u v^2 h sin(u+2 v)-4 v u^2 h sin(2 u-v)-2 v u^2 h sin(2 u+2 v)-2 u v^2 h sin(2 u+2 v)+2 v u^2 h sin(2 u-2 v)-2 u v^2 h sin(2 u-2 v)) f[n+1])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((-8 h u^2 v^2-4 u^3 sin(u) h+2 sin(2 u) h u^3+2 sin(2 v) h v^3-4 v^3 h sin(v)+2 v^3 h sin(2 u+v)+2 u^3 h sin(u-2 v)+2 u^3 h sin(u+2 v)-2 v^3 h sin(2 u-v)-u^3 h sin(2 u+2 v)-v^3 h sin(2 u+2 v)-u^3 h sin(2 u-2 v)+v^3 h sin(2 u-2 v)+4 h u^3 v^2 sin(2 u)+4 h u^2 v^3 sin(2 v)-4 h u^3 v^2 sin(u-v)+4 h u^2 v^3 sin(u-v)-4 h u^3 v^2 sin(u+v)-4 h u^2 v^3 sin(u+v)+4 h u^2 v^2 cos(u)+4 h u^2 v^2 cos(2 u)+4 h u^2 v^2 cos(2 v)+4 h u^2 v^2 cos(v)-4 h u^3 v cos(2 u-v)-2 h u^2 v^2 cos(2 u-v)+2 h u v^3 cos(2 u-v)+4 h u^3 v cos(2 u+v)-2 h u^2 v^2 cos(2 u+v)-2 h u v^3 cos(2 u+v)+2 h u^3 v cos(u-2 v)-2 h u^2 v^2 cos(u-2 v)-4 h u v^3 cos(u-2 v)-2 h u^3 v cos(u+2 v)-2 h u^2 v^2 cos(u+2 v)+4 h u v^3 cos(u+2 v)+4 h u^3 v cos(u-v)+4 h u v^3 cos(u-v)-4 h u^3 v cos(u+v)-4 h u v^3 cos(u+v)+4 u sin(u) v^2 h-2 sin(2 u) h u v^2-2 sin(2 v) h u^2 v+4 v h sin(v) u^2-2 v u^2 h sin(2 u+v)-2 u v^2 h sin(u-2 v)-2 u v^2 h sin(u+2 v)+2 v u^2 h sin(2 u-v)+v u^2 h sin(2 u+2 v)+u v^2 h sin(2 u+2 v)-v u^2 h sin(2 u-2 v)+u v^2 h sin(2 u-2 v)) f[n+2])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((-6 u^2 h^2-6 v^2 h^2-4 cos(2 u) h^2 u^2 v^2-4 cos(2 v) h^2 u^2 v^2+8 v^2 u^2 h^2 cos(u-v)+8 v^2 u^2 h^2 cos(u+v)+8 u sin(u) v^2 h^2+8 sin(2 u) h^2 u v^2+8 sin(2 v) h^2 u^2 v+8 v h^2 sin(v) u^2+4 v u^2 h^2 sin(2 u+v)-4 u v^2 h^2 sin(2 u+v)+8 v u^2 h^2 sin(u-v)-8 sin(u-v) h^2 u v^2-8 v u^2 h^2 sin(u+v)-8 sin(u+v) h^2 u v^2+4 v u^2 h^2 sin(u-2 v)+4 u v^2 h^2 sin(u-2 v)-4 v u^2 h^2 sin(u+2 v)+4 u v^2 h^2 sin(u+2 v)-4 v u^2 h^2 sin(2 u-v)-4 u v^2 h^2 sin(2 u-v)-4 u v h^2 cos(2 u-v)+4 u v h^2 cos(2 u+v)-4 u v h^2 cos(u-2 v)+4 u v h^2 cos(u+2 v)-2 u v h^2 cos(2 u+2 v)+2 u v h^2 cos(2 u-2 v)+8 cos(u-v) h^2 u v-8 cos(u+v) h^2 u v-8 v^2 u^2 h^2+8 h^2 cos(u) u^2-2 cos(2 u) h^2 u^2+6 cos(2 u) h^2 v^2+6 cos(2 v) h^2 u^2-2 cos(2 v) h^2 v^2+8 h^2 cos(v) v^2-4 v^2 h^2 cos(2 u-v)-4 v^2 h^2 cos(2 u+v)-4 u^2 h^2 cos(u-2 v)-4 u^2 h^2 cos(u+2 v)+u^2 h^2 cos(2 u+2 v)+v^2 h^2 cos(2 u+2 v)+u^2 h^2 cos(2 u-2 v)+v^2 h^2 cos(2 u-2 v)) g[n])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((6 u^2 h^2+6 v^2 h^2+4 cos(2 u) h^2 u^2 v^2+4 cos(2 v) h^2 u^2 v^2-8 v^2 u^2 h^2 cos(u-v)-8 v^2 u^2 h^2 cos(u+v)-8 u sin(u) v^2 h^2-8 sin(2 u) h^2 u v^2-8 sin(2 v) h^2 u^2 v-8 v h^2 sin(v) u^2-4 v u^2 h^2 sin(2 u+v)+4 u v^2 h^2 sin(2 u+v)-8 v u^2 h^2 sin(u-v)+8 sin(u-v) h^2 u v^2+8 v u^2 h^2 sin(u+v)+8 sin(u+v) h^2 u v^2-4 v u^2 h^2 sin(u-2 v)-4 u v^2 h^2 sin(u-2 v)+4 v u^2 h^2 sin(u+2 v)-4 u v^2 h^2 sin(u+2 v)+4 v u^2 h^2 sin(2 u-v)+4 u v^2 h^2 sin(2 u-v)+4 u v h^2 cos(2 u-v)-4 u v h^2 cos(2 u+v)+4 u v h^2 cos(u-2 v)-4 u v h^2 cos(u+2 v)+2 u v h^2 cos(2 u+2 v)-2 u v h^2 cos(2 u-2 v)-8 cos(u-v) h^2 u v+8 cos(u+v) h^2 u v+8 v^2 u^2 h^2-8 h^2 cos(u) u^2+2 cos(2 u) h^2 u^2-6 cos(2 u) h^2 v^2-6 cos(2 v) h^2 u^2+2 cos(2 v) h^2 v^2-8 h^2 cos(v) v^2+4 v^2 h^2 cos(2 u-v)+4 v^2 h^2 cos(2 u+v)+4 u^2 h^2 cos(u-2 v)+4 u^2 h^2 cos(u+2 v)-u^2 h^2 cos(2 u+2 v)-v^2 h^2 cos(2 u+2 v)-u^2 h^2 cos(2 u-2 v)-v^2 h^2 cos(2 u-2 v)) g[n+2])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)),u=0,32),polynom),v=0,32),polynom);

Consider the following:

with(Physics):
Setup(anticommutativeprefix = psi):
psi^2,psi__1^2,psi__a^2;   # double underscores

Why does the square of psi__a (double underscore) not vanish as well? Or, perhaps, more to the point: why is psi__a not as well considered Grassmann-odd by Maple? Is this counter-intuitive behaviour intentional, or is it a bug?

I'm trying to numerically solve the differential equation: y' = -2xy + 1. Naturally, I come across the non-elementary integral of e^(x^2). By hand, I used a 2nd degree MacLaurin polynomial to get y = xe^(-x^2) + x^3/3e^(-x^2)+x^3/6e(-x^2). 

How do I use Maple to numerically solve this, with step sizes of h=0.1 and h=0.05 and plot them?