Advanced Search

Maple 18 Questions and Posts

These are Posts and Questions associated with the product, Maple 18
View unanswered posts
Maple 18 Questions and Posts Feed

Why fsolve does not work in this case?...

Asked by: Hanskira 5
rootfinding fsolve incomplete-question
January 26 2023
1  5

I started to use maple recently and I do not know why fsolve does not work in this case:

Sol-NLSE.mw

I am wondering if the conjugate(x[j]) is the reason.

Thanks in advance

Maple Code to solve Perturbation Iteration Solutio...

Asked by: abiodunowoyemi9 30
January 14 2023
1  0

Hi,

May  I please get a Maple code to solve the below Perturbation Iteration Solutions for Volterra Integral Equations? 

I shall be very much grateful if I can get the Maple Code.

And possible plotting

Thanks

Problem with solving Pade approximant algorithm....

Asked by: Saha 75
plotting mtm incomplete-question
January 04 2023
0  2

 How to rectify the Error, (in MTM:-findsym) invalid input: a string/name list is expected for sort method `lexorder`   .Could you please help me to plot the curve.

TC1.mw

How to do this?...

Asked by: Saha 75
plotting
December 27 2022
0  2

SS.mw

How to rectify this error,I want a downfall curve ,Please help.

How to get theta(2),theta(3),...?...

Asked by: Saha 75
syntax homework user-error
December 22 2022
0  5

Thermal_Stress_Distribution.mw

Could you please help me to rectify this error.
> 

restart
> 

with(PDEtools):
> 

"Theta[0]:=1:Theta[1]:=Gamma:delta(0):=1:for w from 1 to 30 do delta(w):=0;od:delta(-1):=0:beta:=2:N[cc]:=1:N[rc]:=0.5:Q[int]:=0.2:gamma:=0.3:Mu:=0.5:Theta[a]:=0.2: "

Error, attempting to assign to `gamma` which is protected.  Try declaring `local gamma`; see ?protect for details.

  
> 

> 

for k from 0 to 8 do Theta(v+2) := solve(sum(delta[w-1]*(v-w+1)*(v-w+2)*Theta[v-w+2], w = 0 .. v)-Q[int]*gamma*Theta[a]*delta[v]-Q[int]*gamma*delta[v-1]*Theta[a]+M*Theta[a]*delta[v]/(1-Theta[a])+M*Theta[a]*delta[v-1]/(1-Theta[a])-M*(sum(delta[w-1]*Theta[v-w], w = 0 .. v))/(1-Theta[a])-M*Theta[v]/(1-Theta[a])+N[rc]*Theta[a]^4*delta[v-1]+Q[int]*gamma*(sum(delta[w-1]*Theta[v-w], w = 0 .. v))+Q[int]*gamma*delta[v]+N[cc]*Theta[a]*delta[v-1]+N[rc]*Theta[a]^4*delta[v]+N[cc]*Theta[a]*delta[v]+(v+1)*(v+2)*Theta[v+2]+Q[int]*delta[v-1]+Q[int]*delta[v]-N[rc]*(sum(sum(sum(sum(delta[l-1]*Theta[w-l]*Theta[Zeta-w]*Theta[w-Zeta]*Theta[v-w], l = 0 .. w), w = 0 .. Zeta), Zeta = 0 .. w), w = 0 .. v))-N[rc]*(sum(sum(sum(Theta[w]*Theta[Zeta-w]*Theta[w-Zeta]*Theta[v-w], w = 0 .. Zeta), Zeta = 0 .. w), w = 0 .. v))-N[cc]*(sum(delta[w-1]*Theta[v-w], w = 0 .. v))+beta*(sum(sum(delta[Zeta-1]*(w-Zeta+1)*Theta[w-Zeta+1]*(v-w+1)*Theta[v-w+1], Zeta = 0 .. w), w = 0 .. v))-N[cc]*Theta[v]+beta*(sum(Theta[v-w]*(w+1)*Theta[w+1], w = 0 .. v))+(v+1)*Theta[v+1]+beta*(sum((v-w+1)*Theta[v-w+1]*(w+1)*Theta[w+1], w = 0 .. v))+beta*(sum(Theta[v-w]*(w+1)(w+2)*Theta[w+2], w = 0 .. v))+beta*(sum(sum(delta[Zeta-1]*Theta[w-m]*(v-w+1)*(v-w+2)*Theta[v-w+2], Zeta = 0 .. w), w = 0 .. v)), Theta(v+2)) end do

