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.

 

 

a1 := 5; b1 := 3; a2 := 3; b2 := 4; a3 := 3; b3 := 7; eq1 := expand((y-2)^2/b1^2+(x-5)^2/a1^2 = 1); 1 2 4 13 1 2 2 eq1 := - y - - y + -- + -- x - - x = 1 9 9 9 25 5 eq2 := expand((y+2)^2/b2^2+(x+1)^2/a2^2 = 1); 1 2 1 13 1 2 2 eq2 := -- y + - y + -- + - x + - x = 1 16 4 36 9 9 Sys := {eq1, eq2}; Sol := [solve(Sys, explicit)]; L := map(proc (t) options operator, arrow; eval([x, y], t) end proc, simplify(`~`[fnormal](evalf(Sol), 9), zero)):I don't want complex roots. Tank you.

I have solved the following linear system for 6 variables on Maple using the following code:

sys := {w = -2*Pi*i*k_2*v + 2*Pi*i*k_2*(4*K^2)/(5*Pi)*u, z = -2*Pi*i*k_1*v + 2*Pi*i*k_1*(4*K^2)/(5*Pi)*u, p*s*x = -2*Pi*i*k_1*u - 4*Pi^2*(k_1^2 + k_2^2)*x + a_1, p*s*y = -2*Pi*i*k_2*u - 4*Pi^2*(k_1^2 + k_2^2)*y + a_2, k_1*z + k_2*w = 0, k_1*x + k_2*y = 0}

solve(sys, [x, y, z, w, v, u])

However, the solution yields z = w = 0, I was wondering if this is correct as I feel that these quantities should not be 0 in the problem which I am studying, is there a way to find a solution which does not involve setting these two to 0?

.

 

 

Hello,

I am looking help for solving this integral equation using the collocation method with 1.5<= x<=3.5 ? 

I have used the successive approximation method and the solution seems to be increasing.

Thanks

Error, (in PDEtools/NumerDenom) invalid input: `PDEtools/NumerDenom` expects its 1st argument, ee, to be of type algebraic, but received {[(s_j*e^sigma*b^m*`σ_m`*y_i-s_j*e^sigma*beta*`σ_i`+s_j*e^sigma*b_ilo+(1/2)*b^(m+2)*`σ_j`*e^sigma*`σ_m`*y_i-b^2*`σ_j`*e^sigma*beta*`σ_i`+(1/2)*b^2*`σ_j`*e^sigma*b_ilo-(1/2)*b^(2*m)*`σ_m`^2*b_j*e^sigma*y_i+(1/2)*b^m*`σ_m`*b_j*e^sigma*beta*`σ_i`-(1/2)*b^m*`σ_m`*b_j*e^sigma*b_ilo+b_j*(s_ilo+(1/2)*(b^2*`σ_i`-b^m*`σ_m`*beta)*_lo)+b_j*(s_jlo+(1/2)*(b^2*`σ_j`-b^m*`σ_m`*beta)*_lo)+s_i*e^sigma*b^m*`σ_m`*y_j-s_i*e^sigma*be... when I use simplify I have this error. please guide me Saberali

I need to take the inverse of a tensor which I have denoted as e[~mu,nu] which is defined by a rather larger Matrix. I had computed this matrix using Mathematica and then simply transferred the resulting matrix by using the calling sequence 

with(MmaTranslator):

which worked swell for transferring said matrix in to Maple. Then using the Physics package I was able to define it as a tensor, with a contravariant and covariant index, respectively. Now, when trying to transfer the inverse of said matrix into Maple to define as a new tensor which I intend to call f[mu,~nu], I get an error saying that the number of free indices on the left hand side does not coincide with the number of free indices on the right hand side. Since, this "new" tensor will really just be the inverse of the matrix which I used to define e[~mu,nu], I was wondering if there was any way in which I can simply compute the inverse of the matrix defining e[~mu,nu] in Maple and then let it be equivalent to f[mu,~nu], afterwhich I would then define it as a new tensor itself. 

 

Research_project.mw

Any help would be greatly appreciated.

I would like to plot a piecewise function. g := x -> piecewise(x < -1, x + 1, x = 1, 4, 2 < x, x^2) with showing the value of 4 using plot command. How to do that? 
As I using the command plot(g, -5 .. 5, -5 .. 5, discont = true, thickness = 3), the value of 4 when x = 1 did not show up.