I want to plot the solutions of the equation (x-y)^2+(1-z)^2=0.

However, implicitplot3d is not able to plot them, at least using the default arguments. Any recommendations?

I know a priori that it is going to be a curve contained in a plane in case that makes this task easier.

The final solution is the sum of 2 terms and I am doing it individually for 1st term. I think I am doing it wrong because answer did not match when comparing with textbook answer. Can anyone teach me or hint me compute the final series for m=1 to 10. An example of final series for m=1 would be helpful. 

Download ODE_sol.mw

Hello,

Attached is a simple worksheet dealing with units and angular motion.   I need to reverse calculate from an acceleration to the uniform motion (rad/sec or RPM) which is used in forces and FEA simulation.

The initial units work fine and get to a velocity (feet/sec) but then fail to convert to rad/sec or RPM.   The attached shows the results of units Hz (1/sec) and the problems with resolving to RPM or rad/sec.   the unit (length / length (radius)) as an answer stil confounds me.

Uniform_circular_motion_units_issue.mw

Any help greatly appreciated,
Bill

Hello, I have the following (nested) summation:

add(binomial(196, j)*0.5^(32.1+j)*add(binomial(109, l)*(-1)^(j+l)*add(0.5^(3.1*h)*GAMMA(-4.1+h)*GAMMA(9.3170731707317073170731707317073170731707317073170+.24390243902439024390243902439024390243902439024390*l+.24390243902439024390243902439024390243902439024390*j+.75609756097560975609756097560975609756097560975610*h)/(factorial(h)*(10.354838709677419354838709677419354838709677419355+.32258064516129032258064516129032258064516129032258*j+h)*GAMMA(.24390243902439024390243902439024390243902439024390*j+.24390243902439024390243902439024390243902439024390*l+15.417073170731707317073170731707317073170731707317+.75609756097560975609756097560975609756097560975610*h)), h = 0 .. infinity), l = 0 .. 109), j = 0 .. 196)

But it takes ages to sum h up to 80, for example. I was wondering whether we could improve to gain time. any ideas?

CorrBasis := {<sqrt(3),2,0>, <0,1/5,1>};
CorrMat := Matrix([op(CorrBasis)]);

Basis := {<3^(0.5), 2., 0.>, <0.,0.2,1.>}; #float
Basis := remove(i -> is(LinearAlgebra:-Rank(<CorrMat|i>) <> nops(CorrBasis)),  Basis);

I need to use matrices/ vectors with different datatypes. Do I need to turn them all into floats in order for Basis to remain the same after calling

Basis := remove(i -> is(LinearAlgebra:-Rank(<CorrMat|i>) <> nops(CorrBasis)),  Basis);

Is there any better way to manage this?

How do you stop Maple during this code execution:

proc()
x:
print(1);
goto(x);
end():

     Hello everyone !

     I have a problem asking for help:

     In the Oxy coordinate plane, for rectangles are limited by straight lines: x=1, x=7, y=1, y=9 and there are 63 points distinguished from coordinates that are integers located on this rectangle.

     These include:

• 7 black points with coordinates are listed in the list:

[[1,1], [2,1], [3,1], [4,1], [5,1], [6,1], [7,1]].

• 7 red points with coordinates are listed in the list:

[[1,2], [2,2], [3,2], [4,2], [5,2], [6,2], [7,2]].

• 8 yellow points with coordinates are listed in the list:

[[1,3], [4,3], [5,3], [7,3], [1,4], [4,4], [5,4], [7,4] ].

• 6 pink points with coordinates are listed in the list:

[[2,3], [3,3], [6,3], [2,4], [3,4], [6,4]].

• 8 brown points with coordinates are listed in the list:

[[1,5], [3,5], [5,5], [7,5], [1,6], [3,6], [5,6], [7,6]].

• 6 purple points with coordinates are listed:

[[2,5], [4,5], [6,5], [2,6], [4,6], [6,6]].

• 9 blue points with coordinates are listed in the list:

