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

How to do this?...

Asked by: Saha 80
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 80
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 80
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 80
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 80
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;

Extracting matrix of 208x2...

Asked by: ijuptilk 100
matrix indexing plotting
October 13 2022
0  2

Please how  can I use getdata to extract only 200x2 of 208x2 matrix from maple to excel.

For instance I have:

Q:= ( seq( seq( plottools:-getdata(ans1[s1,3])[j,3],j=1..3), s1=1));

It returns  208x2,  200x2 200x2 matrices and I can't extract it using 

`<|>`( seq( seq( plottools:-getdata(ans1[s1,3])[j,3],j=1..3), s1=1));

because the matrices have different dimension. Please, how can make 208x2 matrix to 200x2?

How to delete 0*I-2 from a list...

Asked by: ijuptilk 100
complex simplification floating-point
October 11 2022
1  2

Hi, please how do I delete only the 0*I of the complex expression: 0.*I - 3?

I tried this:  remove(`=`, [0.*I-3, 5, 7, -2], 0.*I);

But couldn't work.

How to linearize equation?...

Asked by: Ahmed111 130
linearize
September 22 2022
2  14

How to linearize eq (6) by neglecting all higher order terms i.e., epsilon[1]^2, epsilon[2]^2, epsilon[1]*epsilon[2]... etc? How to do it in maple?

Calculation.mw

First 6 7 8 9 10 11 12 Last Page 8 of 88