MaplePrimes Questions

Hi, 

Let's suppose 

f := piecewise(x < 1, x, -x)

Using 

discont(f, x)

says me x=1 is a discontinuity point.
But can I obtain that f is continuous over  these two intervals?

(-infinity, 1)
[1, +infinity)


Thanks in advance

PS: I'm sure I've seen this type of question here before, but I couldn't find them.


lign1 := S*`union`(T intersect R);
                 lign1 := S (`intersect`(R, T))

lign2 := (S union T) intersect (S union R);
       lign2 := `intersect`(`union`(R, S), `union`(S, T))

verify({lign1}, {lign2});

 

My textbook says the expression lign1 and lign2 should be the same, but it outputs false. I suspect it's because Maple is ignoring the parenthesis in the first expression.

 


I m Stucked with this step

restart; with(plots)

eq1 := 2*n*4^n*eta^((n+1)*(1/2))*(diff(f(eta), `$`(eta, 2)))^(n-1)*(diff(f(eta), `$`(eta, 3)))+4^n*(n+1)*eta^((n-1)*(1/2))*(diff(f(eta), `$`(eta, 2)))^n+4*f(eta)*(diff(f(eta), `$`(eta, 2)))-4*m*(diff(f(eta), eta))^2+m-2*M*(diff(f(eta), eta)) = 0;

2*n*4^n*eta^((1/2)*n+1/2)*(diff(diff(f(eta), eta), eta))^(n-1)*(diff(diff(diff(f(eta), eta), eta), eta))+4^n*(n+1)*eta^((1/2)*n-1/2)*(diff(diff(f(eta), eta), eta))^n+4*f(eta)*(diff(diff(f(eta), eta), eta))-4*m*(diff(f(eta), eta))^2+m-2*M*(diff(f(eta), eta)) = 0

(1)

eq2 := 2*eta*(diff(theta(eta), `$`(eta, 2)))+2*(diff(theta(eta), eta))+Pr*(f(eta)*(diff(theta(eta), eta))-s*(diff(f(eta), eta))*theta(eta))+Nb*(diff(theta(eta), eta))*(diff(phi(eta), eta))+Nt*(diff(theta(eta), `$`(eta, 2)))^2 = 0;

2*eta*(diff(diff(theta(eta), eta), eta))+2*(diff(theta(eta), eta))+Pr*(f(eta)*(diff(theta(eta), eta))-s*(diff(f(eta), eta))*theta(eta))+Nb*(diff(theta(eta), eta))*(diff(phi(eta), eta))+Nt*(diff(diff(theta(eta), eta), eta))^2 = 0

(2)

eq3 := 2*eta*(diff(phi(eta), `$`(eta, 2)))+2*(diff(phi(eta), eta))+Sc*(f(eta)*(diff(phi(eta), eta))-s*(diff(f(eta), eta))*phi(eta))+Nb*(2*eta*(diff(theta(eta), `$`(eta, 2)))+2*(diff(theta(eta), eta)))/Nt = 0;
``

2*eta*(diff(diff(phi(eta), eta), eta))+2*(diff(phi(eta), eta))+Sc*(f(eta)*(diff(phi(eta), eta))-s*(diff(f(eta), eta))*phi(eta))+Nb*(2*eta*(diff(diff(theta(eta), eta), eta))+2*(diff(theta(eta), eta)))/Nt = 0

(3)

eq4 := 2*eta*(diff(chi(eta), `$`(eta, 2)))+2*(diff(chi(eta), eta))+Lb*(f(eta)*(diff(chi(eta), eta))-s*(diff(f(eta), eta))*chi(eta))-Pe*(2*eta*chi(eta)*(diff(phi(eta), `$`(eta, 2)))+2*chi(eta)*(diff(phi(eta), eta))+2*eta*(diff(chi(eta), `$`(eta, 2)))*(diff(phi(eta), `$`(eta, 2)))) = 0;

