Hi all,

I have the attached code, it returns a "unable to parse" error after the line with the "FOC__1D1" assignment. I can't find out why.

Thank you for your help in advance,

JTamas

1_2_test.mw

I'm am pretty new to maple, coming from a mathcad background, so sorry in advance if this is a dumb question.  My original need was to plot text objects with sold backgrounds on top of other plotted objects.  I saw in other posts that Maple doesn't natively support this.  So instead I'm trying to create a composite plot of objects by plotting text objects over polygons or rectangles.  However I can't seem to make a given plotted object "cover" another plotted object.  

Below is a simple example.  The easy analogy is just that I want to plot these objects in “layer order”, with L1 being the top layer.  So I would like the polygon to opaquely obscure the “underlying” contour plot, and then in turn, the text object to behave as a “top-most layer” with the polygon acting as a background for the text. 

L1 := textplot([2, 2, "Polygon"], color = white);

L2 := polygon([[0, 0], [3, 4], [3, 1]], color = red);

L3 := contourplot(x^2 + y^2, x = 1 .. 2, y = 1 .. 2);

display(L3, L2, L1); 

Why does "table" appear right under my equation? Like

Eq1 := .9*(diff(f(eta), eta, eta, eta))+.1*We*(diff(f(eta), eta, eta))*(diff(f(eta), eta, eta, eta))-(diff(f(eta), eta))^2+f(eta)*(diff(f(eta), eta, eta))^2+0.6e-1-.5*(diff(f(eta), eta))+.1*(table([w = 1.1]))(eta)*(1+beta1*(table([w = 1.1]))(eta))+.1*phi(eta)*N*(1+.1[c]*phi(eta)) = 0

Hi, 

I am trying to solve a system of 3 ODE of order 4. 

Is there any template to do it? 

When I try it, I get a message:

Error, (in dsolve/numeric/bvp/convertsys) too few boundary conditions: expected 10, got 8
What is wrong with it?

> restart;
> with*plots; blt := 5;
> lambda := .1; delta := .1;
                                     0.1
                                     0.1
> gamma2 := .1;
                                     0.1
> lambda2 := .1; lambda[T] := .7; lambda[C] = .1;
> Pr := 2;
> Nb := 0.1e-1;
> Nt := 0.1e-1;
> Sc := 0.1e-1;
> M := .5; E = -.1;
> N := 0.1e-1; n = .1;
> sigma[1] := .1;
> alpha[0] := 0; omega := 0;
                                      0
                                      0
> Eq1 := diff(f(eta), eta, eta, eta)+gamma2*((diff(f(eta), eta, eta))*(diff(f(eta), eta, eta))-f(eta)*(diff(f(eta), eta, eta, eta, eta)))-(1+lambda2)*((diff(f(eta), eta))*(diff(f(eta), eta))-f(eta)*(diff(f(eta), eta, eta)))-(1+lambda2)*M*sin(omega)^2*(diff(f(eta), eta))+(1+lambda2)(lambda*(1+lambda[T]*theta(eta))*theta(eta)+lambda*N*(1+lambda[C]*phi(eta))*phi(eta))*cos(alpha[0]) = 0;
                                                           2
/  d   /  d   /  d         \\\       /  d   /  d         \\ 
|----- |----- |----- f(eta)||| + 0.1 |----- |----- f(eta)|| 
\ deta \ deta \ deta       ///       \ deta \ deta       // 

                                                                           2
                /  d   /  d   /  d   /  d         \\\\       /  d         \ 
   - 0.1 f(eta) |----- |----- |----- |----- f(eta)|||| - 1.1 |----- f(eta)| 
                \ deta \ deta \ deta \ deta       ////       \ deta       / 

                /  d   /  d         \\          
   + 1.1 f(eta) |----- |----- f(eta)|| + 1.1 = 0
                \ deta \ deta       //          
