MaplePrimes Questions

I need to find the intersection points of two circles (x-7)^2+(y-2)^2=100 and (x-11)^2+(y-5)^2=75 using the modified Newton secant method. 

I have already done this by using these two equations to make one in terms of x. But, i am looking for suggestions on how to make a multiple variable version of the secant method.

Given a sequence for example: 8, 32,128,512,2048 and so on. Suppose we do not know it is 2^(2n +1). How to make Maple to find the principle?

Letter_Lrotate2D.mws

   I've asked a similar question before, and got the letters to rotate in 3D using the rotate command.  The resulting size in a 3D plot is relatively small..  The program gives a plotted letter L.  I wish to get some animation of it being rotated about the y-axis - but remaining in 2D.  ie looking at the letter it appears to be rotating.  To do thisI have multiplied the x-coords by the cosine of the angle of rotation.   

   In the for loop towards the end of the program there is a variable c1 which is incremented by one in the loop - yet it remains at 2, regardless of the angle increments.   Why is that?

  Again, any help or comments much appreciated.

 

Download test.mw
Udp: document attached.

Q__gr := Typesetting[delayDotProduct](Vector[row](8, {(1) = 0, (2) = 5, (3) = 10, (4) = 15, (5) = 20, (6) = 25, (7) = 30, (8) = 35}), Unit('m'^3/'day'), true)

Vector[row](%id = 18446746861178193550)

(1)

`η__gr` := Vector[row](8, {(1) = 0, (2) = 9.324, (3) = 17, (4) = 23.232, (5) = 27.192, (6) = 29.6, (7) = 29.516, (8) = 24.592})

Vector[row](%id = 18446746861214593262)

(2)

`pointsη` := [seq([Q__gr[i], `η__gr`[i]], i = 1 .. 8)]

[[0, 0], [5*Units:-Unit(m^3/d), 9.324], [10*Units:-Unit(m^3/d), 17], [15*Units:-Unit(m^3/d), 23.232], [20*Units:-Unit(m^3/d), 27.192], [25*Units:-Unit(m^3/d), 29.6], [30*Units:-Unit(m^3/d), 29.516], [35*Units:-Unit(m^3/d), 24.592]]

(3)

CurveFitting[PolynomialInterpolation](`pointsη`, q)

0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d)

(4)

eta := proc (q) options operator, arrow; 0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d) end proc

proc (q) options operator, arrow; 0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d) end proc

(5)

eta(12*Unit('m'^3/'day'))

.8415811058*Units:-Unit(m^3/d)^7/Units:-Unit(m^3/d)^7-8.900090268*Units:-Unit(m^3/d)^6/Units:-Unit(m^3/d)^6+36.88243199*Units:-Unit(m^3/d)^5/Units:-Unit(m^3/d)^5-75.97578241*Units:-Unit(m^3/d)^4/Units:-Unit(m^3/d)^4+80.63761922*Units:-Unit(m^3/d)^3/Units:-Unit(m^3/d)^3-45.58150401*Units:-Unit(m^3/d)^2/Units:-Unit(m^3/d)^2+31.86194286*Units:-Unit(m^3/d)/Units:-Unit(m^3/d)

(6)

simplify(19.76619849)

Error, (in Units:-Standard:-+) the units `m^21/s^7` and `m^18/s^6` have incompatible dimensions

 

``


 

Download test.mw

 

Hi!

I got a function (from CurveFitting), that produce a polynomial with the some units inside:

 

As i can see, it can be easy simplified, but i get:


 

Is where any trick to do it, without stripping units with convert(unit_free) or something?

Thank you!

 

Dear, I want to plot the graph for differerent values of F=10,20,30 one after the other  in the title automatic it should take the values of F and shown in the following sample here is my codes,  and sample graph

restart:
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
lambda:=Int(K1,z=0..1):
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plot( [seq(eval(lambda, Nb=j), j in [0.1,0.2,0.3])], delta2=0.02..0.1);
 

Letter_L_rotate.mws

The attached program is of the letter L rotating.  about an axis which is the outer side of the 'upright' of the L.  I want to give the illusion that the letter L is rotating about a y-axis, where the x, y axes are in the plane of the screen, with the x-axis going from left to right.  The program at present 'works', in the sense that the letter L is rotating about an axis.  ...but it's not quite what I want.  My eventual aim is to plot the letters A and L side by side, AL, and as the letter L rotates the leg of the L would bump into the A,  so I want the A to rise in sync with movement of the letter L.   ...but that's a long way off!  For the moment I just want to know how to orientate the letter L, so that it is in its usual L state

   I think I have two options: either to change the coordinate positions of L, (long and tedious..) or quicker might be to change the axis of rotation.  Any suggestions on making the animation smoother would also be appreciated.

Thanks, David

 .

