Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I just wanted to let everyone know that the Call for Papers and Extended Abstracts deadline for the Maple Conference has been extended to June 14.

The papers and extended abstracts presented at the 2019 Maple Conference will be published in the Communications in Computer and Information Science Series from Springer. We welcome topics that fall into the following broad categories:

  • Maple in Education
  • Algorithms and Software
  • Applications of Maple

You can learn more about the conference or submit your paper or abstract here: 

https://www.maplesoft.com/mapleconference/Papers-and-Presentations.aspx

Hope to hear from you soon!

I describe here a finite difference scheme for solving the boundary value
problem for the heat equation

"(&PartialD; u)/(&PartialD; t)= ((&PartialD;)^)/((&PartialD;)^( )x^)(c(x)(&PartialD; u)/(&PartialD; x)) + f(x,t)   a<x<b,   t>0"

for the unknown temperature u(x, t)subject to the boundary conditions

u(a, t) = alpha(t), u(b, t) = beta(t), t > 0

and the initial condition

"u(x,0)=`u__0`(x),    a < x < b."

 

This finite difference scheme is designed expressly with the goal of avoiding

differentiating the conductivity c(x), therefore c(x) is allowed to be

nonsmooth or even discontinuous. A discontinuous c(x) is particularly
important in applications where the heat conduction takes place through layers
of distinct types of materials.

 

The animation below, extracted from the worksheet, demonstrates a solution 

corresponding to a discontinuous c(x).  The limit of that solution as time goes to

infinity, which may be calculated independently and exactly, is shown as a gray
line.

Download worksheet: heat-finite-difference.mw

 

i wnat to take this export numbers from for loops in the matrix
restart;
with(LinearAlgebra);
f(x) :=  3*x^2+1 :
for i by .5 to 3.5 do print(i, f(i)) end do;

same this pic

I want to get solutions (approximate solutions) of the inequality g(x) > 0. I tried. 
 

restart; 
fprime := x-> (x-1)*(x-2)^2*(x-3)*(x-4); 
f := unapply(simplify(int(fprime(x), x)), x); 
g := unapply(expand(f(3*x+1)-x^3+3*x), x); 
solve(g(x) > 0, x);

I got

Please, I need to use Maple to solve Euler-Bernoulli Beam on Pasternak Foundation using Homotopy Perturbation Method.

The governing equation is 

initial conditions are 

the boundary condition is 

The governing equation represents Euler-Bernoulli beam on a generalized Pasternak viscoelastic foundation under an
arbitrary distributed dynamic load. in which E, I , ρ ,A are the parameters of the beam, representing Young’s modulus of elasticity, moment of inertia, density and area of cross section, respectively. K,C and Gp are spring stiffness, damping coefficient, and shear coefficient of the foundation. Moreover, y(x, t) and F(x, t) are defined as the vertical deflection of the beam and the generic arbitrary dynamic loads, respectively, where the loads distribute along the x-axis and t is time.

 

I will appreciate anyone who can help me with a Maple solution.

 

Thank you.

Would anybody be able to explain to me why I had to rewrite this one baby-procedure that has an if-statement inside my larger procedure, so that I could get a later on if-statement to work? As best I can tell, they mean exactly the same thing.


 

Information that may be pertinent:

 - The entire (overall) procedure worked without throwing any warnings prior to my adding this new if-statement.
 - The if-statement itself worked perfectly fine on another worksheet, though I hadn't tested it inside a procedure.
 - Commenting out the if-statement but retaining the interior commands led to the sheet compiling again, showing that it was in fact the presence of the if-statement causing the issue. 
 - When I commented out the offending baby-procedure that started throwing the warning, the sheet once again compiled.

The added if statement:

if eval_size > 4 then Search(max(seq(eval_cont[i], i = 4 .. eval_size)), eval_cont); evals_pairlist := evals_pairlist[[1 .. 3, %]] fi:

The error I get, and its "positioning": 

Error, invalid expression for inline function

 

Hello,

How I can remove some problems in my code.

I wrote it before in  software, but in some lines, I have some problems with it.

How I can change it in maple to run well?


 

restart

M := 2; B := 1; C := 3; K := 6; f := 2*exp(-t); q0 := 1; qd0 := -1; qd20 := 1; alpha := 1/2; beta := 1/12; gammaa := 1/4; h := .1; Tf := 2

