Hello! Hope everything fine with you. I am try to solve the three-point differential by numerical method in attached file but failed. Please see the attachement and solve my problem. I am very thankful your kind effort. Please take care.

three-point.mw

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

What code to use in a maplet that would open the open file window? 

ie when I push a button the open file window opens allowing the user to choose a file.

Download maple_soft.mw

I found the solution of P(r) at P(0)=0, but could obtain the result of v(r) at v(1014030)=-0.4283, v(r) may have a graph such that i can goes from -0.4283 to 0.

I've posted  a code for a procedure, but it seems it does not execute with other versions of Maple the way it does with my version (5.4). Here is the link. I am curious to identify within the procedure code the cause of this divergence. Note: in 5.4 the code is preceded with readlib(eliminate).
 

I Have a problem with opening a document

When i open the document it shows a blanc document and this message 

no matter wich i choose, maple just shows a plain document.

Is there any solution to this problem?

Hi guys, i'm working on a problem using the runge kutta method (4). The problem is, when I run this program, with the given starting values, it perfectly calculates every value of x[w-i] except x[0] and i have no idea why. Can anyone help? Thanks!

Hello.

In maple you can use the interactive menu and choose what axes you want to look at the plot, ie. I choose Z axes and see the plot from above, can i do this by using syntax as a plot option?

Download stokes_start.mw

I have a vector x of this type:

x :=Vector[column]([A__11*u__1+A__12*u__2+...+A__1m*u__m,

A__21*u__1+A__22*u__2+...+A__2m*u__m,

...,

A__n1*u__1+A__n2*u__2+...+A__nm*u__m]);

If I define u:=Vector[column]([u__1,u__2,...,u__m]), then it is clear that the equation has the form x=A*u.

I want to extract the matrix A, for the given vectors x and u.

IMPORTANT: I know I could create a loop (i=1 to m) and set u__i=1 and all other u__j=0 (for all j not equal i) and then reconstruct each column by this method, but it seems to be a overkill for such an easy problem.

I would be glad, if someone could show me a method how one can achive this in maple.

I am trying to reproduce some systems for multibody systems by using maple.

I want to define a rotation matrix

Rot__alpha:=Matrix([[ cos(alpha) , -sin(alpha) , 0 ], [ sin(alpha) , cos(alpha) , 0 ], [ 0 , 0 , 1 ]]);

Now I want to differentiate this matrix with respect to time t. I have tried the following but it does not work:

map(implicitdiff, Rot__alpha,alpha);

I actually want something like this:

d/dt (Rot(alpha))=Matrix([[ -sin(alpha)*a , -cos(alpha)*a , 0 ], [ cos(alpha)*a , -sin(alpha)*a , 0 ], [ 0 , 0 , 0 ]]);

where a is d(alpha)/dt. So alpha is not an explicit function of t.

It would be great if someone could provide me a method for performing these kind of calculations. 
           

In C++, we use continue statement to continue with the loop. What is maple's equivalent statement of continue as used in C++?

Thanks:

Knowing that log[12](27) = a. How can I find log[36](24) in term a in Maple?

I need yours hepl.  I work with the physics paсkage and I set:

with(Physics)

Setup(mathematicalnotation = true)

 Coordinates(X)

Setup(Dgammarepresentation = standard)

Setup(spaceindices = uppercaselatin)

Define(m, m5, y, p, mm, pp)

I try to square the next value: 

W := Dgamma[mu]*d_[mu]+M+Psigma[A]*aa[A]-mm*Dgamma[0]-m5*Dgamma[0]*Dgamma[5]+I*Dgamma[5]*Psigma[B]*pp[B]+I*Dgamma[5]*y

("*" is multiplication)

W*W

And after that I want to simplify it:

Simplify(W*W)

I guess that matter is owing to d_[`~mu`]. If I remove this term:

E:=Psigma[A]*aa[A]-mm*Dgamma[0]-m5*Dgamma[0]*Dgamma[5]+I*Dgamma[5]*Psigma[B]*pp[B]+I*Dgamma[5]*y

And if i do:

E*E

Then next error emerges:

What is it?

here is my homework, i am new to maple and need some help!

Consider the logistic equation:

,

with seasonally varying population constraint

M(t) = 20 + 1.54sin(pi t / 6)

use maple to plot the solution curve for P(0) = 25

Given 3 surfaces:

,(the xy base plane) and 

To plot these I suggested to use cylindrical coordiates knowing x=r*cos(t) and y=r*sin(t)

Which leads to 

However i got problems knowing how to plot this object and dearly ask for help.

plothelp.mw

Hi!

I am an error with the use of the function "Analytic" of the packpage RootFinding. These are the procedures:

CreaCos := proc (C, n, m, t) local k, F; F := C[1][1]+(C[1][2]-C[1][1])*t; for k to n-1 do F := F, C[k+1][1]+((1/2)*C[k+1][2]-(1/2)*C[k+1][1])*(1-cos(Pi*m^k*t)) end do; return F end proc;

Then, for k=50, 100, 150... the instruction

works correctly. However, for higher values of k (for instance, k=250) returns the below error. Some idea or suggets about occurs this error?

Many thanks for your time! 

Error, (in RootFinding:-Analytic) unable to evaluate `@`(evalf, proc (x) option remember; table( [( 0.524900000000000000000000000000e-1+Float(undefined)*I ) = Float(undefined)+Float(undefined)*I ] ) 31250*Pi*sin(62500*Pi*x)/(7/18-(1/2)*cos(62500*Pi*x)) end proc) at the value 0.524900000000000000000000000000e-1+Float(undefined)*I. The expression to be solved was probably not analytic.