Here is the minimal working example: 

with(Physics):
Setup(noncommutativeprefix={W,V});
simplify(W^(-1)*(-W^(-1)-W*V^(-2)));

It gives the result  below, which is wrong. Or am I misinterpreting? (Maple 2018)

(-V^2-W^2)*(1/V^2)*(1/W^2)

how i can solve a system of integral equations? thanks for the help.
 

Download fsolve.mw

Hello,

I have a procedure, SoS that returns the value for sos. After compiling the procedure, I get float(undefined) for the return value. How can this be fixed?

Thanks so much for your time.

Download Compile_proc.mw

Hey, to all!

First post here and after searching the web for answers I thought I'll ask it here.

haven't used Maple for a long time so I'm quite new in it.

I've solved a PDE in a numeric and in a analytic way and I want to the find the difference between the results.

I've come up with this:

1.2.mw

Had few tries with a few commands but for no success.

I seem to be able to find the specific value of the analytical solution but not the numeric one although it is solved.

Also I can't subtract to find the difference.

I will appreciate help in any kind,

Thank you!

I'm not sure if it's possible to find a formula with predictive qualities, but I thought I'd ask anyway ...

IOD.mw

         I'd like to draw a graph with its  vertices  of 3d style.  Just like following image：

  I did not find the desired choices for this graph.

Hello All,

I have an autonomous system of ordinary nonlinear equations. In order to investigate it near the fixed points, I would like to reduce it to the normal or hyper-normal form. I see in the literature that some authors developed unique algorithms for the very specific differential systems (e.g., the systems with two/three eigenvalues), which is not my case.  

Maybe you know, if Maple has any specific commands that may conduct such type of reduction?

Thanks in advance,

Dmitry

Hi,

I’m trying to animate the volume generated by rotating the function curve around the Ox axis. the Solidofrevolution command with animation output generate a disk volume

Thanks

QuestionAnimationVolumeOfRevolution.mw

I am solving coupled PDEs under various BCs and ICs on the (x,t) domain, to obtain momentum p((x(t),t) solution plots of the within a MAPLE worksheet. 

However, I am getting JAVA memory overload (> 2GB), which brings MAPLE to a halt when plots are displayed in the worksheet.  Is there a way of writing the graphical plots directly to a file on my computer, so that the plots can later be view separately on my computer screen, thus avoiding JAVA memory overload within MAPLE.

If this is possible, please could someone send me MAPLE code to create external plots...perhaps using simple examples in MAPLE?

Melvin .

with(DEtools); with(plots); with(plots, implicitplot);

LC := [diff(x(t), t) = x(t)*(1-x(t)^2-y(t)^2)+y(t)*((1-x(t))^2+y(t)^2), diff(y(t), t) = y(t)*(1-x(t)^2-y(t)^2)-y(t)*((1-x(t))^2+y(t)^2)]

p1 := DEplot(LC, {x(t), y(t)}, t = 0 .. 100, [x(0) = 1, y(0) = 1])

Error, (in DEtools/DEplot/CheckInitial) too few initial conditions: [x(0) = 1]

Download mm.mwmm.mw

Hello,

I am trying to generate a compiled c function

dc_test := proc( a :: Array(datatype=float[8]), b :: Array(datatype=float[8]), c :: Array(datatype=float[8]) ) 

c[..,..] := a . b;

return 0

end:

cdc_test:=Compiler:-Compile(dc_test, optimize):

a := Array(1..2, 1..2, fill=1.0, datatype=float[8]);
b := Array(1..2, 1..2, fill=1.0, datatype=float[8]);
c := Array(1..2, 1..2, datatype=float[8]);
                           
dc_test(a, b, c);

cdc_test(a, b, c);

This doesn't work because it doesn't recognize the "." symbol. Also if I were to take the transpose, for example Transpose(a), then it doesn't recognize that either. Are there commands for this that I'm unware of or do you have to program them from scratch?

Hi,

How to animate a plot with Shadebetween command?

Thanks

AnimationShadebetween.mw

Hi, I am currently working with Maple on school, and I wanted to use the diff(); function to find the partial derivative with respect to x. Anyways; I noticed when I was not using the multiplication sign I got the answer 0 as seen below:

diff(2*sqrt(xyz), x) = 0

But when I used the multiplication sign as seen (x*y*z instead), I got the correct answer.

diff(2*sqrt(x*y*z), x) = y*z/sqrt(x*y*z)

Why so?

Thanks,

Adrian

hi daer 

i am trying to run the following but error came 

acually the boundry conditions are periodic

x=linspace(0,1,50);
t=linspace(0,1,50);
m=0;
eqn=@(x,t,u,dudx)transistorPDE(x,t,u,dudx,C);
ic=@(x)transistorIC(x,C);
sol=pdepe((m,eqn,ic,transistorBC,x,t);
function [c,f,s]=transistorPDE(x,t,u,dudx,C)
a=0.1;
c=1;
f=0;
s=-a*dudx;
end
function u0=transistorIC(x,C)
L=1;
u0=exp(-(x-L/2)^2);
end
function [pl,ql,pr,qr]=transistorBC(xl,ul,xr,ur,t)
pl=ul;
ql=0;
pr=pl;
qr=0;
end
 

Hello Everyone

I have solved the equation below numerically and now I want to integral p(x)

ode := diff((b+5)*((1/6)*(diff(p(x), x))*h^3+(1/2)*c*h^2)/(b+2)+De^2*((1/10)*(diff(p(x), x))^3*h^5+c^2*(diff(p(x), x))*h^3+(1/2)*c*(diff(p(x), x))^2*h^4+c^3*h^2)/(b+2), x)+diff(h, x) = 0

ics := p(0) = 0, p(1) = 0;
                   ics := p(0) = 0, p(1) = 0
solnum := dsolve({ics, ode}, numeric, method = bvp, maxmesh = 600);
                solnum := proc(x_bvp)  ...  end;
odeplot(solnum, 0 .. 1)

but when I write evalf(int(p,0..1)) I get this result:. What should I do to have the value?