I am new to Maple. I have been plotting 3-dimensional graphs, and would like to plot points on the graph, but I don't see how to do that.

Thanks.

Hi all

for a given two parameter in the attached code ( G and PA) to compute the coefficient Cr

I would like to display the text to describe the corresponding situation ( maybe we one add a table )

Many tahnks for your help

testprint.mw

I want to write a code for 3D animation of solutions of PDEs.  Specially, I wonder 3D animation of wave equations etc.

Example: I can write a code for 2d animation but not write 3d animation.

a:=-2:
b:=2:
alpha:=0.1:
beta:=1:
pde:=diff(uu(x,t),t)=alpha*diff(uu(x,t),x,x)+beta*uu(x,t)*(1-uu(x,t));
IC:=uu(x,0)=(sech(10*x))^2;
BC:=uu(a,t)=0,uu(b,t)=0;  

s:=pdsolve(pde,{IC,BC},numeric ,timestep=1/100, spacestep=1/100);
s :- animate( uu(x,t) ,t=0..tbas, frames=40, labels=["x", "u(x,t)"], labelfont=[TIMES,ROMAN,14]); 
g3:=s :- plot3d(uu(x,t), x=-xbas..xbas,t=0..tbas, shading=zhue, axes=boxed, labels=["x","t","u(x,t)"], labelfont=[TIMES,ROMAN,16]);

 But how to write a code for 3d animation?

Thanks.

plot([v[Iterations-4], v[Iterations-3], v[Iterations-2], v[Iterations-1], v[Iterations-1]], t = 1 .. 2.9, color = [red, green, blue, black, purple]); g11 := plot([v[Iterations-4], v[Iterations-3], v[Iterations-2], v[Iterations-1], v[Iterations-1]], t = 1 .. 2.9, color = [red, green, blue, black, purple], numpoints = 200);

with(plots); p11 := display({g11}, axes = boxed, thickness = 2); dat1 := `~`[plottools:-getdata]([p11]); for i while i <= 5 do A[i] := dat1[i, 3]; Y[i] := A[i][1 .. -1, 2] end do; X := A[1][1 .. -1, 1]; MM1 := `<|>`(X, Y[j2]$j2 = 1 .. 5);

        [[                                                        
        [[                                                        
dat1 := [["curve", [1. .. 2.9, 1.25843241253069330813854745573 .. 
        [[                                                        

                                    [ 200 x 2 Matrix       ]]  
                                    [ Data Type: sfloat    ]]  
  3.64143345197220537035927195897], [ Storage: rectangular ]], 
                                    [ Order: Fortran_order ]]  

  [                                                              
  [                                                              
  ["curve", [1. .. 2.9, 1.26033290400709038909782882918 .. 6.355\
  [                                                              

                               [ 200 x 2 Matrix       ]]  [
                               [ Data Type: sfloat    ]]  [
  77353302842195250192992449], [ Storage: rectangular ]], [
                               [ Order: Fortran_order ]]  [

  "curve", [1. .. 2.9, 1.26039280181499972068892216922 .. 9.3496\

                              [ 200 x 2 Matrix       ]]  [
                              [ Data Type: sfloat    ]]  [
  4246576526181629314962273], [ Storage: rectangular ]], [
                              [ Order: Fortran_order ]]  [

  "curve", [1. .. 2.9, 1.26039397590711135163918382958 .. 10.112\

                              [ 200 x 2 Matrix       ]]  [
                              [ Data Type: sfloat    ]]  [
  7233044121210212769952674], [ Storage: rectangular ]], [
                              [ Order: Fortran_order ]]  [

  "curve", [1. .. 2.9, 1.26039397590711135163918382958 .. 10.112\

                              [ 200 x 2 Matrix       ]]]
                              [ Data Type: sfloat    ]]]
  7233044121210212769952674], [ Storage: rectangular ]]]
                              [ Order: Fortran_order ]]]
Error, 1st Matrix dimension size (1..1) is too small to fit source size (1..200)
Error, number of indices exceeds rank
ExportMatrix("C:/Users/DELL/Dropbox/Common Data (UH)/Hamid made model/For Hamid Thesis/P4 Picard for ODEs/Tecplot//P5 case 1.dat", MM1)
Error, invalid input: ExportMatrix expects its 2nd argument, M, to be of type {Matrix, list(Matrix)}, but received Vector[row](6, {(1) = X, (2) = Y[1], (3) = Y[2], (4) = Y[3], (5) = Y[4], (6) = Y[5]})
 