Error, invalid left hand side in assignment

  
> 

``

Download Thermal_Stress_Distribution.mw

"How do I solve this''?...

Asked by: Saha 75
syntax user-error
December 18 2022
0  2

PROBLEM-IP.mw

What to do to get an upward curve in f(eta) and downward curve in theta(eta)?I get the curve wrong.Please Help.
> 

restart``
> 

with(PDEtools):
> 

ODE1 := diff(f(x), x, x, x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[f(x)])*(eta*cos(omega)+f(x)*`sinω`)*(diff(f(x), x, x))+(1-phi)^2.5*M*sin(alpha)^2*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[f(x)])*Gr[x]*`cosγ`*theta = 0;

diff(diff(diff(f(x), x), x), x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[f(x)])*(eta*cos(omega)+f(x)*`sinω`)*(diff(diff(f(x), x), x))+(1-phi)^2.5*M*sin(alpha)^2*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[f(x)])*Gr[x]*`cosγ`*theta = 0

 (1)
> 

ODE2 := K[nf]*(diff(theta(x), x, x))/K[f]+(1/2)*Pr*(eta*cos(omega)+`sinω`*f)*(diff(theta(x), x)) = 0;

K[nf]*(diff(diff(theta(x), x), x))/K[f]+(1/2)*Pr*(eta*cos(omega)+`sinω`*f)*(diff(theta(x), x)) = 0

 (2)
> 

Implementation*of*Differential*Transformation*Method*(Fe[3]*O[4]-Water);

Implementation*of*Differential*DifferentialGeometry:-Transformation*Method*(Fe[3]*O[4]-Water)

 (3)
> 

rho[s] := 5200:
> 

for k from 0 to 15 do F(k+3) := -((1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[f])*(sum((m+1)*(m+2)*F(m+2)*(cos(omega)*delta(k-m+1)+sin(omega)*F(k-m)), m = 0 .. k))+(1-phi)^2.5*M*sin(alpha)^2*(1-(k+1)*F(k+1))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[f])*Gr[x]*cos*gamma*Theta(k))/((k+1)*(k+2)*(k+3)) end do:
> 

for k from 0 to 15 do Theta(k+2) := -(1/2)*Pr*(sum((m+1)*Theta(m+1)*(cos(omega)*delta(k-m+1)+sin(omega)*F(k-m)), m = 0 .. k))/((k+1)*(k+2)*K[nf]) end do:
> 

f(eta) := simplify(sum(F(i)*eta^i, i = 0 .. 5), symbolic):
> 

indets(f(eta));

{A, B, M, alpha, cos, eta, omega, phi, Gr[x], (1.-phi)^(1/2), sin(alpha), sin(omega)}

 (4)
> 

Theta(eta) := simplify(sum(Theta(i)*eta^i, i = 0 .. 5), symbolic):
> 

indets(Theta(eta));

{A, B, M, Pr, alpha, cos, eta, omega, phi, Gr[x], (1.-phi)^(1/2), sin(alpha), sin(omega)}

 (5)
> 

> 

U := subs(M = 1, phi = .1, sin(alpha) = 1, cos(omega) = 1, sin(omega) = 0, cos*gamma = 1/sqrt(2), Gr[x] = 0, Pr = 6.2, `ρβ`[f] = 20939.1, K[nf] = .6842, f(eta));

(A-0.3201806132e-1*eta^2-0.1772797452e-1*eta^3-.1280722452*eta+0.6403612263e-1*eta^2*A)*eta^2

 (6)
> 

V := subs(Pr = 6.2, phi = .1, sin(alpha) = 1, M = 1, cos(omega) = 1, sin(omega) = 0, Gr[x] = 0, cos(gamma) = 1/sqrt(2), `ρβ`[f] = 20939.1, K[nf] = .6842, Theta(eta));

1.000000000+1.000000000*B*eta

 (7)
> 

e1 := subs(eta = 0, U) = 0;

0 = 0

  

diff*(0, 0) = 0

  

-14.09573093+12.91590125*A = 1

  

1.000000000 = 1

  

1.000000000+3.000000000*B = 0

  

{1.000000000 = 1, -14.09573093+12.91590125*A = 1, 1.000000000+3.000000000*B = 0}

  