"N:=size(M,1)  T:=0:h:Tf;  lt:=length(T);  F:=subs(f,t,T);  F:=subs(t=T,f);    qd30:=M\(F(1)-B*qd20-C*qd0-K*q0);   q:=zeros(N,lt)  qd:=zeros(N,lt)  qd2:=zeros(N,lt);  qd3:=zeros(N,lt);  q(:,1):=q0  qd(:,1):=qd0  qd2(:,1):=qd20  qd3(:,1):=qd30"

Error, invalid function arguments

"N:=size(M,1)  T:=0:h:Tf;  lt:=length(T);  F:=subs(f,t,T);  F:=subs(t=T,f);    qd30:=M\(F(1)-B*qd20-C*qd0-K*q0);  q:=zeros(N,lt)  qd:=zeros(N,lt)  qd2:=zeros(N,lt);  qd3:=zeros(N,lt);  q(:,1):=q0  qd(:,1):=qd0  qd2(:,1):=qd20  qd3(:,1):=qd30"

 

"A:=M+alpha*h*B+gammaa*(h^2)*C+beta*(h^3)*K;  for n:=1:lt-1 do      D:=F(n+1)-B*(qd2(:,n)+(1-alpha)*h*qd3(:,n))-...          C*(qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n))-...          K*(q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n));      qd3(:,n+1):=A\D;      qd2(:,n+1):=qd2(:,n)+(1-alpha)*h*qd3(:,n)+alpha*h*qd3(:,n+1);      qd(:,n+1):=qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n)+gammaa*(h^2)*qd3(:,n+1);      q(:,n+1):=q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n)+...          beta*(h^3)*qd3(:,n+1);  end do"

Error, invalid function arguments

"A:=M+alpha*h*B+gammaa*(h^2)*C+beta*(h^3)*K;  for n:=1:lt-1 do      D:=F(n+1)-B*(qd2(:,n)+(1-alpha)*h*qd3(:,n))-...          C*(qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n))-...          K*(q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n));      qd3(:,n+1):=A\D;      qd2(:,n+1):=qd2(:,n)+(1-alpha)*h*qd3(:,n)+alpha*h*qd3(:,n+1);      qd(:,n+1):=qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n)+gammaa*(h^2)*qd3(:,n+1);      q(:,n+1):=q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n)+...          beta*(h^3)*qd3(:,n+1);  end do"

 

plot(T, q, 'r-o', 'linewidth', 1.5)

Error, (in plot) unexpected options: [r-o, , 1.5]

 

``


 

Download

 

How should wirte a while loop for solve nonlinear equations by newton raphson method

 Hi guys,
I am trying write a code for homotopy perturbation, i have already generated the polynomial as you can see, i have also  solve for concentration equation since is not couple. But i have a lot of error massages for temperature, velocity and induced magnetic field. can some one please go through the code?
 

NULL

restart

PDEtools[declare](f(x),theta(x),u(x),w(x), prime=x):

f(x)*`will now be displayed as`*f

 

theta(x)*`will now be displayed as`*theta

 

u(x)*`will now be displayed as`*u

 

w(x)*`will now be displayed as`*w

 

`derivatives with respect to`*x*`of functions of one variable will now be displayed with '`

(1)

N := 4:

NULL

NULL

f(x):=sum((p^(i))*f[i](x),i=0..N);

f[0](x)+p*f[1](x)+p^2*f[2](x)+p^3*f[3](x)+p^4*f[4](x)

(2)

theta(x) := sum(p^i*theta[i](x), i = 0 .. N);

theta[0](x)+p*theta[1](x)+p^2*theta[2](x)+p^3*theta[3](x)+p^4*theta[4](x)

(3)

``

u(x) := sum(p^i*u[i](x), i = 0 .. N);

u[0](x)+p*u[1](x)+p^2*u[2](x)+p^3*u[3](x)+p^4*u[4](x)

(4)

``

w(x) := sum(p^i*w[i](x), i = 0 .. N);

w[0](x)+p*w[1](x)+p^2*w[2](x)+p^3*w[3](x)+p^4*w[4](x)

(5)

HPMEq := (1-p)*(diff(f(x), `$`(x, 2)))+p*(diff(f(x), `$`(x, 2))-k1*(diff(f(x), x))-k2*f(x));