Hi, in the 'assign' command below, a certain row from matrix 'R' is assigned to the variables RA, MA .. yB. The command works well. However, I wonder if it could not be better/ shorter. So, without having to write out R [1], R [2] .. R [8].

I would be very happy if someone could advise me on this.

Here is the algorithm:

restart:with(linalg):
R:=Matrix(2,8,[1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18]);

aFshear:=proc(nr,R)
assign(('RA','MA','thetaA','yA','RB','MB','thetaB','yB')
= (R[1],R[2],R[3],R[4],R[5],R[6],R[7],R[8]));
end:

aFshear(2,row(R,2));
print(thetaB);

I am solving a wave propagation problem and after lots of derivations I ended up with the following differential equations. Any idea how to decouple this system so I can get two decoupled equations of the following forms:

equation (1): = function of (Q)      Equation (2): = function of (w)

System of Equations:

When writing a post on MaplePrimes it is possible to access the source code, which allows to introduce/edit many attributes. Is it possible to observe the source of another post?

I want to evaluate a function of two variables at different values for these variables and collect the outputs in a matrix form. This is my best attempt at figuring it out -- using the string tools to evaluate a matrix at multiple values and spit out data. How can give each output labeled uniquely (not each line labeled A, but A,B,C... for n number of lines)? And How can I take the inner brackets out of the sequence output, to make each row of output data an array? Then I could refer to each data point individually with A[1],A[2],...,B[1],B[2],..., etc.

Please let me know if there is a better way to do this! Thank you, Jeremy

Also here is my code copied in text format for reference:

restart;
Zres:=5;
data:=[];
M:=Matrix(eval~(z*x,z=~[seq(n,n=1..Zres)]));
Array(eval~(z*x,z=~[seq(n,n=1..Zres)]));
with(LinearAlgebra):
for x from 1 to 5 by 0.5 do A:=[((op(data),seq(Column(M,n),n=1..Zres)))] od;

I am able to generate random polynomials with non-zero coefficents, and define sets of all the positive divisors of the leading coefficient and the constant terms. My question is this, how may I apply the rational zeros theorem to generate the set of all possible rational zeros of the polynomial. I basically need to form all the possible quotients (positive and negative) with numerator in one set and denominator in the other set, ignoring duplicates. The attached worksheet has what I've done so far.rational_zeros.mw
 

Download rational_zeros.mw

 Dear ...

Let us consider $F$ be a parametric polynomial set. How can decide that $F$ contains a constant or parameters, automatically? For example, if $F=/{a*x^3-1, b*y+x, -1, c/}$ where $a, b$ and $c$ are parameters and $x$ and $y$ are variables then the output of procedure is true since of $F$ contains $c$ and also $-1$ and the outputs for $G=\{\a*x^3-1, b*y+x, x+y}$ is false. Is there any command or simple procedure for this conclusion?

I am looking forward to hearing from you

Sincerely yours.

The command eval allows to simplify a complicated expression

in a more compact form for a later output in LATEX

In the example which follows I was able to insert A and B but non C in the expression.

There is already a post on this kind of topic but I failed to understand the details.

Perhaps a Maple worksheet of answer on this topic  would be useful !

bye Lorenzo
 

Download maple_primes_eval.mw

maple does not work at all

it displays this error

Error, (in StringTools:-FormatMessage) unknown option MAPLE
 

"From the list of positive integers 1; 2; 3; 4; : : :, remove every second number, leaving
1; 3; 5; 7; 9; : : : Since 3 is the flrst surviving number above 2 that has not been used
as the \killer", we remove every 3rd number from the remaining numbers, yielding
1; 3; 7; 9; 13; 15; 19; 21; : : :"

 I have tried using this:

A:= proc(n::integer)::list;
   local B,C,j;
   B:={seq(i,i=1..n)};
   for j from 2 by 2 to n do
       C:= {B minus {B[j]}};
   end do;
 end proc;

But I'm not getting the desired result.

Hi, I am trying to get back into the habit of making maplets for the sake of uncomplicating things, I never had trouble with loading this package before, but strangely now when i enter:

with(Maplets[Elements]);
maplet1 := Maplet(["Hello World!", Button("OK", Shutdown())]);
 

I get the error message:

Error, `Elements` does not evaluate to a module