[[1,7], [2,7], [7,7], [1,8], [2,8], [7,8], [1,9], [2,9], [7,9]].

• 6 green points with coordinates are listed:

[[3,7], [5,7], [3,8], [5,8], [3,9], [5,9]].

• 6 orange points with coordinates are listed in the list:

[[4,7], [6,7], [4,8], [6,8], [4,9], [6,9]].

     Help me find the integer coordinates of the 63 points when arranging them on the rectangle knowing that their HorizontalCoord has not changed, and the VerticalCoord of the points of the same color is always different with the Maple command.

     Thank you so much for your help!

Here's a procedure using GraphTheory by morphing one four letter word into another by changing only one letter at a time.  This is my initial working version.  I've commented out the DrawGraph portion as it takes a long time (5 minutes or so) to produce.  Using the Neighbors command from the GraphTheory package the graph can be shrunk to only include the relevant paths and will take a shorter time to draw.  It's an initial version so there is room for improvements.
 

Download Word_Morph_3.mw

Word_Morph.maple

I would show you what I have but I actually don't know where to start this problem..

If n people (numbered 1 to n) stand in a circle and someone starts going around the circle and eliminating every other person till only one person is left, the number J(n) of the person left at the end is given by 

    J(n) = 1                           if n = 1
    J(n) = 2*J(n/2) - 1          if n > 1 and n is even
    J(n) = 2*J((n-1)/2) + 1   if  n > 1 and n is odd

(i) Write a recursive procedure to compute J. [As a check the first 16 values (starting with 1) of J(n) are 1,1,3,1,3,5,7,1,3,5,7,9,11,13,15,1]. 
(ii)Compute the value of J(10000). 
(iii) Can you explain why this is so much faster than our recursive procedure to compute the n-th Fibonacci number?

Why this error message ?
with(geometry);
vartheta := (2*Pi)/17;
x || k := cos(k*theta);
y || k := sin(k*theta);
 xk := cos(k theta): yk := sin(k theta)
.M||k:=(point,x||k,y||k);
 with(plots):
Points := pointplot([seq(M || k, k = 1 .. 16)], symbol = solidcircle, color = red, symbolsize = 10):
Error, (in plots:-pointplot) cannot convert data to an Array of datatype float[8]; Thank you for your answer..

Hello, I have a question.

The exercise consists of defining the image in relation to a chosen base. For example, I choose {x^4,x^3*y,x^2*y^2,x*y^3,y^4} as a base. Then I want to calculate T(x^4), T(x^3*y),T(x^2*y^2),T(x*y^3),T(y^4) through maple but but maple won't calculate it.

Do you know what I'm doing wrong?
Thanks in advance

 I came across this question while trying to verify the equality of expressions containing elliptic integrals.

Download EllpticK(1).mw

I have no experience with elliptic integrals. Can I assume in this case that infinity is correct?

Download simulasi_I(t)_kurang_dari_1.mw

I am using the Calculus Study Guide and the mathematical notation is half Maple notation and then half typeset notation. I have tried to tweak with Tools>Options>Display settings but it stays the same

For example:

{a[n]}[n = n[0]]^infinity;, where 
                             "a[n]"

 is the general term for the sequence, and 
                             "n[0]"

 is the starting index.
The notation is often shortened to {a[n]}; to save printing costs.

...and the math notation past the word "where" is in typeset notation but changed to Maple Notation when I cut and pasted.

Thank You in advance.

I'd like to plot the graph in polar coordinates of the function r=f(t)=2/cos(t), which should be a vertical line parallel to the y-axis. I've tried it with two different commands, and below you can see the results

plot(2/cos(theta), coords = polar)

What is going on here?

I then tried using polarplot from the plots package.

with(plots):

polarplot(2/cos(t), t = 0 .. 2*Pi)

Shouldn't the graph be a vertical line at 2?

Notice that if I plot two similar plots

polarplot([1/cos(t), 2/cos(t)], t = 0 .. 2*Pi)

The same result is obtained

What is going on here?