(1-p)*(diff(diff(f[0](x), x), x)+p*(diff(diff(f[1](x), x), x))+p^2*(diff(diff(f[2](x), x), x))+p^3*(diff(diff(f[3](x), x), x))+p^4*(diff(diff(f[4](x), x), x)))+p*(diff(diff(f[0](x), x), x)+p*(diff(diff(f[1](x), x), x))+p^2*(diff(diff(f[2](x), x), x))+p^3*(diff(diff(f[3](x), x), x))+p^4*(diff(diff(f[4](x), x), x))-k1*(diff(f[0](x), x)+p*(diff(f[1](x), x))+p^2*(diff(f[2](x), x))+p^3*(diff(f[3](x), x))+p^4*(diff(f[4](x), x)))-k2*(f[0](x)+p*f[1](x)+p^2*f[2](x)+p^3*f[3](x)+p^4*f[4](x)))

(6)

HPMEr := (1-p)*(diff(theta(x), `$`(x, 2)))+p*(diff(theta(x), `$`(x, 2))-k11*(diff(theta(x), x))+k12*(diff(u(x), x))^2+k13*(diff(w(x), x))^2+k14*theta(x));

(1-p)*(diff(diff(theta[0](x), x), x)+p*(diff(diff(theta[1](x), x), x))+p^2*(diff(diff(theta[2](x), x), x))+p^3*(diff(diff(theta[3](x), x), x))+p^4*(diff(diff(theta[4](x), x), x)))+p*(diff(diff(theta[0](x), x), x)+p*(diff(diff(theta[1](x), x), x))+p^2*(diff(diff(theta[2](x), x), x))+p^3*(diff(diff(theta[3](x), x), x))+p^4*(diff(diff(theta[4](x), x), x))-k11*(diff(theta[0](x), x)+p*(diff(theta[1](x), x))+p^2*(diff(theta[2](x), x))+p^3*(diff(theta[3](x), x))+p^4*(diff(theta[4](x), x)))+k12*(diff(u[0](x), x)+p*(diff(u[1](x), x))+p^2*(diff(u[2](x), x))+p^3*(diff(u[3](x), x))+p^4*(diff(u[4](x), x)))^2+k13*(diff(w[0](x), x)+p*(diff(w[1](x), x))+p^2*(diff(w[2](x), x))+p^3*(diff(w[3](x), x))+p^4*(diff(w[4](x), x)))^2+k14*(theta[0](x)+p*theta[1](x)+p^2*theta[2](x)+p^3*theta[3](x)+p^4*theta[4](x)))

(7)

HPMEs := (1-p)*(diff(u(x), `$`(x, 2)))+p*(diff(u(x), `$`(x, 2))-R*(diff(u(x), x))-A-k8*w(x)-k7*u(x)+k5*theta(x)+k6*f(x));

(1-p)*(diff(diff(u[0](x), x), x)+p*(diff(diff(u[1](x), x), x))+p^2*(diff(diff(u[2](x), x), x))+p^3*(diff(diff(u[3](x), x), x))+p^4*(diff(diff(u[4](x), x), x)))+p*(diff(diff(u[0](x), x), x)+p*(diff(diff(u[1](x), x), x))+p^2*(diff(diff(u[2](x), x), x))+p^3*(diff(diff(u[3](x), x), x))+p^4*(diff(diff(u[4](x), x), x))-R*(diff(u[0](x), x)+p*(diff(u[1](x), x))+p^2*(diff(u[2](x), x))+p^3*(diff(u[3](x), x))+p^4*(diff(u[4](x), x)))-A-k8*(w[0](x)+p*w[1](x)+p^2*w[2](x)+p^3*w[3](x)+p^4*w[4](x))-k7*(u[0](x)+p*u[1](x)+p^2*u[2](x)+p^3*u[3](x)+p^4*u[4](x))+k5*(theta[0](x)+p*theta[1](x)+p^2*theta[2](x)+p^3*theta[3](x)+p^4*theta[4](x))+k6*(f[0](x)+p*f[1](x)+p^2*f[2](x)+p^3*f[3](x)+p^4*f[4](x)))

(8)

HPMEt := (1-p)*(diff(w(x), `$`(x, 2)))+p*(diff(w(x), `$`(x, 2))-R*(diff(w(x), x))+k9*u(x)-k10*w(x))