> Eq2 := 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))) = 0;
/  d   /  d             \\            /  d             \
|----- |----- theta(eta)|| + 2 f(eta) |----- theta(eta)|
\ deta \ deta           //            \ deta           /

                                                                       2    
          /  d             \ /  d           \        /  d             \     
   + 0.02 |----- theta(eta)| |----- phi(eta)| + 0.02 |----- theta(eta)|  = 0
          \ deta           / \ deta         /        \ deta           /     
> Eq3 := diff(phi(eta), eta, eta)+Sc*f(eta)*(diff(phi(eta), eta))+Nt*(diff(theta(eta), eta, eta))/Nb-Sc*sigma[1]*(1+delta*theta(eta))^n*exp(-E/(1+delta*theta(eta))) = 0;
       /  d   /  d           \\               /  d           \
       |----- |----- phi(eta)|| + 0.01 f(eta) |----- phi(eta)|
       \ deta \ deta         //               \ deta         /

                        /  d   /  d             \\
          + 1.000000000 |----- |----- theta(eta)||
                        \ deta \ deta           //

                                      n    /          E         \    
          - 0.001 (1 + 0.1 theta(eta))  exp|- ------------------| = 0
                                           \  1 + 0.1 theta(eta)/    
> bcs1 := f(0) = 0, (D(f))(0) = 1, (D(f))(blt) = 0, (D(D(f)))(blt) = 0, theta(0) = 1, theta(blt) = 0, Nb*(D(phi))(0)+Nt*(D(theta))(0) = 0, phi(blt) = 0;
   f(0) = 0, D(f)(0) = 1, D(f)(5) = 0, @@(D, 2)(f)(5) = 0, theta(0) = 1, 

     theta(5) = 0, 0.01 D(phi)(0) + 0.01 D(theta)(0) = 0, phi(5) = 0
> NULL;
> L := [.1, .5, 1.0, 5.0, 10.0];
                         [0.1, 0.5, 1.0, 5.0, 10.0]
> 
> for k to 5 do R := dsolve(eval({Eq1, Eq2, Eq3, bcs1}, M = L[k]), [f(eta), theta(eta), phi(eta)], numeric, method = bvp[midrich], maxmesh = 4096, abserr = 0.1e-1, output = listprocedure); X1 || k := rhs(R[3]); X2 || k := rhs(R[4]); X3 || k := rhs(R[5]); Y1 || k := rhs(R[6]); Y2 || k := -rhs(R[7]); Z1 || k := rhs(R[8]); Z2 || k := -rhs(R[9]) end do;
Error, (in dsolve/numeric/bvp/convertsys) too few boundary conditions: expected 10, got 8
> R;
                                      R

> print([(X1 || (1 .. 5))(0)]);
                  [X11(0), X12(0), X13(0), X14(0), X15(0)]
> print([(X2 || (1 .. 5))(0)]);
                  [X21(0), X22(0), X23(0), X24(0), X25(0)]
> 
> print([(Y1 || (1 .. 5))(0)]);
                  [Y11(0), Y12(0), Y13(0), Y14(0), Y15(0)]
> print([(Y2 || (1 .. 5))(0)]);
                  [Y21(0), Y22(0), Y23(0), Y24(0), Y25(0)]
> print([(Z1 || (1 .. 5))(0)]);
                  [Z11(0), Z12(0), Z13(0), Z14(0), Z15(0)]
> print([(Z2 || (1 .. 5))(0)]);
                  [Z21(0), Z22(0), Z23(0), Z24(0), Z25(0)]
> NULL;
> NULL;
> plot([X1 || (1 .. 5)], 0 .. blt, labels = [eta, (D(f))(eta)], thickness = 1, color = black);
Warning, unable to evaluate the functions to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

> plot([Y1 || (1 .. 5)], 0 .. blt, labels = [eta, theta(eta)], thickness = 1, color = black);
Warning, unable to evaluate the functions to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