2*eta*(diff(diff(chi(eta), eta), eta))+2*(diff(chi(eta), eta))+Lb*(f(eta)*(diff(chi(eta), eta))-s*(diff(f(eta), eta))*chi(eta))-Pe*(2*eta*chi(eta)*(diff(diff(phi(eta), eta), eta))+2*chi(eta)*(diff(phi(eta), eta))+2*eta*(diff(diff(chi(eta), eta), eta))*(diff(diff(phi(eta), eta), eta))) = 0

(4)

bcs := (D(f))(a) = 0, f(a) = 2*s*a*(D(phi))(a)/Sc, theta(a) = 1, phi(a) = 1, chi(a) = 1, (D(f))(10) = 1/2, theta(10) = 0, phi(10) = 0, chi(10) = 0;

(D(f))(a) = 0, f(a) = 2*s*a*(D(phi))(a)/Sc, theta(a) = 1, phi(a) = 1, chi(a) = 1, (D(f))(10) = 1/2, theta(10) = 0, phi(10) = 0, chi(10) = 0

(5)

params := {Lb = .1, M = .1, Nb = .6, Nt = .2, Pe = 5, Pr = 6.2, Sc = .1, a = 0.1e-1, m = 1/3, n = 1, s = .1};

{Lb = .1, M = .1, Nb = .6, Nt = .2, Pe = 5, Pr = 6.2, Sc = .1, a = 0.1e-1, m = 1/3, n = 1, s = .1}

(6)

sol := dsolve(eval([eq1, eq2, eq3, eq4, bcs], params), numeric, output = listprocedure, maxmesh = 1024)

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

 

odeplot(sol, [[eta, f(eta)][eta, phi(eta)], [eta, theta(eta)], [eta, chi(eta)]], eta = eval(c, params) .. 10)

Error, invalid subscript selector

 

NULL


 

Download Pblm2.mw

Hello. I want to solve a 2X2 system using Maple. I have written this code but I do not get expected results.

restart;
with(DEtools);
with(plottools);
with(plots);
with(Algebraic);

F := -1 - y - exp(x);
G := x^2 + y*(exp(x) - 1);

sol1 := eliminate({F, G}, {x, y});

sol1 := [{x = RootOf(-exp(_Z)^2 + _Z^2 + 1), y = -exp(RootOf(-exp(_Z)^2 + _Z^2 + 1)) - 1}, {}]

But, the solutions are x=0 and y=-2 (I have solved the system by hand).

 

Any help?!

Hi,how to separate the letters with commas, to put them in a set like { a,b,c,d} instead of {abcd} ?

Thanks

 CAP13CombinatoireQuestion.mw?

When I try to solve this problem above error displays

Here is the file
 

Download Main.mw

Hi,