(1-p)*(diff(diff(w[0](x), x), x)+p*(diff(diff(w[1](x), x), x))+p^2*(diff(diff(w[2](x), x), x))+p^3*(diff(diff(w[3](x), x), x))+p^4*(diff(diff(w[4](x), x), x)))+p*(diff(diff(w[0](x), x), x)+p*(diff(diff(w[1](x), x), x))+p^2*(diff(diff(w[2](x), x), x))+p^3*(diff(diff(w[3](x), x), x))+p^4*(diff(diff(w[4](x), x), x))-R*(diff(w[0](x), x)+p*(diff(w[1](x), x))+p^2*(diff(w[2](x), x))+p^3*(diff(w[3](x), x))+p^4*(diff(w[4](x), x)))+k9*(u[0](x)+p*u[1](x)+p^2*u[2](x)+p^3*u[3](x)+p^4*u[4](x))-k10*(w[0](x)+p*w[1](x)+p^2*w[2](x)+p^3*w[3](x)+p^4*w[4](x)))

(9)

for i from 0 to N do equ[1][i] := coeff(HPMEq, p, i) = 0 end do;

diff(diff(f[0](x), x), x) = 0

 

diff(diff(f[1](x), x), x)-k1*(diff(f[0](x), x))-k2*f[0](x) = 0

 

diff(diff(f[2](x), x), x)-k2*f[1](x)-k1*(diff(f[1](x), x)) = 0

 

diff(diff(f[3](x), x), x)-k2*f[2](x)-k1*(diff(f[2](x), x)) = 0

 

diff(diff(f[4](x), x), x)-k1*(diff(f[3](x), x))-k2*f[3](x) = 0

(10)

for i from 0 to N do equa[1][i] := coeff(HPMEr, p, i) = 0 end do;

diff(diff(theta[0](x), x), x) = 0

 

diff(diff(theta[1](x), x), x)-k11*(diff(theta[0](x), x))+k12*(diff(u[0](x), x))^2+k13*(diff(w[0](x), x))^2+k14*theta[0](x) = 0

 

diff(diff(theta[2](x), x), x)+2*k13*(diff(w[0](x), x))*(diff(w[1](x), x))-k11*(diff(theta[1](x), x))+2*k12*(diff(u[0](x), x))*(diff(u[1](x), x))+k14*theta[1](x) = 0

 

diff(diff(theta[3](x), x), x)+k12*(2*(diff(u[0](x), x))*(diff(u[2](x), x))+(diff(u[1](x), x))^2)+k14*theta[2](x)+k13*(2*(diff(w[0](x), x))*(diff(w[2](x), x))+(diff(w[1](x), x))^2)-k11*(diff(theta[2](x), x)) = 0

 

diff(diff(theta[4](x), x), x)+k12*(2*(diff(u[0](x), x))*(diff(u[3](x), x))+2*(diff(u[1](x), x))*(diff(u[2](x), x)))-k11*(diff(theta[3](x), x))+k14*theta[3](x)+k13*(2*(diff(w[0](x), x))*(diff(w[3](x), x))+2*(diff(w[1](x), x))*(diff(w[2](x), x))) = 0

(11)

for i from 0 to N do equat[1][i] := coeff(HPMEs, p, i) = 0 end do;

diff(diff(u[0](x), x), x) = 0

 

diff(diff(u[1](x), x), x)-R*(diff(u[0](x), x))-A-k7*u[0](x)+k5*theta[0](x)+k6*f[0](x)-k8*w[0](x) = 0

 

diff(diff(u[2](x), x), x)-R*(diff(u[1](x), x))-k7*u[1](x)+k6*f[1](x)-k8*w[1](x)+k5*theta[1](x) = 0

 

diff(diff(u[3](x), x), x)-R*(diff(u[2](x), x))+k6*f[2](x)-k7*u[2](x)+k5*theta[2](x)-k8*w[2](x) = 0

 

diff(diff(u[4](x), x), x)-R*(diff(u[3](x), x))+k5*theta[3](x)+k6*f[3](x)-k7*u[3](x)-k8*w[3](x) = 0

(12)

``

for i from 0 to N do equati[1][i] := coeff(HPMEt, p, i) = 0 end do;

diff(diff(w[0](x), x), x) = 0

 

diff(diff(w[1](x), x), x)-R*(diff(w[0](x), x))-k10*w[0](x)+k9*u[0](x) = 0

 

diff(diff(w[2](x), x), x)-k10*w[1](x)+k9*u[1](x)-R*(diff(w[1](x), x)) = 0

 

