Hi,

I need to find out how to print the number of row or the numbers on right based only on the combination of the two numbers in first two columns. The external file looks like this:

I have prepared possible combinations of numbers and I need to read the numbers on the right which belongs to this combination. There is no such thing as the sum of two numbers = number of row, so it can't be done this easy.

Can somebody help please? Thanks a lot

Hello,

I would like to understand how i can construct module with submodules.

In the maple help, this chapter should answer perfectly to my need : https://www.maplesoft.com/support/help/Maple/view.aspx?path=ProgrammingGuide/Chapter11

At this page, there is a package called Shapes which should be very useful for me.

However, i don't understand 1 point at the moment about the package architecture of this example of package.

Why there is a submodule also called Shapes inside the module Shapes ? In other words, why the different submodules point, segment, circle, square, triangle have not be constructed directly under the module Shapes but under the submodule Shapes?

I thank you in advance for your help.

Hi everyone, i'm working on a project and i'm having a problem as such:

I got t1 already (t1=26,88) by solving s(t)=0 and when I try to calculate s(t) with t=t1, i type s(t1) Maple give me 1.10^(-6) as result but the result should be 0.22144. Here are my codes:

hi, is there any option to make inequality in animate with filled/marked area?

restart;
with(plots);
display(animate(implicitplot, [abs(x + y*I - a) <= 2*(1 - abs(x + y*I)), x = -2 .. 2, y = -2 .. 2], a = -1.75 .. 1.75, frames = 100), implicitplot(x^2 + y^2 = 1, x = -2 .. 2, y = -2 .. 2, color = black));

Hey everyone,

I see I can use SimpleLieAlgebraData to create Lie algebras of types A, B, C, D, and also G2 and F4. Is there a built-in way to generate E6, E7, and E8? If not, is there any plan to add these?

Thanks!

These two animations do not work under Maple 2020. How do I make them work? What am I doing wrong?

plots[animate]( plot3d, [sin(A)*(x^2+y^2), x=-2..2, y=-2..2], A=0..2*Pi );

plots[animate]( plot, [A*sin(x), x=0..10], A=0..2 );

Thank you!

mapleatha

system_of_PDE.pdf

How to solve the system of partial differential equations in Maple.  I have attached a pdf file, please check it. Kindly help me.

Thank you

[Edit: I've also atached the OP's worksheet below.--Carl Love]

not_solve_sys_pde.mw

Hi everyone!

Is it possible to evaluate int(diff(u(x,y,t),x)*v(x,y,t)+diff(v(x,y,t),x)*u(x,y,t),x) directly as u(x,y,t)*v(x,y,t) in Maple?

Thanks.

hi, do you know how to use discont in this plot? (is it even possible?)

restart;
with(plots);
animate(implicitplot, [y^2 = x^3/(2*a - x), x = -5 .. 5, y = -5 .. 5], a = -5 .. 5);
 

  restart:
  ODES := (1+E)*(diff(f(eta), eta, eta, eta))+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2-E*delta*(diff(f(eta), eta, eta))^2*(diff(f(eta), eta, eta, eta))+M*(diff(f(eta), eta)-1)+Upsilon*(1-(diff(f(eta), eta))^2) = 0,
          diff(theta(eta), eta, eta)+(3/4)*R*(diff((1+(B-1)*theta(eta))^3*(diff(theta(eta), eta)), eta))+Pr*f(eta)*(diff(theta(eta), eta))+Pr*Nb*(diff(theta(eta), eta))*(diff(phi(eta), eta))+Pr*Nt*(diff(theta(eta), eta))*(diff(theta(eta), eta))+Pr*Df*(diff(phi(eta), eta, eta))+Pr*Ec*(diff(f(eta), eta, eta))^2+Q*theta(eta) = 0,
           diff(phi(eta), eta, eta)+Sc*(diff(phi(eta), eta))+Nt*(diff(theta(eta), eta, eta))/Nb-Sc*kappa*phi(eta)+Sc*Sr*(diff(theta(eta), eta, eta)) = 0:
  bcs:= f(0) = 0,
        D(f)(0) = 1,
        D(f)(10) = 0,
       
        
        (D(theta))(0) = -.2*(1-theta(0)),
        theta(10) = 0,
        phi(0) = 1,
        phi(10) = 0:
  params:=[E=EV,Nb = 3.9, lambda = .01, Sc = 1, Nt = .5, delta = 0.9, Pr = 7., M = mv, R =5., B = 6., Df = 1., Ec  =5., kappa = 1., Sr = 3., Upsilon =- 3., Q = 3.]:
  EVals:=[-1, -3, -7]:
  Mvals:= [3, 5, 7]:
  ans:=Matrix( numelems(Mvals)*numelems(EVals)+1, 5):
  ans[1,..]:= < 'M' | 'E' | expr1 | expr2 |expr3>:
  for k from 1 by 1 to 3 do
      mv:= Mvals[k]:
      for j from 1 by 1 to 3 do
          EV:=EVals[j]:
          sol:=dsolve( eval([ODES, bcs], params), numeric, output=listprocedure);
          ans[3*(k-1)+j+1,..]:= `<|>`([ mv,
                                        EV,
                                        R__e^(-0.5)*sh= eval( -diff(phi(eta), eta),
                                                              [sol[], params[]]
                                                            )(0),
                                        R__e^(-0.5)*NU= eval( -diff(theta(eta), eta),
                                                              [sol[], params[]]
                                                            )(0),
                                        R__e^(0.5)*C__f=eval( (1+0.5)*(diff(f(eta), eta,eta)-1/3*diff(f(eta), eta)^3),[sol[], params[]])(0)]);od:od:ans;
Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging
               Matrix(%id = 18446744074195107958)
 for k from 1 by 1 to 3 do
      plot( [ seq( [ seq( [ans[j,1], rhs(ans[j,2+k]) ], j=i..10,3 ) ], i=2..4 ) ],
            color=[red, green, blue],
            labels=[typeset(M), typeset( lhs(ans[2,2+k]) )],
            labelfont=[times, bold, 20],
            legend=[typeset(E=0.01),typeset(E=0.1),typeset(E=0.2)],
            legendstyle=[font=[times, bold, 20]],
            title=typeset( ans[1,2+k], " versus ", M, " parameterized by ", E),
            titlefont=[times, bold, 24]
          )
 od;
