MaplePrimes Questions

Already searched and browsed multiple different threads and still cannot find a solution.

Apologizing the noob nature of this question.

 

 

In this code below, Why is the factor command not working?

f := a^2+x^2-2*ax;

a^2+x^2-2*ax

(1)

factor(f);

a^2+x^2-2*ax

(2)

expand((x-a)*(x-a));

a^2-2*a*x+x^2

(3)

``

NULL

NULL


Download factor_polynomial_2_multivariable.mwfactor_polynomial_2_multivariable.mw

 

Hello people in maple primes

I have a question, which is about the matrix shown in http://www.mapleprimes.com/questions/217852-HOW-I-Convert-Root-Of-In-To-Another-Common-Form

Why can't C below be shown with beta?

A := Matrix(3, 3, [[-a, a, 0], [0, 0, -sqrt(l*b*c*(j+k))/(j+k)], [2*j*sqrt(l*b*c*(j+k))/((j+k)*l), 2*k*sqrt(l*b*c*(j+k))/((j+k)*l), -c]]);
B:=subs(l*b*c*(j+k)=alpha,A);
C:=subs(j*alpha^(1/2) = beta,B);
e:=subs(alpha^(1/2) = gamma,B);

Best wishes.

taro

 

 

From differential expression I am trying to extract coefficients corresponding to given indets, but the coeffs command is giving me error as "invalid argument to coeffs"

For example for worksheet titled "[852] Determining Equations using Prolongation" the error is showing for expression (9).

Whereas the same command in worksheet titled "[411] Determining Equations for Variable Coeff KdV(Reviewed)" is working okay and is giving desired result at (16).

I am running both worksheets in Maple 18.

 

[852]_Determining_Equations_Using_Prolongation.mw

[411]_Determing_Equations_for_Variable_coefficient_KdV_(reviewed).mw

Regards

Can anybody please tell, where I can find seminar/workshops on Maple training ?

Regards

Dear All,

 

I am a new Maple user and I am still unaware of a lots of fancy features of Maple. I have a problem of simultaneous fitting polynomials. I wish that I could have help from you. Say, we have two polynomials of two variables,

f1(x,y)=a1+a2*x+a3*y+(a4+a5)*x2+(a4-a5)*y2;

f2(x,y)=b1+b2*x+b3*y+(a4-a5)*x2+(a4+a5)*y2.

Note that a4 and a5 are shared by the two polynomials. I would like to fit the two polynomials against their respective data set. Is there anyway I can do it using Maple? Any of your help is highly appreciated!

 

Best regards,

 

Toby

HI

HOW I convert root of in to another form

thanks

lorenz_5.mw

I'd like to find out what number Maple can still calculate the isprime for, and what number will make it go "No thanks, this is madness - ERROR".

I tried isprime (10100,000,000 + 2), and it can still calculate the result. No prime of course, that's the whole point of the +2.

I tried isprime(1010ˆ9 + 2) and it crashed.

But what is the tipping number?

I decided I should write a loop, and let maple write away the last number till it crashes. Is this possible? I think I'll need something like

> for n from 5 to ??? do isprime (10ˆn + 2) and ??? end do

This isn't suitable at all, but I cant figure out how to program "to error" or "write away this number"...

 

I hope the question is clear and look forward to ideas on how to make this work.

 

Dear Community,

 

I have difficulties evaluating a rather complicated infinite integral (see attachment). The result should be something like 12.319. Any idea how to go get this?

Tx for the kind help in advance,

best regards

Andras

Infinite_Integral_2.mw

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

Am here again, pls help me check out this adm code, cant get a result.

below is the attached file

 adm_2_method.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?

hi...how i can pdsolve this equation numerically or analyticlly?

this equation is time-fractional  equation with generalized Cattaneo model

where

 

 is the fractional derivative operator considered in the
Caputo sense.

 

FRACTION.mw

restart

k := 1; -1; rho := 1; -1; h := 1; -1; alpha := 2-Upsilon; -1; 0 < Upsilon and Upsilon <= 1

0 < Upsilon and Upsilon <= 1

(1)

k*(diff(T(z, t), z, z)) = rho*(diff(T(z, t), [`$`(t, alpha)]))

diff(diff(T(z, t), z), z) = diff(T(z, t), [`$`(t, 2-Upsilon)])

(2)

k*(diff(T((1/2)*h, t), z)) = 1:

k*(diff(T((-h)*(1/2), t), z)) = 0:

T(z, 0) = 0

T(z, 0) = 0

(3)

NULL



Download FRACTION.mw

 

Dear Community,

I would like to calculate a rather sophisticated function (see it in the attachmenet) Once I try to evaluate it at certain values I get the error message:

"Error (in pD) illegal use of an object as a name" 

What object Maple means? I guess its a minor mistake, yet I'm lost here, and the online help says "There is no help page available for this error"

Tx for the kind help in advance,

Andras

 

Bessel_plot.mw 


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

Dear Pros, I'm a biginer so I have a question about my program.

I have a lot of arrays which are result from 2 while loop. Now, I want creat a matrix from them but i can't. So, could you help me to do it.

For detail: 

V[1]:=[ 1 2 3]

V[2]:=[2 3 4]

V[3]:=[3 4 5]

V[4]:=[2 6 7]

V[5]:=[7 8 9]

...

V[n]

with type of V[i] is a array.

I searched and found a solution by manual to create a matrix as follow:

V_matrix:=<V[1],V[2],V[3]>

but in this case i can't but manual with n=100

please help me to have a Matrix.

Thank a lots.

First 1072 1073 1074 1075 1076 1077 1078 Last Page 1074 of 2428