diff(diff(w[3](x), x), x)-k10*w[2](x)+k9*u[2](x)-R*(diff(w[2](x), x)) = 0

 

diff(diff(w[4](x), x), x)+k9*u[3](x)-R*(diff(w[3](x), x))-k10*w[3](x) = 0

(13)

con[1][0] := f[0](-1) = 1, f[0](1) = 1:

-.5000000000*k2+0.3435019841e-1*k2^4+.5000000000*k2*x^2-.2500000000*k2^2*x^2+.2083333333*k2^2+0.4166666667e-1*k2^2*x^4-0.2083333333e-1*k2^3*x^4+.1041666667*k2^3*x^2+0.4166666667e-1*k1^2*k2+0.5952380952e-3*k2^3*k1*x^7-0.8472222222e-1*k2^3+0.2480158730e-4*k2^4*x^8+0.4166666667e-2*x^6*k1^2*k2^2+0.8333333333e-2*k2*x^5*k1^3-0.9722222222e-2*k1*k2^3*x^5-0.3472222222e-1*k1^2*k2^2*x^4-0.2777777778e-1*k2*x^3*k1^3+0.5046296296e-1*k1*k2^3*x^3+0.6805555556e-1*k2^2*k1^2*x^2-0.4133597884e-1*k1*k2^3*x+0.1944444444e-1*k2*k1^3*x+1.+0.1388888889e-2*k2^3*x^6+0.1666666667e-1*k1*k2^2*x^5+0.4166666667e-1*k2*x^4*k1^2-.1111111111*k2^2*k1*x^3+0.9444444444e-1*k1*k2^2*x-0.8333333333e-1*k2*k1^2*x^2+.1666666667*k2*k1*x^3-0.3750000000e-1*k1^2*k2^2-0.6944444444e-3*k2^4*x^6+0.8680555556e-2*k2^4*x^4-0.4236111111e-1*k2^4*x^2-.1666666667*k1*k2*x

 

1-(1/2)*k2+(277/8064)*k2^4+(1/2)*k2*x^2-(1/4)*k2^2*x^2+(5/24)*k2^2+(1/24)*k2^2*x^4-(1/48)*k2^3*x^4+(5/48)*k2^3*x^2+(1/24)*k1^2*k2+(1/1680)*k2^3*k1*x^7-(61/720)*k2^3+(1/40320)*k2^4*x^8+(1/240)*x^6*k1^2*k2^2+(1/120)*k2*x^5*k1^3-(7/720)*k1*k2^3*x^5-(5/144)*k1^2*k2^2*x^4-(1/36)*k2*x^3*k1^3+(109/2160)*k1*k2^3*x^3+(49/720)*k2^2*k1^2*x^2-(125/3024)*k1*k2^3*x+(7/360)*k2*k1^3*x+(1/720)*k2^3*x^6+(1/60)*k1*k2^2*x^5+(1/24)*k2*x^4*k1^2-(1/9)*k2^2*k1*x^3+(17/180)*k1*k2^2*x-(1/12)*k2*k1^2*x^2+(1/6)*k2*k1*x^3-(3/80)*k1^2*k2^2-(1/1440)*k2^4*x^6+(5/576)*k2^4*x^4-(61/1440)*k2^4*x^2-(1/6)*k1*k2*x

 

2.400000000*k2+0.3589208394e-1*k2^4+1.104000000*k2^2+2.713333334*k1*k2+1.904000000*k1^2*k2+0.115520003e-1*k2^3+.939244445*k1*k2^2+.3973226666*k1^2*k2^2+0.1412642116e-1*k1*k2^3+.9218444444*k1^3*k2

(14)

NULL

"cond[1][0]:=theta[0](-1)=0.1, theta[0](1)=1,w[0](-1)=0, w[0](1)=0,u[0](-1)=0, u[0](1)=0:  for j from 1 to N do:  cond[1][j]:=theta[j](-1)=0, theta[j](1)=0,w[j](-1)=0, w[j](1)=0,u[j](-1)=0, u[j](1)=0:  end do:    for i from 0 to N do:  dsolve({equa[1][i],cond[1][i]},theta[i](x));  theta[i](x):=rhs(`%`):    end do:    theta(x):=evalf(simplify(sum(theta[n](x),n=0..N)));  convert(theta(x),'rational'); "