In the fraction 3^665/2^x=y, where y is more than 1 and less than 2, I want to know the value of x. What input should I write, please?

  1. (k-2)*(k^2+5)*(k^3-k^2+7*k+8)/(6*k*(k^2-3*k+8))
  2. (k+2)*(k^2+5)*(k^3-5*k^2+13*k-8)/(6*k*(k^2-3*k+8))
  3. k^3+3*k^2+11*k-3
  4. k^2+2*k+9

How to draw the following function containing an imaginary number
 

NULL

u := proc (x, t) options operator, arrow; 1/2+I*sqrt(2)/(exp(x-t)+2*exp(-x+t))+(1/2)*(exp(x-t)-2*exp(-x+t))/(exp(x-t)+2*exp(-x+t)) end proc

proc (x, t) options operator, arrow; 1/2+I*sqrt(2)/(exp(x-t)+2*exp(-x+t))+(1/2)*(exp(x-t)-2*exp(-x+t))/(exp(x-t)+2*exp(-x+t)) end proc

 

``

 

 

 

``

 

 

 

 

 

 

``

 

 

``

NULL


 

Download plot33.mwDownload plot33.mw

I'm attempting to solve the complex-valued differential equation

restart:
assume(t::real):

pde := diff(A(t),t) - I * conjugate(A(t))*A(t)^2 = 0:
dsolve(pde);

However, it seems Maple attempts to solve it using separation of variables and gives

t-Intat(-I/conjugate(_a)/_a^2,_a = A(t))+_C1 = 0

Unless I'm mistaken, a complex integral such as this isn't even defined without a contour.

Working the integral out by hand, I know the solution to be

_C1*exp(I*abs(_C1)*t)

Is there a trick I'm missing to get Maple to find this solution? Or is this outside the scope of what Maple can handle on it's own?

Thanks!

I believe I am having memory issues which is causing the MAPLE kernal to terminate.  If I assign m a value then it seems to work, but I would like to leave m unassigned so that S1 can later be evaluated for any arbitrary m.  Is there a way around this?


 

T := M*tau;

`assuming`([simplify(expand(combine(S1)))], [m::integer])

``


 

Download MAPLE_crashing.mw

Hello everybody

I am currently having trouble 'manually simplifying' an equation that I differentiated using maple. Basically, the equation that I get involves hessian-matrices. Since the resulting equation is rather long, I would like to replace every such matrix by a sign, say H.
This is easy to do for non-matrix equations, i.e.

subs(sin(x)sin(y) = z, sin(x)sin(y) + xy)

will give me z + xy, which looks a bit easier.

However, I couldn't figure out how to do a similar thing with matrices. For example, the following code

subs(Matrix(3, 3, [[x, y, z], [y, z, x], [z, x, y]]) = A, Matrix(3, 3, [[x, y, z], [y, z, x], [z, x, y]])+B)

 will not work.
I already thought about converting the matrices first into lists (because for some reason it works for lists). However, I would also like to do the same as above for for functions of matrices, i.e. set F(A) = G, where A is a matrix (this also fails, probably for the same reason).

Also, the same thing seems to fail for vectors and general arrays, so I guess the actual problem might be the array type.

I also tried alternative ways, such as eval or algsubs, but the later even gives an error since it cannot handle matrices at all.

This feels a bit like a noob question, but I spent almost 2 days now searching for an answer or a workaround, so my apologies if I missed something trivial.

All the best
Adrian

 

 

Hi, I try to solve the below integral. when I press enter key maple dosen't show answer and show the integral again.

int(r*r[bc]*r[tc], r = r[bc] .. r[tc]);

but when I write intgral this way and use " i " as subscript ,maple solve it.

int(r[i]*r[bc]*r[tc], r[i] = r[bc] .. r[tc]);

I just want to know why?

what is difference between first and second integral?

and also is there any way (or any packages) to solve these integrals?

(I read https://www.maplesoft.com/applications/view.aspx?SID=6846&view=html article befor)

thanks.

Good day everyone,

I'm relatively new to Maple, and I study engineering (not maths), so maybe my problem is related to a bad approach of the equations.

The problem is that I'm asking Maple to do a double integral of a "complex" expression and it only solve the first integral, but not the second.

Here is the code:

https://www.dropbox.com/home/C%C3%81LCULO%20FUERZAS

Thanks.

Dears!

Hope everyone should be fine. I am face to import excel sheet in maple. I saved execl sheet with name "Employees.xlsx" at desktop. When I use the following command
S := Import("Employees.xlsx", 1, "A1:B101");

I got the following error. 

"Error, (in ExcelTools:-Import) Could not open the file"

Please help me to fix this problem. I am waiting positive response. 

Special request to:

@Carl Love
@Preben Alsholm
@acer

First 863 864 865 866 867 868 869 Last Page 865 of 2428