{A = 1.168771009, B = -.3333333333}

 (8)
> 

U1 := subs(A = 1.168771009, B = -.3333333333, U);

(1.168771009+0.4282550234e-1*eta^2-0.1772797452e-1*eta^3-.1280722452*eta)*eta^2

 (9)
> 

V1 := subs(A = 1.168771009, B = -.3333333333, V);

1.000000000-.3333333333*eta

 (10)
> 

p1 := M = 1, phi = .1, A = 1.168771009*(1/2), B = -.3333333333, sin(alpha) = 1, cos(omega) = 1, sin(omega) = 0, cos*gamma = 1/sqrt(2), Gr[x] = 0, Pr = 6.2:
> 

s1 := subs(p1, f(eta)):
> 

d1 := diff(s1, `$`(eta, 1)):
> 

plot([d1, d2, d3], eta = 0 .. 1);

Warning, expecting only range variable eta in expression (.1080744102e-1*eta-.7183594026e-1*eta^2-.1280722452-.3441379740e-1*cos+.5735632896e-2*eta*cos-.3966707070e-2*eta^2*cos)*eta^2+2*(.5843855045+.540372051e-2*eta^2-.2394531342e-1*eta^3-.1280722452*eta-.3441379740e-1*eta*cos+.2867816448e-2*eta^2*cos-.1322235690e-2*eta^3*cos)*eta to be plotted but found name cos

  

 
> 

p1 := M = 1, phi = .1, A = 1.168771009*(1/2), B = -.3333333333, sin(alpha) = 1, cos(omega) = 1, sin(omega) = 0, cos*gamma = 1/sqrt(2), Gr[x] = 0, Pr = 6.2:
> 

s1 := subs(p1, Theta(eta)):
> 

plot([s1, s2, s3], eta = 0 .. 1);

Warning, expecting only range variable eta in expression 1.000000000-.3333333333*eta+.6611393405e-1*eta^4-.8693603699e-2*eta^5-.2336024612e-2*eta^5*cos to be plotted but found name cos

  

 
> 

``

Download PROBLEM-IP.mw

''How to plot a curve ?"...

Asked by: Saha 75
homework plotting duplicate-question
December 12 2022
0  1

Please Help! I am a beginner in using maple, I am getting  error in theta(k+2).How to get  m1,m2,m3 values from theta(k+2) and how to plot a curve theta(eta) like,f(eta). pp_DTM.mw

Download pp_DTM.mw

''How to rectify this error?''...

Asked by: Saha 75
plotting
December 09 2022
0  2

Please help! How to solve this?

pp-velocity.mw

How to resolve of this ODE System Error...

Asked by: mpopopopops 5
ode dsolve numeric
December 02 2022
1  1

How to resolve of this ODE System Error, the proble is "Error, (in dsolve/numeric/process_input) input system must be an ODE system, got independent variables {t, [S, E, In, R]}

Download TKO.mw

how can I calculate this Fourier transforms analy...

Asked by: kiraMou 45
integration fourier-transform
December 01 2022
1  6

FourTrans.mwI am wondering if I can calculate this Fourier Transfrom analytically using Maple:

Simplification of polynomials on radicals....

Asked by: ContrapuntoBrowniano 40
polynomial radical simplification
November 28 2022
0  1

Hi! I need to simplify a polynomial over an arbitrary field F on a radical variable ν, such that for some power j, ν is in F.

the polynomial takes the form:

k0+a1ν1/j+...+anνn/j

but, since νj is in F, this divides every element νi/j into a cotient group of order j, so this can be rewritten as a radical extension of F, in the form: 

k1+(a1+...)ν1/j+...+(ak+...)ν(j-1)/j

where ki, ai are in F. I feel like this is a very straightforward technique for handwritten algebra, but i can't see a command for this on Maple 18. Nor the simplify(..., radicals) Or the combine(...,radicals) seems to help here. Maybe there's an special command for this? Must i do it myself? Please, any help is aprecciated.

"unable to compute solution for t>HFloat(0.0):" Ho...

Asked by: 666 jvbasha 130
pde numeric pdsolve
November 17 2022
0  7

Hi maple users 

I am working with the PDE solver.

i am receiving a following error

"unable to compute solution for t>HFloat(0.0):" 

Kindly do the needful how to rectify this error.

dumm.mw

How can we determine the highest (or peak) value?...