Error, (in dsolve) found the following equations not depending on the unknowns of the input system: {u[0](-1) = 0, u[0](1) = 0, w[0](-1) = 0, w[0](1) = 0}

 

theta[0](x)+theta[1](x)+theta[2](x)+theta[3](x)+theta[4](x)

 

theta[0](x)+theta[1](x)+theta[2](x)+theta[3](x)+theta[4](x)

(15)

``

"condi[1][0]:=theta[0](-1)=0.1, theta[0](1)=1,w[0](-1)=0, w[0](1)=0,u[0](-1)=0, u[0](1)=0,f[0](-1)=1, f[0](1)=1:  for j from 1 to N do:  condi[1][j]:=theta[j](-1)=0, theta[j](1)=0,w[j](-1)=0, w[j](1)=0,u[j](-1)=0, u[j](1)=0, f[j](-1)=0, f[j](1)=0:  end do:    for i from 0 to N do:  dsolve({equat[1][i],condi[1][i]},u[i](x));  u[i](x):=rhs(`%`):    end do:    u(x):=evalf(simplify(sum(u[n](x),n=0..N)))"

Error, (in dsolve) found the following equations not depending on the unknowns of the input system: {f[0](-1) = 1, f[0](1) = 1, w[0](-1) = 0, w[0](1) = 0, theta[0](-1) = 1/10, theta[0](1) = 1}

 

u[0](x)+u[1](x)+u[2](x)+u[3](x)+u[4](x)

(16)

``

"condit[1][0]:=theta[0](-1)=0.1, theta[0](1)=1,w[0](-1)=0, w[0](1)=0,u[0](-1)=0, u[0](1)=0,f[0](-1)=1, f[0](1)=1:  for j from 1 to N do:  condit[1][j]:=theta[j](-1)=0, theta[j](1)=0,w[j](-1)=0, w[j](1)=0,u[j](-1)=0, u[j](1)=0, f[j](-1)=0, f[j](1)=0:  end do:    for i from 0 to N do:  dsolve({equati[1][i],condit[1][i]},w[i](x));  w[i](x):=rhs(`%`):    end do:    w(x):=evalf(simplify(sum(w[n](x),n=0..N)))"

Error, (in dsolve) found the following equations not depending on the unknowns of the input system: {f[0](-1) = 1, f[0](1) = 1, u[0](-1) = 0, u[0](1) = 0, theta[0](-1) = 1/10, theta[0](1) = 1}

 

w[0](x)+w[1](x)+w[2](x)+w[3](x)+w[4](x)

(17)

NULL

``

``


 

Download completecode.mw

Hello,

I got some problems with command solve .

I have to solve 6 equations with 6 unknowns. The problem is that Maple solves this system in a wrong way - the values calculated by the command are bad. I would like to ask if there is any command in Maple while solving the equations that are responsible for possible approximations when solving equations.
In the attachment I will insert the code I used.

https://www.pastiebin.com/5ce6a1943f548

Hi freind,

Anybody knows what happen to this coding?

why this error happens? how can i fix it?

thank you!

Hi everyone,

I was wondering how I could modify the thickness of edges in a graph displayed as a Maple plot through DrawGraph. The point is, the graph comprises 100 vertices and 1000 edges. By default, edge thickness is set as 2 but due to the high number of edges I would like to set edge thickness to1 or even 0. How can I do so?

Here is an example with a random graph:

with(GraphTheory) : with(RandomGraphs) : G := RandomGraph(100, 1000) : DrawGraph(G)

 

Thank you very much

where

and .

Also,

 

MY TRY: You can download the code.mw

restart:
with(LinearAlgebra):
interface(rtablesize=20):
with(plots):
k:=2:M:=2:
Tm:=(t,m)-> 2*t*Tm(t,m-1)-Tm(t,m-2);
Tm(t,0):=1:
Tm(t,1):=t:

Tm2:=(tt,m)->subs(t=tt,Tm(t,m)): 
alpha:=(m)->piecewise(m=0,1/sqrt(Pi),sqrt(2)/sqrt(Pi));
psii:=(n,m,x)->piecewise((n-1)/(2^(k-1)) <= x and x <= n/(2^(k-1)),
alpha(m)*(2^(k/2))*Tm2(2^k*x-2*n+1,m), 0); 
psi:=(t)->Array([seq(seq(psii(i,j,t),j=0..M-1),i=1...2^(k-1))] ):
for i from 1 to ((2^(k-1))*M) do
r[i]:=evalf(psi((2*i-1)/((2^k)*M))):
end do:
m:=M*(2^(k-1)):
xi:=(i,n)->((i+1)^(n+1)-2*i^(n+1)+(i-1)^(n+1));
PB_n:=Matrix(m,m)

 

