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

Why does Maple give me an answer like this?

How do I force Maple to automatically multiply the exponents in the denominator and eliminate the factor y^2?

Simplify will not do it.

Thank you!

mapleatha

PrimeCounting(6469693230)

It never finishes computing, but smaller values work.

Is there any way to speed it up for larger values(I'm looking to do quite large values.

https://primes.utm.edu/nthprime/index.php#piofx

returns instantly.