Asked by: WA573 80
maximize plotting
November 07 2022
2  3

How can we determine the highest (or peak) value for every hump and at what value of (x, t) in the plot?

peak.mw

finding the dispersion relation...

Asked by: WA573 80
homework pde pdsolve
October 28 2022
1  2

How to find the dispersion relation (3) in Maple? I try it by equating the coefficients of cos(kx + wt) and its derivative to zero, but couldn't find it.

> 

restart
> 

with(LinearAlgebra):
> 

alias(E = E(x, t), P = P(x, t), N = N(x, t), Q = Q(x, t));

E, P, N, Q

 (1)
> 

d := Matrix([[c*(diff(E, t))+diff(E, x), alpha*P, 0, 0], [N[0]*E, diff(P, t), -sigma[2]*Q[0]*N/N[0], sigma[2]*omega[0]*Q], [0, sigma[1]*sigma[2]*Q[0]*P/N[0], diff(N, t), 0], [0, -omega[0]*P, 0, diff(Q, t)]]) = 0

d := (Matrix(4, 4, {(1, 1) = c*(diff(E(x, t), t))+diff(E(x, t), x), (1, 2) = alpha*P(x, t), (1, 3) = 0, (1, 4) = 0, (2, 1) = N(x, t)[0]*E(x, t), (2, 2) = diff(P(x, t), t), (2, 3) = -sigma[2]*Q(x, t)[0]*N(x, t)/N(x, t)[0], (2, 4) = sigma[2]*omega[0]*Q(x, t), (3, 1) = 0, (3, 2) = sigma[1]*sigma[2]*Q(x, t)[0]*P(x, t)/N(x, t)[0], (3, 3) = diff(N(x, t), t), (3, 4) = 0, (4, 1) = 0, (4, 2) = -omega[0]*P(x, t), (4, 3) = 0, (4, 4) = diff(Q(x, t), t)})) = 0

 (2)
> 

> 

re := -c*omega^3-kappa*omega^2+(c*omega[0]*(N[0]^2*sigma[2]+Q[0]^2*sigma[1])/N[0]^2-alpha*N[0])*omega+kappa*omega[0]*(N[0]^2*sigma[2]+Q[0]^2*sigma[1])/N[0]^2 = 0;

-c*omega^3-kappa*omega^2+(c*omega[0]*(Q[0]^2*sigma[1]+N[0]^2*sigma[2])/N[0]^2-alpha*N[0])*omega+kappa*omega[0]*(Q[0]^2*sigma[1]+N[0]^2*sigma[2])/N[0]^2 = 0

 (3)
> 

``

``

Download d_relation.mw

why can not evaluate the integration...

Asked by: bashar27 55
integration differential-equations incomplete-question
October 13 2022
0  1

restart;
alias(u = u(x, y, z, t), f = f(x, y, z, t));
                              u, f
u := (diff(f, y, x, x, x, x))*f^4-(diff(f, x, x, x, x))*(diff(f, y))*f^3-4*(diff(f, y, x, x, x))*(diff(f, x))*f^3+8*(diff(f, x, x, x))*(diff(f, x))*(diff(f, y))*f^2-4*(diff(f, x, x, x))*(diff(f, y, x))*f^3+12*(diff(f, y, x, x))*(diff(f, x))^2*f^2-36*(diff(f, x, x))*(diff(f, x))^2*(diff(f, y))*f+24*(diff(f, x, x))*(diff(f, x))*(diff(f, y, x))*f^2-6*(diff(f, x, x))*(diff(f, y, x, x))*f^3+6*(diff(f, x, x))^2*(diff(f, y))*f^2-24*(diff(f, x))^3*(diff(f, y, x))*f+24*(diff(f, x))^4*(diff(f, y))+2*(diff(f, y, x, t))*f^4-2*(diff(f, y, x))*(diff(f, t))*f^3-2*(diff(f, x, t))*(diff(f, y))*f^3+4*(diff(f, x))*(diff(f, y))*(diff(f, t))*f^2-2*(diff(f, x))*(diff(f, y, t))*f^3-3*(diff(f, z, x, x))*f^4+3*(diff(f, x, x))*(diff(f, z))*f^3+6*(diff(f, x))*(diff(f, z, x))*f^3-6*(diff(f, x))^2*(diff(f, z))*f^2;

5 6 7 8 9 10 11 Last Page 7 of 87