> 
> plot([Z1 || (1 .. 5)], 0 .. blt, labels = [eta, phi(eta)], thickness = 1, color = black);
Warning, unable to evaluate the functions to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

> 

Hello, 

For a few days Maple crashs everytime i try to use the command "plot3d()". 

I had'nt this problem befor and I have no idea what the reason could be. It ist irrelevant what Funktion I try to visualize,  the window just get closed evertime.

I hope someone can help me.

Thank you!

Tom

I have the following simple dsolve code:

ode := diff(y(x), x, x) + (n*pi)^2*y(x) = A^3*sin(n*pi*x)^3;
dsol1 := dsolve({ode, y(0) = 0, y(1) = 0}, y(x));

How to declare that n is an integer here?

Thanks!

The question reads as follows: 

a) Use Maple to find the linear function f(x,y) = c +mx + ny, whose graph passes through the points P=(5,1,0), Q=(-1,7,1), and R = (2,1,4). 

b) Use Maple to graph the linear function whose formula you found in part a.

c) Use Maple to find the (exact) area of the triangle having P,Q and R as vertices. 

Thank you so much. I am new to Maple and just trying to figure it all out.  

I am trying to write a program for the representation of a function F with three variables i, j, x, where i, j are non-negative integers and x is real. I have tried with different indexings like F(i,j,x), F[i](j,x), F[i,j](x), but did not get the desired result. My objective is to get F as a function of i, j, x.

Here is my maple code with detailed description of my problems in each step.

test1702.mw

Hi

I am plotting a system of equations as follows

restart; with(plots);
sys := {p+x+.6*y-15, p+.3*x+.2*y-10, p+.5*x+y-14};
solve(sys, {p, x, y});
implicitplot3d(sys, x = 0 .. 10, y = 0 .. 10, p = 0 .. 10, shading = zhue, orientation = [-160, 50]);

All is well, but I would llike to see the solution clearly, something along the lines of https://www.maplesoft.com/applications/view.aspx?SID=4093&view=html

How do I set up the spacecurve bit?

Also, I would like to plot the plane formed by the solutions how do I combine that into the plot I already have.

thanks!

Need some help building a worksheet where I try to populate a listbox with items from a package / library.

