I posted a similar question about typesetting:-delayDotProduct. I used two different methods to solve the same vector calculus problem and came up with two seperate problems, though both are using the with(Physics[Vectors]) package. Below is my code showing the error.

Thanks for any help you all can give on this.

Cheers,

Dave
 

Download Different_basis.mw

I ran into this problem when trying to do some symbolic vector calculus using the with(Physics[Vectors]) package. I have uploaded my document (at least I hope it uploads correctly). I use 1-D math input because I learned maple 13 years ago.

Thanks for any help you all can provide.

Cheers,

Dave
 

Download Typesetting-delayDotProduct.mw

ode1a := diff(y1(tt),tt) = round(rhs(odeparm1[1][1]))*y1(tt) + round(rhs(odeparm1[1][2]))*y2(tt) + round(rhs(odeparm1[1][3]))*y3(tt);
ode2a := diff(y2(tt),tt) = round(rhs(odeparm1[1][4]))*y1(tt) + round(rhs(odeparm1[1][5]))*y2(tt) + round(rhs(odeparm1[1][6]))*y3(tt);
ode3a := diff(y3(tt),tt) = round(rhs(odeparm1[1][7]))*y1(tt) + round(rhs(odeparm1[1][8]))*y2(tt) + round(rhs(odeparm1[1][9]))*y3(tt);
sys := subs(y3(tt)=1,[ode1a,ode2a]);
print(DEplot(sys, [y1(tt), y2(tt)], tt = 0 .. 16, y1 = -16 .. 16, y2 = -16 .. 16, color = magnitude, title = `Stable Limit Cycles`, arrows = curve, dirfield = 800, axes = none));

how to mirror the vector field graph mathematically?

mirror the graph about x=0 this line,

so that the graph looked flip

i find curl can do, but how to do ?

I am relatively new to maple and am using it for a multivariate calculus class.

I want to define a function g(x,y) which is the derivative of f(x,y).

I am trying:

which is my function, and

which should be the derivative.

When I try g(x,y) I get

but when I try to put numerical values in the first argument, (eg: f(1,y) I get

"Error, (in g) invalid input: diff received 1, which is not valid for its 2nd argument"
Any tips?

Why do I get different answers for the same command?

 different_answers.mw

Download different_answers.mw

convert(arctan(y, x), abs);
                             / y + I x \
                        -I ln|---------|
                             \|y + I x|/

The output has x and y swapped. Further, for complex x and y it's wrong even if x and y are swapped back, because Maple's definition of arctan has sqrt(x^2+y^2) in the denominator, which isn't the same as abs(x+I*y).

Hello everybody

I'm using discrete distributions from the Statistics package and I found a rather strange result.

In short the theoritical values of some statistics of a NegativeBinomial(1, P) Random Variable (P being the probability of success equal to 1e-4) are correctly computed, but their empirical estimators computed from a sample of this RV are roughly wrong.

For NegativeBinomial(1, P) is similar to Geometric(P) I asked Maple to compute the theoritical values of some statistics of Geometric(P) and next to assess their empirical values from a sample of Geometric(P).
Some discrepancies still remain but they can be explained by statistical fluctuations.

Could you please look to the attached file (an error on my part is still possible) and help me to fix this ?

Thanks in advance


PS : the histogram of Sample(NegativeBinomial(K, P), AnySizeYouWant) is obviously wrong (it should look like a decreasing exponential) 

Download NegativeBinomial.mw

hi..i have a problem with dsolve.

please help me thanks

1.mw
 

Download 1.mw

Hi Maple folks!

I am trying to solve the following PDE in maple:

pde:= diff(c(x,t),t)+c(x,t)*diff(c(x,t),x)=0

bc:=c(x,0)=piecewise(x < 0, 0, x > 0, 1)

I tried the following:

pdsolve({pde,bc}, c(x, 0))

But it gives no solution, i also tried a numeric solution, but i couldn't make it work. Please help.

Thanks! :)

I have written a Maple code which extracts terms in a matrix into two new matrices dependent on certain markers. I am running into problems with this particular code as it requires enormous amounts of RAM when going to large matrix sizes (32GB+); and will take 4-5 days to complete for larger matrix sizes, or simply crashes.

Each matrix element is an algebraic expression containing A and B and the markers hh and ss. The example in the attached Maple file (3 x 3 matrix) is a minimal working example and the actual expressions are much longer and the matrices much larger (5000 x 5000 or larger). I utilise the symmetry of the matrix to ease the process, and use the tril command of MTM to extract only the lower triangular matrix; as this is all that is needed for the next stage (I have yet to test how efficient this command is, but use it for now).

The two markers hh and ss are targeted and if found in an expression that term will be extracted accordingly. The matrix elements are expanded and a procedure run over the matrix and will extract terms with hh into one matrix and terms with ss into another matrix. I think the map(expand,...) part is what causes the memory issue; but I believe this needs to be done to correctly extract terms.