When I do large computations in a worksheet and that I manually save it, I get a message in a pop-up window which asks me if I want to save these calculations too:

  1. What is it that is really saved if I answer "YES" 
  2. As I always answer "NO" how can I force MAPLE to ApplyGlobally this answer (I wasn't capable to fin it out in the Preferences>Interface menu).
  3. is this preference accounted for during AutoSaving?

Always in the goal to save only the necessary informations (in a mw format, export to mpl format could of course be a solution), how can I remove all the outpouts (and more specifically all the plots) in the current worksheet?
Presently I'm doing this and, before my last  saving, I turn TF to false and rexecute all _display commands ... not very practical, I agree

restart

_display := proc(tf, graph)
if tf then print(graph) end if:
end proc;

TF := true: # or false if no plot is required

# for each plot:
p := display( something ):
_display(TF, p)



Thanks in advance

How to represent this function in maple?

h(x) is sin(x)/x

I use the product() but do not how to represent the step here like 1,3,5,7.... And also the integral term not working. Making Maple frozen

 

Thank you!

Dear all

I have two matrices represent ( f(x[i],y[j])  and  g(x[i],y[j]) of a  given two dimensional functions.

How can I use thoses values of matrices to compute the largest distance between the two curves.

Many thanks

 

 

 


When i try to solve the problem i got error like this 

Please help me get the solution

restart

with(plots):

eq1 := (2*(eta*(diff(f(eta), `$`(eta, 3)))+diff(f(eta), `$`(eta, 2))))/((1-phi[1])^2.5*(1-phi[2])^2.5)+((1-phi[2])*(1-phi[1]+phi[1]*rho[s1]/rho[f])+phi[2]*rho[s2]/rho[f])*f(eta)*(diff(f(eta), `$`(eta, 2)))-sigma[hnf]*M*(diff(f(eta), eta))/(sigma[f]*(m^2+1)) = 0;

2*(eta*(diff(diff(diff(f(eta), eta), eta), eta))+diff(diff(f(eta), eta), eta))/((1-phi[1])^2.5*(1-phi[2])^2.5)+((1-phi[2])*(1-phi[1]+phi[1]*rho[s1]/rho[f])+phi[2]*rho[s2]/rho[f])*f(eta)*(diff(diff(f(eta), eta), eta))-sigma[hnf]*M*(diff(f(eta), eta))/(sigma[f]*(m^2+1)) = 0

(1)

eq2 := k[hnf]*(eta*(diff(theta(eta), `$`(eta, 2)))+diff(theta(eta), eta))/k[f]+(1/2)*Pr*((1-phi[2])*(1-phi[1]+phi[1]*`&rho;C`[p][s1]/`&rho;C`[p][f])+phi[2]*`&rho;C`[p][s2]/`&rho;C`[p][f])*f(eta)*(diff(theta(eta), eta)) = 0;

k[hnf]*(eta*(diff(diff(theta(eta), eta), eta))+diff(theta(eta), eta))/k[f]+(1/2)*Pr*((1-phi[2])*(1-phi[1]+phi[1]*`&rho;C`[p][s1]/`&rho;C`[p][f])+phi[2]*`&rho;C`[p][s2]/`&rho;C`[p][f])*f(eta)*(diff(theta(eta), eta)) = 0

(2)

eq3 := 2*(1-phi[1])^2.5*(1-phi[2])^2.5*(eta*(diff(chi(eta), `$`(eta, 2)))+diff(chi(eta), eta))+Sc*f(eta)*(diff(chi(eta), eta))-beta*Sc*(1+delta*theta(eta))^n*exp^(-E/(1+delta*theta(eta)))*chi(eta) = 0;

2*(1-phi[1])^2.5*(1-phi[2])^2.5*(eta*(diff(diff(chi(eta), eta), eta))+diff(chi(eta), eta))+Sc*f(eta)*(diff(chi(eta), eta))-beta*Sc*(1+delta*theta(eta))^n*exp^(-E/(1+delta*theta(eta)))*chi(eta) = 0

(3)

bcs := f(c) = (1/2)*lambda*c, (D(f))(c) = (1/2)*lambda, theta(c) = 1, chi(c) = 1, (D(f))(10) = (1-lambda)*(1/2), theta(10) = 0, chi(10) = 0;

f(c) = (1/2)*lambda*c, (D(f))(c) = (1/2)*lambda, theta(c) = 1, chi(c) = 1, (D(f))(10) = 1/2-(1/2)*lambda, theta(10) = 0, chi(10) = 0

(4)

params := {E = .5, M = .1, Pr = 6.2, Sc = .5, beta = .5, c = .1, delta = .5, lambda = -1, m = .5, n = .1, k[f] = .6130, k[hnf] = 1.039626683, phi[1] = .1, phi[2] = .1, rho[f] = 997.1, rho[s1] = 5180, rho[s2] = 1800, sigma[f] = 0.5e-1, sigma[hnf] = 0.8888888337e-1, `&rho;C`[p][f] = 4.1668809*10^6, `&rho;C`[p][s1] = 3470600, `&rho;C`[p][s2] = 1290600};

{E = .5, M = .1, Pr = 6.2, Sc = .5, beta = .5, c = .1, delta = .5, lambda = -1, m = .5, n = .1, k[f] = .6130, k[hnf] = 1.039626683, phi[1] = .1, phi[2] = .1, rho[f] = 997.1, rho[s1] = 5180, rho[s2] = 1800, sigma[f] = 0.5e-1, sigma[hnf] = 0.8888888337e-1, `&rho;C`[p][f] = 4166880.900, `&rho;C`[p][s1] = 3470600, `&rho;C`[p][s2] = 1290600}

(5)

ODEs := [eq1, eq2, eq3]

[2*(eta*(diff(diff(diff(f(eta), eta), eta), eta))+diff(diff(f(eta), eta), eta))/((1-phi[1])^2.5*(1-phi[2])^2.5)+((1-phi[2])*(1-phi[1]+phi[1]*rho[s1]/rho[f])+phi[2]*rho[s2]/rho[f])*f(eta)*(diff(diff(f(eta), eta), eta))-sigma[hnf]*M*(diff(f(eta), eta))/(sigma[f]*(m^2+1)) = 0, k[hnf]*(eta*(diff(diff(theta(eta), eta), eta))+diff(theta(eta), eta))/k[f]+(1/2)*Pr*((1-phi[2])*(1-phi[1]+phi[1]*`&rho;C`[p][s1]/`&rho;C`[p][f])+phi[2]*`&rho;C`[p][s2]/`&rho;C`[p][f])*f(eta)*(diff(theta(eta), eta)) = 0, 2*(1-phi[1])^2.5*(1-phi[2])^2.5*(eta*(diff(diff(chi(eta), eta), eta))+diff(chi(eta), eta))+Sc*f(eta)*(diff(chi(eta), eta))-beta*Sc*(1+delta*theta(eta))^n*exp^(-E/(1+delta*theta(eta)))*chi(eta) = 0]

(6)

``

sol := dsolve(eval([ODES, bcs], params), numeric, output = listprocedure)

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

 

``


 

Download 1.mw
 

 

 

 

I am doing the first interactive plot in Tutorial 9:  Dynamic Applications.    The first part is Exploration Assistant.   This is mostly going OK but I don't want to have it rescale the y-axis to the data extents.   So I uncheck the box and type in my min & max values as shown below.    But when I drag the slider for coefficient a or b, the plot reverts back to using data extents.   How do I override this so that the axis range stays fixed when I drag the sliders?

I'd like to generate a string where the first part is changeable, and the second one needs to include apostrophs.

NODESteelProfiles_I:-Property("allmembers")

The section until til left bracket can be stored in one variable, but how can I add "allmembers") into the string?

Hi,

The evalDG command included in the LieAlgebras via the DifferentialGeometry package allows summations such as evalDG(e1+e2) where e1 and e2 are two generators of a given Lie algebra.

Now, let L be a list of such summations, e.g.,

L= [e1+a*e3,2*e2+b*e4,e1+3*e5]

where a and b are symbols for variable names in the real domain.
Then, while the code

evalDG(add(eval(cat(e,i)),i=1..n));

works fine (n be a given integer), the code

evalDG(add(eval(L[i]),i=1..n));

does not, and an error message results due to an "invalid subscript selector".

What is the right code to realize this summation?

Thanks.

Jaqr

Dear all,
I would like to get the Null Vectors of some matrix B but not those which are returned by NullSpace(B), i would like to get specific structure of the Kernel, i mean, i want to get the following structure of the Kernel when the matrix B is 10*10 for example:

Kernel= {
Vector[column](10, [0, 0, 0, 0, 0, 0, 0, 0, 0, alpha]),
Vector[column](10, [0, 0, 0, 0, 0, 0, 0, 0, alpha, beta]),
Vector[column](10, [0, 0, 0, 0, 0, 0, 0, alpha,beta, gamma]),
Vector[column](10, [0, 0, 0, 0, 0, 0, alpha, beta, gamma, delta])
              }
Is it possible to do that with Maple ? Have you some ideas ? 
In other words,  could you help me please ?

Best regards

What do folks think about adding \, space automatically after Maple's Latex output for sqrt command?

This is all using the new Latex command. Not the orginal latex() command.

For me, it looks much better. I did some tests, and see no negative effect by always adding \, at the end of \sqrt{}

The reason is this: If there is a symbol after sqrt(.....)<HERE> then the math looks like the symbol is almost inside the integral, when it should be out. By adding \, this makes it more clear and better looking.

It is possible to write macro in Latex to do this in preamble I suppose. But if it can be done in the source, then that will be much better as it will always be there without the need to write complicated macro to change it later in Latex. 

Here are 4 expressions to compare

One can see the current problem area and the effect of adding \, automatically. The space was added for each case above in the second column, even for those where it was needed (second and third expressions) to show that adding space has no negative side-effect even when not needed.

Here is the above raw latex used, compiled on Tex Live 2020 on Linux,  and the Maple worksheet below that.

\documentclass[12pt]{article}
\usepackage{amsmath}
\usepackage{longtable}
\usepackage{mleftright}  
\mleftright

\begin{document}

\begin{longtable}{|p{3in}|p{3in}|}\hline
Default. No space after & With space added after\\\hline 
\[
\sqrt{\left(3+x \right)^{2}+1} a +\sin \left(x \right)
\]
&
\[
\sqrt{\left(3+x \right)^{2}+1}\, a +\sin \left(x \right)
\]\\\hline 
\[
\frac{\sqrt{\left(3+x \right)^{2}+1}}{5}+b
\]
&
\[
\frac{\sqrt{\left(3+x \right)^{2}+1}\,}{5}+b
\]\\\hline 
\[
\frac{\sqrt{x}}{5}+b
\]
&
\[
\frac{\sqrt{x}\,}{5}+b
\]
\\\hline 
\[
\frac{\sqrt{x^{2}+6} c}{5}+b
\]
&
\[
\frac{\sqrt{x^{2}+6}\, c}{5}+b
\]\\\hline 
\end{longtable}

\end{document}


 

interface(version);
Physics:-Version();

`Standard Worksheet Interface, Maple 2020.2, Windows 10, November 11 2020 Build ID 1502365`

`The "Physics Updates" version in the MapleCloud is 882. The version installed in this computer is 881 created 2020, November 21, 11:11 hours Pacific Time, found in the directory C:\Users\me\maple\toolbox\2020\Physics Updates\lib\`

restart;

Latex:-Settings(UseImaginaryUnit=i,
      UseColor = false,
      powersoftrigonometricfunctions= computernotation,
      leavespaceafterfunctionname = true,
      cacheresults = false
);

`* Partial match of  '`*UseImaginaryUnit*`' against keyword '`*useimaginaryunit*`' `

`* Partial match of  '`*UseColor*`' against keyword '`*usecolor*`' `

[powersoftrigonometricfunctions = computernotation, leavespaceafterfunctionname = true, cacheresults = false, useimaginaryunit = i, usecolor = false]

expr1:=sqrt((3+x)^2+1)*a+sin(x);
Latex(expr1)

((3+x)^2+1)^(1/2)*a+sin(x)

\sqrt{\left(3+x \right)^{2}+1} a +\sin \left(x \right)

expr2:=sqrt((3+x)^2+1)/5+b;
Latex(expr2)

(1/5)*((3+x)^2+1)^(1/2)+b

\frac{\sqrt{\left(3+x \right)^{2}+1}}{5}+b

expr3:=sqrt(x)/5+b;
Latex(expr3)

(1/5)*x^(1/2)+b

\frac{\sqrt{x}}{5}+b

expr4:=(sqrt((3+x^2)+3)*c)/5+b;
Latex(expr4)

(1/5)*(x^2+6)^(1/2)*c+b

\frac{\sqrt{x^{2}+6} c}{5}+b

 


(typo, should be sqrt, not int)

Download space_after_int.mw

First 446 447 448 449 450 451 452 Last Page 448 of 2402