Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hello, 

I have a PDE system. When I use pdsolve it gets me the messege " pdsolve->Warning: System is inconsistent". Is there a way I can see which equations breaks the system down? 
For this system, it's difficult to see from ayeball where the problem is. 
Thank you! 

test.mw

When I use the Determinant function on a matrix with (single variable) polynomial entries with real coefficients I often get an incorrect answer. I know the answers are incorrect because they have a higher degree or a lower lowest degree than is possible given the matrix elements.

However, when I replace the coefficients in the polynomials with rational numbers or I put in the option method=minor, I get the correct answer.

The problem seems to be roundoff error. However, the important error is not simply small changes in the resulting polynomial. The important error is the presence of entirely incorrect powers of the variable and not with very small coefficients.

How does this happen and why does the help page for Determinant( ) not warn of this behavior? In particuiar, why does the help page not say that using Gaussian elimination (i.e., the default) will often give incorrect answers in such cases, but using method=minor will work? Is this behavior known? I cannot find any reference to it on the internet.

 

Hello all. Is there any solution for the indefinite integralBadIntegral.mw

``

``

``

``

int((B*alpha[1]+(1/8)*`Νu`)*HeunT((1/16)*alpha[1]*(8*B*alpha[1]+`Νu`)*3^(2/3)*(2*alpha[2]+1)^2/(alpha[2]*(alpha[1]*alpha[2]*(8*B*alpha[1]+`Νu`))^(1/3)), 0, (1/2)*alpha[1]*(8*B*alpha[1]+`Νu`)*3^(1/3)/(alpha[1]*alpha[2]*(8*B*alpha[1]+`Νu`))^(2/3), (1/3)*3^(2/3)*(alpha[1]*alpha[2]*(8*B*alpha[1]+`Νu`))^(1/6)*y), y)

int((B*alpha[1]+(1/8)*`Νu`)*HeunT((1/16)*alpha[1]*(8*B*alpha[1]+`Νu`)*3^(2/3)*(2*alpha[2]+1)^2/(alpha[2]*(alpha[1]*alpha[2]*(8*B*alpha[1]+`Νu`))^(1/3)), 0, (1/2)*alpha[1]*(8*B*alpha[1]+`Νu`)*3^(1/3)/(alpha[1]*alpha[2]*(8*B*alpha[1]+`Νu`))^(2/3), (1/3)*3^(2/3)*(alpha[1]*alpha[2]*(8*B*alpha[1]+`Νu`))^(1/6)*y), y)

(1)

``

``

``

``

``

``

``

 

Download BadIntegral.mw

?

 

Thanks

i copy maple code from notepad to maple in maple window,

there is no error

my function in the code

explicit define parameters are Local type

for example

appendto("...");

func1(aaa)

Local aaa;

 

but when i run cmaple to read the code text file in window 8

it return error

missing operator, syntax error

at Local aaa;

originally 

i have defined

Local aaa, ii;

for ii from 1 to nops(aaa) do

etc.

but it has error too,

then i change to one by one

Local aaa;

Local ii;

still have error at Local aaa;

if i m working in maple 13 i have to solve a non linear integral equation then what will be the steps to use the do loop.

 

 

Hello all,

 

Is here anyone who could help me with this Integral.

I want to determine P(y=2) in which the Nu is a constant and also is unknown.

In deed I want to derive the P which only has one unknown parameter(Nu). But I cannot solve this indefinite integral although the \alpha_1 , \alpha_2 and B are known constants. The Nu parameter mustn't be valued preferely.Indefinite_INtegral.mw

 


Here is my code. I am deriving soil compression line for oedometric testing with porosity, below is a part of the derivation.

I would like to differentiate the (Sst/Vst)/(Spt/Vpt) by dsm first then integrate it with dsm ranges from 0 to dmax to get my final answer...

Hope someone can help me on writing the code... Thank you very much!!!

restart

``

Sst := `βss`*[dsx(dsm)^(2-Ds)-dsm^(2-Ds)]/(2-Ds);

`βss`*[dsx(dsm)^(2-Ds)-dsm^(2-Ds)]/(2-Ds)

(1)

Spt := `βps`*[dpx(dsm)^(2-Dp)-dpm(dsm)^(2-Dp)]/(2-Dp);

`βps`*[dpx(dsm)^(2-Dp)-dpm(dsm)^(2-Dp)]/(2-Dp)

(2)

Vst := `βsv`*[dsx(dsm)^(3-Ds)-dsm^(3-Ds)]/(3-Ds);

`βsv`*[dsx(dsm)^(3-Ds)-dsm^(3-Ds)]/(3-Ds)

(3)

Vpt := `βpv`*[dpx(dsm)^(3-Dp)-dpm(dsm)^(3-Dp)]/(3-Dp);

`βpv`*[dpx(dsm)^(3-Dp)-dpm(dsm)^(3-Dp)]/(3-Dp)

(4)

``

(dpm(dsm)/dpx(dsm))^(2-Dp) = (dsm/dsx(dsm))^(2-Ds);

(dpm(dsm)/dpx(dsm))^(2-Dp) = (dsm/dsx(dsm))^(2-Ds)

(5)

dpx := proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds) end proc;

proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds) end proc

(6)

dpm := proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds) end proc;

proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds) end proc

(7)

eval(Sst*Vpt/(Vst*Spt));