Could you help me to complete the code?

Best regards.

Good day,

I was recently using Maple 2019 for work on a project, and ran into an error. This error (which will be copied and pasted below for others to test) occurs when making assumptions across multiple lines (whether using the additionally function or not) while using IntegrationTools[Change]. It seems that, if during the process a variable that was within both the assumptions is subtracted from itself, the subtraction fails to happen and leaves what effectively equals 0 in the workings, making further workings impossible.

I'm wondering if anyone else is able to reproduce this error? I know the fix for it is to not disjoint the assumptions, but I am curious if others can easily reproduce it or if others have experiences with it!

As promised, below you will find my workings in order to reproduce this error!

Base Error:

restart;
assume(a>0,b>0,b>a,c>0,t>0);
interface(showassumed=0);
F := Int(sqrt(d-a*c^2*t),d=0...infinity);
assume(b>a);
IntegrationTools[Change](F,-a*t*c^2+d=-y,y)


Simple Fix:

restart;
assume(a>0,b>0,b>a,c>0,t>0);
interface(showassumed=0);
F := Int(sqrt(d-a*c^2*t),d=0...infinity);
IntegrationTools[Change](F,-a*t*c^2+d=-y,y)


Error Without Interface Change:

restart;
assume(a>0,b>a,c>0,t>0);
F := Int(sqrt(d-a*c^2*t),d=0...infinity);
assume(b>a);
IntegrationTools[Change](F,-a*t*c^2+d=-y,y)


Error When Using Additionally:

restart;
assume(a>0,b>a,c>0,t>0);
interface(showassumed=0);
F := Int(sqrt(d-a*c^2*t),d=0...infinity);
additionally(b>a);
IntegrationTools[Change](F,-a*t*c^2+d=-y,y)

 

restart;
lambda := 1;
                               1
mu := 1;
                               1
alpha := 1;
                               1
v := 2;
                               2
delta := 1;
                               1
m := 1;
                               1
d := 3;
                               3
l := 1;
                               1
omega := -(1/2)*a*lambda^2-a*m^2+2*a*mu*v-delta*l*m-(1/2)*delta*lambda^2+2*delta*mu*v-l^2*mu-(1/2)*lambda^2*mu+2*mu^2*v-2*a*mu-2*delta*mu-2*mu^2;
                            1      
                            - a + 1
                            2      
a[0] := ((2*d*v-2*d-lambda)*(1/2))*sqrt(2)*sqrt(gamma*(a+delta+mu))/gamma;
                    (1/2)                (1/2)
                 5 2      (gamma (a + 2))     
                 -----------------------------
                            2 gamma           
a[1] := 0;
                               0
a[2] := -(2*(1/2))*sqrt(2)*sqrt(gamma*(a+delta+mu))*(d^2*v-d^2-d*lambda+mu)/gamma;
                     (1/2)                (1/2)
                  7 2      (gamma (a + 2))     
                - -----------------------------
                              gamma            
Omega := lambda^2-4*mu*v+4*mu;
                               -3
H := (-lambda+sqrt(-Omega)*tan((1/2)*sqrt(-Omega)*xi))/(2*(v-1));
                  1   1  (1/2)    /1  (1/2)   \
                - - + - 3      tan|- 3      xi|
                  2   2           \2          /
u := a[0]+a[1]*(d+H)+a[2]/(d+H);
              (1/2)                (1/2)
           5 2      (gamma (a + 2))     
           -----------------------------
                      2 gamma           

                       (1/2)                (1/2)    
                    7 2      (gamma (a + 2))         
              - -------------------------------------
                      /5   1  (1/2)    /1  (1/2)   \\
                gamma |- + - 3      tan|- 3      xi||
                      \2   2           \2          //

eta := -l*y-m*x+omega*t;
                               /1      \  
                      -y - x + |- a + 1| t
                               \2      /  