Starting point could be the "Simply Supported Beam Design" example from MapleSoft Application Center (https://www.maplesoft.com/applications/view.aspx?SID=154265).

In this example, profiles are defined from the AISC Steel Shapes Database.

1. I would like to have a list of the defined properties of the package / library.
2. I would like to have a list of all available profiles.

Plot the Vector Field f(x,y,z) = (-x,y,cos z).

How do I go about plotting this? I have tried everything and nothing is working for me. 

I'm not sure what I am doing wrong. I was given 2 matrices to find B^T*(M^-1)*B where M is a 10x10 matrix and B is a 10x5 Matrix. 

Matrix(10, 10, [[1, 0, 0, 0, 1/2, 0, 0, 0, 0, 0], [0, 1/2, 0, 0, 0, 0, 0, 1/3, 0, 0], [0, 0, 1/2, 0, 0, 0, 0, 0, 1/3, 0], [0, 0, 0, 1/3, 0, 0, 0, 0, 0, 0], [1/2, 0, 0, 0, 1/3, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 1/3, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 1/4, 0, 0, 0], [0, 1/3, 0, 0, 0, 0, 0, 1/4, 0, 0], [0, 0, 1/3, 0, 0, 0, 0, 0, 1/4, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 1/4]]);

M := rtable(1 .. 10, 1 .. 10, [[1, 0, 0, 0, 1/2, 0, 0, 0, 0, 0], [0, 1/2, 0, 0, 0, 0, 0, 1/3, 0, 0], [0, 0, 1/2, 0, 0, 0, 0, 0, 1/3, 0], [0, 0, 0, 1/3, 0, 0, 0, 0, 0, 0], [1/2, 0, 0, 0, 1/3, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 1/3, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 1/4, 0, 0, 0], [0, 1/3, 0, 0, 0, 0, 0, 1/4, 0, 0], [0, 0, 1/3, 0, 0, 0, 0, 0, 1/4, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 1/4]], subtype = Matrix); 'M';
                               M
Matrix(10, 5, [[0, 0, 1/3, 0, 0], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [0, 1/4, 0, 0, 0], [0, 0, 1/4, 0, 0], [0, 0, 0, 1/4, 0], [1/2, 1/2, 1, 0, 0], [1, 1/2, 1/2, 1, 0], [0, 1, 1/2, 1/2, 1], [0, 0, 1, 1/2, 1/2]]);

B := rtable(1 .. 10, 1 .. 5, [[0, 0, 1/3, 0, 0], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [0, 1/4, 0, 0, 0], [0, 0, 1/4, 0, 0], [0, 0, 0, 1/4, 0], [1/2, 1/2, 1, 0, 0], [1, 1/2, 1/2, 1, 0], [0, 1, 1/2, 1/2, 1], [0, 0, 1, 1/2, 1/2]], subtype = Matrix); 'B';
                               B
`~`[`*`](`~`[`*`](B^%T, 1/M), ` $`, B);
Error, dimension bounds must be the same for all container objects in an elementwise operation

We can say that prime p is a partition prime of n if there is at least one prime partition of n having p as least part. Example 8=3+5 so 3 is a partition prime, but 5 is not.
 

Furthermore, say that p is a singular partition prime if there is one and only one partition of n with p as least part. I am trying to find numbers n for which the set Q(n) of singular partition primes is {phi}. That is to say, if we take any prime partition of n, then there are at least two partitions associated with its least member. I find so far only two examples: 63 and 161. Clearly no such n can be prime because then n is a singular partition prime of itself (Incidentally, primes having only themselves as singular partition primes are: 2,3,7,13,23,31,41,79,101,107,149..).

I am asking for a code to compute more terms for the case Q(n)={phi}. 

Thanks in advance

David. 

Dear all, 

I am new to Maple, bit I am getting familiar with it. My current interest is in Lie Symmetry analysis of systems of PDE's such as mass, momentum, and energy balance equations in fluid dynamics, for example.

I have carried out a Lie Symmetry analysis for a single PDE's -- attached herewith. I know it to be correct because it is already published. What I want to do now is how to write down more than 1 PDE's and carry out a similar Lie Symmetry analysis to this system of PDE's. I have been unable to find an example on the web.

Can you advise please - especially how to enter a systems of PDE's as apposed to just one PDE.

Thanks

Nadeem Malik

Sheet9.mw

Hi

i have vector like

Lambda := Vector[row](4, {(1) = *x[1]+.987329667009504589705605634543*x[2]-.436101517853804204529756811668*x[3]+.129526898189138290691985450717*x[4]-u[1], (2) = -1.90060788691970327587464611304*x[1]+.73919873843577037910167897814*x[2]+.20740624192749301981208286864*x[3]-0.45997093443560123039115733769e-1*x[4]-u[2], (3) = -.90583425825293904132909638808*x[1]-.63522563126813377056284195404*x[2]+.28771415757341360851847098373*x[3]+.25334573194765920337346735830*x[4]-u[3], (4) = -1.09745135827817665891430946698*x[1]+.3569463183370574834390092104*x[2]-1.0463476282509211343139483201*x[3]+.78685266819204030978924857651*x[4]-u[4]})

i want to set 2 column till 4 equal zero, i mean

-1.90060788691970327587464611304*x[1]+.73919873843577037910167897814*x[2]+.20740624192749301981208286864*x[3]-0.45997093443560123039115733769e-1*x[4]-u[2]=0

-.90583425825293904132909638808*x[1]-.63522563126813377056284195404*x[2]+.28771415757341360851847098373*x[3]+.25334573194765920337346735830*x[4]-u[3]=0...

how can i do it

thanks