`βss`*[dsx(dsm)^(2-Ds)-dsm^(2-Ds)]*(3-Ds)*(2-Dp)*`βpv`*[(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds))^(3-Dp)-(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds))^(3-Dp)]/((2-Ds)*`βsv`*[dsx(dsm)^(3-Ds)-dsm^(3-Ds)]*`βps`*[(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds))^(2-Dp)-(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds))^(2-Dp)]*(3-Dp))

(8)

``

``

``

``


Download GEO_Assignment_3.mw

n := 5:
z1 := exp(2*3.14*I*k1/n)*cosh(z)^(2/n);
z2 := exp(2*3.14*I*k2/n)*sinh(z)^(2/n);
xx := Re(z1);
yy := Re(z2);
uu := cos(alpha)*Im(z1) + sin(alpha)*Im(z2);

i find that the 3d graph has many intersection points to itself

how to find these intersection points of calabi yau ?

 

 

I need to solve an ode of the type ay''+by'+cy=f(x) using cubic b spline.

can any one help me with the code or algorithm. Thanks

I need  some examples s.t. the computation of their lexicographic Groebner basis is heavy?

Thank you so much.

Plz help me! How to accelerate the calculation speed the following the program?

restart;
with(PDEtools); declare((u, W)(x, t)); KN := 10;
AFP := proc (C1, C2, C3, C4, H, KN, N) local ADM1, n, lambda, F, i, A0, A, U, W, u, L, R, NL, w, PDE1, IC1, d, Eq1, Eq2, Eq3, LI, trL, tr1, trN, Apr, AprK, trSol, Sol, AD; declare((u, W)(x, t)); ADM1 := proc (n) options operator, arrow; convert(subs(lambda = 0, value((Diff(F(Sum(lambda^i*U[i], i = 0 .. n)), `$`(lambda, n)))/factorial(n))), diff) end proc; A0[0] := F(U[0]); for n to KN do A0[n] := ADM1(n) end do; for n from 0 to KN do A[n] := convert(C1*(diff(subs({seq(U[i] = W[i](x, t), i = 0 .. KN)}, expand(subs(F(U[0]) = U[0]*exp(U[0]), A0[n]))), x)), diff)+convert(C2*subs({seq(U[i] = Diff(W[i](x, t), x), i = 0 .. KN)}, expand(subs(F(U[0]) = U[0]^2, A0[n]))), diff) end do; L := proc (w) options operator, arrow; diff(w(x, t), t) end proc; R := proc (w) options operator, arrow; C3*(diff(w(x, t), x)) end proc; NL := proc (w) options operator, arrow; C1*(diff(w(x, t)*exp(w(x, t)), t))-C2*(diff(w(x, t), x))^2 end proc; PDE1 := proc (w) options operator, arrow; L(w)-R(w) = -NL(w) end proc; IC1 := u(x, 0) = sum(2*(int(sin((d+1/2)*Pi*x/H), x = 0 .. H))*exp(-C4*(d+1/2)^2*Pi^2*t/H^2)*sin((d+1/2)*Pi*x/H)/H, d = 0 .. N); LI := proc (w) options operator, arrow; Int(w(x, t), t = 0 .. t) end proc; tr1 := u-rhs(IC1); Eq1 := LI(lhs(PDE1(u))) = LI(rhs(PDE1(u))); Eq2 := simplify(subs(lhs(Eq1) = tr1, Eq1)); trL := u = add(u[j](x, t), j = 0 .. KN); trN := LI(NL(u)) = Int(Sum(A[i], i = 0 .. KN), t = 0 .. t); Eq3 := subs(trL, lhs(Eq2)) = subs(trN, rhs(Eq2)); Apr[0] := u[0](x, t) = rhs(IC1); AprK := u[k+1](x, t) = -(Int(AD[k], t = 0 .. t)); for i from 0 to KN do Apr[i+1] := value(subs({seq(Apr[m], m = 0 .. i)}, subs({seq(W[m] = u[m], m = 0 .. i)}, subs(k = i, AD[i] = A[i], AprK)))) end do; trSol := {seq(Apr[i], i = 0 .. KN)}; value(subs(trSol, trL)) end proc;


 

Download AFP.mw

Ok i have used a diverse range of programming IDE now, but i am of course alot more comfortable the maple interface simply because of the amount of time i have spent on it, also it was my first "real" introduction to programming, but i was keen to trial and review of all the major CAS out on the market, but for a common purposed package within them they provide, for example the number theory package.

 

So the question is quite broad, im looking for suggestions for the most unbias way to go about this.

it would need to be a coauthored review with a representative writer for the seasoned user for each one right?


Is anyone running the 2016.1a patch on ubuntu 16.04?

A fresh installation starts up and works, but after I install the patch the installation seems broken:

cd /usr/local/maple2016/bin
./maple
/usr/local/maple2016/bin.X86_64_LINUX/mserver: error while loading shared libraries: libicuucmpl.so.56: cannot open shared object file: No such file or directory
maple: could not start mserver

it looks like my installation does not have libicuucmpl.so.56. I can see these however:
./bin.X86_64_LINUX/libicuucmpl.so.49.1.1
./bin.X86_64_LINUX/libicuucmpl.so
./bin.X86_64_LINUX/libicuucmpl.so.49

running xmaple has maple open, then stay unresponsive.

Here is the upgrade file, i've been trying to use:
md5sum ./Maple2016.1LinuxX64Upgrade.run
7502caaa65cc623d5d2574823eee9343  ./Maple2016.1LinuxX64Upgrade.run


Has anyone else encountered this problem?
thanks

First 133 134 135 136 137 138 139 Last Page 135 of 2097