Any help or tips are greatly appreciated on how to increase efficiency and improve this method.
 

Sub_matrix_extraction.mw

-Yeti

Hi, I have a system of 15 polynomial equations (of 15 variables) of 2nd degree. Meaning, every monomial has at most mulplications of 2 different variables or a single variable squared.

I tried to solve them in Maple but it's been 2 hours and Maple is still thinking.. I don't have any indication if Maple is just stuck or it should return a solution at some point.  

Below I wrote the Maple commands that I tried to run (sorry but I'm new here so I didn't know how to write them nicely). I also wrote one of the solutions that I'm expecting to get (15 values for the 15 variables). and you can see that by typing eq7 for example I got -2.5=-2.5 (and the same for the rest of the equations) which means that the values are correct.

I know there is an algorithm for solving any system of numerical polynomial equations of any degree. I don't know if Maple implemented it.

Anyone knows if Maple is able to solve the equations below? and if not, is there any other way that can be done to solve them?

Thanks

David

Here's the commands for Maple (including one of the solutions I'm expecting to get).

The variables are: r11, r12, r13, r21, r22, r23, r31, r32, r33, t1, t2, t3, s1, s2, s3

eq1:=r11^2+r21^2+r31^2 = 1;
eq2:=r12^2+r22^2+r32^2 = 1;
eq3:=r13^2+r23^2+r33^2 = 1;
eq4:=r11*r12+r21*r22+r31*r32 = 0;
eq5:=r11*r13+r21*r23+r31*r33 = 0;
eq6:=r12*r13+r22*r23+r32*r33 = 0;
eq7:=-30*r13-.79382581863774e-1*s1*r11-.95259098236529e-1*s1*r12+.992282273297173*s1*r13 = -.83717247687439e-1*t1;
eq8:=-30*r13+.79382581863774e-1*s2*r11+.95259098236529e-1*s2*r12+.992282273297173*s2*r13 = .76364294519742e-1*t2;
eq9:=-30*r13-.86165283952334e-1*s3*r11+.103398340742801*s3*r12+.990900765451843*s3*r13 = -.81460429387834e-1*t3;
eq10:=-30*r23-.79382581863774e-1*s1*r21-.95259098236529e-1*s1*r22+.992282273297173*s1*r23 = -.107930827800543*t1;
eq11:=-30*r23+.79382581863774e-1*s2*r21+.95259098236529e-1*s2*r22+.992282273297173*s2*r23 = .60269029165473e-1*t2;
eq12:=-30*r23-.86165283952334e-1*s3*r21+.103398340742801*s3*r22+.990900765451843*s3*r23 = .105021268850622*t3;
eq13:=-30*r33-.79382581863774e-1*s1*r31-.95259098236529e-1*s1*r32+.992282273297173*s1*r33 = .990627255252918*t1-30;
eq14:=-30*r33+.79382581863774e-1*s2*r31+.95259098236529e-1*s2*r32+.992282273297173*s2*r33 = .995256820446840*t2-30;
eq15:=-30*r33-.86165283952334e-1*s3*r31+.103398340742801*s3*r32+.990900765451843*s3*r33 = .991128009660183*t3-30;

r11 := 1;
r12 := 0;
r13 := 0;
x := .523598775598299;
r21 := 0;
r22 := cos(x);
r23 := -sin(x);
r31 := 0;
r32 := sin(x);
r33 := cos(x);
t1 := 29.862424638395044;
t2 := 32.737813080348531;
t3 := 30.689747387623886;
s1 := 31.493054789905649;
s2 := 31.493054789905628;
s3 := 29.014005238849720;
eq1;
eq2;
eq3;
eq4;
eq5;
eq6;
eq7;
eq8;
eq9;
eq10;
eq11;
eq12;
eq13;
eq14;
eq15;

unassign('r11', 'r12', 'r13', 'r21', 'r22', 'r23', 'r31', 'r32', 'r33', 't1', 't2', 't3', 's1', 's2', 's3');
solve({eq1, eq10, eq11, eq12, eq13, eq14, eq15, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, {r11, r12, r13, r21, r22, r23, r31, r32, r33, s1, s2, s3, t1, t2, t3});

I did dummy indices implication using add command as below. is it OK or there are mistakes?

Is there another way to imply dummy index summation in maple instead of using add command?

any suggestion???

Download Mesal.mw

 ODE:= -2 sin(1/2 theta(t)) cos(1/2 theta(t)) (diff(theta(t),t)^2-9.8000 sin(theta(t))-(150+4 sin(1/2 theta(t))^2) (diff(theta(t),t,t)=0:

 ICS:=  theta(0) = Pi/6, D(theta)(0) = DthetaZero:

When I use functions from the DocumentTools to display information in a Maple document, they display a weird bracked along the left hand side of the display. I would like them to display in the same way they do in the help pages. I have attached a document showing the results I get when I repeat some of the examples from the help pages.

display-problems.mw