u := a[0]+a[1]*(d+H)+a[2]/(d+H);
              (1/2)                (1/2)
           5 2      (gamma (a + 2))     
           -----------------------------
                      2 gamma           

                       (1/2)                (1/2)    
                    7 2      (gamma (a + 2))         
              - -------------------------------------
                      /5   1  (1/2)    /1  (1/2)   \\
                gamma |- + - 3      tan|- 3      xi||
                      \2   2           \2          //

f := diff(u, xi);
                                  /                      2\
       (1/2)                (1/2) |3   3    /1  (1/2)   \ |
    7 2      (gamma (a + 2))      |- + - tan|- 3      xi| |
                                  \4   4    \2          / /
    -------------------------------------------------------
                                                 2         
                  /5   1  (1/2)    /1  (1/2)   \\          
            gamma |- + - 3      tan|- 3      xi||          
                  \2   2           \2          //          
S := diff(u, xi);
                                  /                      2\
       (1/2)                (1/2) |3   3    /1  (1/2)   \ |
    7 2      (gamma (a + 2))      |- + - tan|- 3      xi| |
                                  \4   4    \2          / /
    -------------------------------------------------------
                                                 2         
                  /5   1  (1/2)    /1  (1/2)   \\          
            gamma |- + - 3      tan|- 3      xi||          
                  \2   2           \2          //          
xi := x+y-(-2*alpha*m-delta*l-delta*m-2*l*mu)*t;
                          x + y + 6 t
eq := (-omega-a*m*m-delta*l*m-mu*l*l)*u-(gamma*u*u)*u+(a+delta+mu)*S;
            /   (1/2)                (1/2)
/  3      \ |5 2      (gamma (a + 2))     
|- - a - 3| |-----------------------------
\  2      / |           2 gamma           
            |                             
            \                             

                 (1/2)                (1/2)          \         /
              7 2      (gamma (a + 2))               |         |
   - ------------------------------------------------| - gamma |
           /5   1  (1/2)    /1  (1/2)              \\|         |
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||         |
           \2   2           \2                     ///         \

     (1/2)                (1/2)
  5 2      (gamma (a + 2))     
  -----------------------------
             2 gamma           

                 (1/2)                (1/2)          \     
              7 2      (gamma (a + 2))               |     
   - ------------------------------------------------|^3 + 
           /5   1  (1/2)    /1  (1/2)              \\|     
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||     
           \2   2           \2                     ///     

                                                    /          
                          1                         |          
  ------------------------------------------------- |7 (a + 2) 
                                                  2 \          
        /5   1  (1/2)    /1  (1/2)              \\             
  gamma |- + - 3      tan|- 3      (x + y + 6 t)||             
        \2   2           \2                     //             

                              /                                 2
   (1/2)                (1/2) |3   3    /1  (1/2)              \ 
  2      (gamma (a + 2))      |- + - tan|- 3      (x + y + 6 t)| 
                              \4   4    \2                     / 

  \\
  ||
  ||
  //
value(%);
            /   (1/2)                (1/2)
/  3      \ |5 2      (gamma (a + 2))     
|- - a - 3| |-----------------------------
\  2      / |           2 gamma           
            |                             
            \                             

                 (1/2)                (1/2)          \         /
              7 2      (gamma (a + 2))               |         |
   - ------------------------------------------------| - gamma |
           /5   1  (1/2)    /1  (1/2)              \\|         |
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||         |
           \2   2           \2                     ///         \

     (1/2)                (1/2)
  5 2      (gamma (a + 2))     
  -----------------------------
             2 gamma           

                 (1/2)                (1/2)          \     
              7 2      (gamma (a + 2))               |     
   - ------------------------------------------------|^3 + 
           /5   1  (1/2)    /1  (1/2)              \\|     
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||     
           \2   2           \2                     ///     

                                                    /          
                          1                         |          
  ------------------------------------------------- |7 (a + 2) 
                                                  2 \          
        /5   1  (1/2)    /1  (1/2)              \\             
  gamma |- + - 3      tan|- 3      (x + y + 6 t)||             
        \2   2           \2                     //             

                              /                                 2
   (1/2)                (1/2) |3   3    /1  (1/2)              \ 
  2      (gamma (a + 2))      |- + - tan|- 3      (x + y + 6 t)| 
                              \4   4    \2                     / 

  \\
  ||
  ||
  //

simplify(%);
Error, (in simplify/tools/_zn) too many levels of recursion

1 2 3 4 5 6 7 Last Page 2 of 1544