Error, invalid input: rhs received 0, which is not valid for its 1st argument, expr

I trying to display factorials in a matrix as e.g 6!4!/5! not 144.
 

Download How_to_display_factorial_unevaluated.mw

Dear Users!

For a given piecewise function S(x)

I want to generate the following functions phi[1], phi[2],... for n = 9

I shall be very thankful for your positive response. Please take care and thanks in advance.

I am not familiar with mathematica syntax. I have this code which is written .np and I would like to translate it efficiently to Maple syntax.I tried FromMmaNotebook but it is not helpful for me. I would be grateful if you would give me a little help me translate the code. I attached the code as diagon.txt (you can change .txt with .nb)

Thank you in advance!diagon.txt

Hi everyone

I wonder if the following elliptic PDEs system can be caculated by Maple 2018.

The PDEs system is

where

Maple sketched solutions of the above system, but unfortunately main portions (sine portion) of all the solutions seemed to overlap together.

Is it possible to output all the solutions filling the target region (or at least not overlapping) by revising the code?

I am not the sure if the overlap is caused by the boundary conditions(BC) or caused by the elliptic PDEs system??

Particularly, I do not understand how (D[2](x))(xi, 0) = 0 and (D[2](y))(xi, 0) = 0  work in the BCs. Is there any paper or textbook explaining this setup??

Thanks in advance.

The code is as follows (or the code can be downloaded from attatched mw QUESTION_20210423_001.mw):

restart;
with(plots);


alpha04 := (diff(x(xi, eta), eta))^2+(diff(y(xi, eta), eta))^2;

beta04 := (diff(x(xi, eta), xi))*(diff(x(xi, eta), eta))+(diff(y(xi, eta), xi))*(diff(y(xi, eta), eta));

gamma04 := (diff(x(xi, eta), xi))^2+(diff(y(xi, eta), xi))^2;

PDE04 := {alpha04*(diff(x(xi, eta), xi, xi))-2*beta04*(diff(x(xi, eta), xi, eta))+gamma04*(diff(x(xi, eta), eta, eta)) = 0, alpha04*(diff(y(xi, eta), xi, xi))-2*beta04*(diff(y(xi, eta), xi, eta))+gamma04*(diff(y(xi, eta), eta, eta)) = 0};


IBC04 := {x(0, eta) = eta, x(1, eta) = eta, x(xi, 0) = sin((2.0*Pi)*xi), y(0, eta) = eta, y(1, eta) = eta, y(xi, 0) = 100*xi, (D[2](x))(xi, 0) = 0, (D[2](y))(xi, 0) = 0};

pds04 := pdsolve(PDE04, IBC04, numeric, spacestep = 1/50);


xCurve04 := plots[display]([seq(pds04:-plot(x, eta = i), i = 0 .. 10)]);

yCurve04 := plots[display]([seq(pds04:-plot(y, eta = (1/10)*i), i = 0 .. 10)]);

How can maple be trusted app when window 10 device guard enabled?