I am attempting to combine the sum of two sinusoids having the same frequency in Maple (e.g. A*cos(10*x) + B*sin(10*x + Pi/3) = M*cos(10*x + phi) )

The long way through this is to break each term up into cartesian components and add them (x = A + B*cos(Pi/3 + 3*Pi/2) , y = B*sin(Pi/3 + 3*Pi/2) ). Following the summation, they can be converted back to polar (M = sqrt( x^2 + y^2 ) and phi = arctan(y/x)).

This is a long process for me and if my graphing calculator...

Hi there,

I am currently working on an ODE model with a diffusion term and I have to analyse it analytically (if possible) with Maple12.

With the dsolve function I have a solution but based on BesselJ and BesselY functions.

My questions are:

• Do you know if it's possible to have a solution without these functions?
• When I add the IBCs conditions, the system is inconsistent, do you understand why?

Hi I'm having trouble generating the mandlebrot set using maple... i only have a very basic knowledge of this programme.
>plot3d([-2, 0.7], [0.7, 3] orientation = [-90,0], grid = [250,250], style = patchnogrid, scaling = constrained, colour = mandelbrot);

The above code is giving me trouble, it keeps coming up with an error missing operator or ';'
I don't see anything wrong with it though....

Hello

Is there a way to change FPS via options in an animated curve?

Budgie

The orginial equation to solve is : x" + F(x) = sin 4t

initial conditions : x(0) = 0, x'(0) =1;

F(x) = { bx, x>= 0 and ax, x<0}

I need to do three different values for b = 1, a = 4, then b = 64, a = 4, finially  b = 36, a = 25,

Thanks for any insite I am new to the program.

Courtney

I'm having a few problems with differentiating in Maple. I have a potential function U given by:

U[c] := (1/2)*r^2+M[1]/r[1]+M[2]/r[2];

U[r] := r^2*(M[1]*M[2]-3)/(2*c^2)+((x(tau)+diff(y(tau), tau))^2+(y(tau)-(diff(x(tau), tau)))^2)^2/(8*c^2)+3*(M[1]/r[1]+M[2]/r[2])*((x(tau)+diff(y(tau), tau))^2+(y(tau)-(diff(x(tau), tau)))^2)/(2*c^2)-(M[1]/r[1]+M[2]/r[2])^2/(2*c^2)-M[1]*M[2]*(1/r[1]+(1/r[1]-1/r[2])*(1-3*mu-7*x(tau)-8*(diff(y(tau), tau)))+y(tau)^2*(M[2]/r[1]^3+M[1]/r[2]^3))/(2*c^2);

Hi All,

Is using the NAG package the right way to go about creating the delaunay triangulation of a set of points in 3 dimensions. My points all lie on a surface in R^3.

Thanks. 

Hi All,

Recently I upgraded to Maple 14 from Maple 13, and I've encountered a strange behavior in Maple 14 that doesn't appear in Maple 13. By the definition of Spherical harmonic functions in Maple,  

SphericalY(2*l, -2*l, theta, phi) = SphericalY(2*l, 2*l, theta, phi).

And this is true in Maple 13, however not Maple 14. Can you reproduce this behavior in Maple 14? Or have I done something else that might be causing this?

Thank you for any help you can provide.

Hello!

I'm trying to solve numerically an ODE system with piecewise. And this piecewise is very important for this task.

This system describes behavior of a pulley with friction. There are some constants: m, c, g, mu and J. Values of this constants are not important.

> sys := m*a(t) = piecewise(a(t) < a0, F0*time, a(t) >= a0, 0), v(t) = diff(x(t), t), a(t) = diff(v(t), t);
> m := 5; F0 := 10; a0 := 5;
> initialconditions := x(0) = 0, v(0) = 0;

What does it mean when I get the message, "Warning, model is not of full rank"?

Thank you.
Janice

Dear all,

Find here a part of my precedent post "2D finite element method" which represent a worksheet for 2D triangular Mesh.A part of the mesh can rotate with regard to the other by modifying theta.

restart: with(geometry): with(plots):

ms := 8: n := 6: mr := 5: theta := (1/180)*(23*Pi):

wrs := (1/1000)*<30.8, 33, 36, 39, 40, 43, 46, 49, 53.25>:

wrr := (1/1000)*<20,23,26,28,30,30.8>:

wt := (1/180)*Pi*<0,10,20,30,40,50,60>:

Good Morning.

I have problem when i want to graph my equations. I received this error: 

Error, (in plots/odeplot) curve is not fully specified in terms of the ODE solution, found additional unknowns {t}

I would like to know how i can solve it

Thanks in advance 
 

Hi, my professor wants our graphs to be labeled (which is reasonable). How can I label the graph I made this way:

>inequal({x[1] >= 0, x[2] >= 0, x[1]-3*x[2] >= -3, 2*x[1]+3*x[2] <= 6}, x[1] = 0 .. 3.5, x[2] = 0 .. 2.5)

Thanks.

Hi all,

I noticed that there are not many applications which deal with the finite element method with Maple. I attached in the file below a code for the magnetostatic probems and which allows the calculation of the magnetic field in a permanent magnet synchronous machine.

mafem.zip

I hope this work will help many users of maple in the numerical analyzis with the finite element method

I posted...

The AddVertex command (in GraphTheory) does not seem to copy any attributes that were set in the original graph.

with(GraphTheory):
G := CompleteGraph(5):
SetGraphAttribute(G,"graphAtt" = 1);
SetVertexAttribute(G,1,"vertexAtt" = 1);
SetEdgeAttribute(G,{1,2},"edgeAtt" = 1);
H := AddVertex(G,6):
GetGraphAttribute(H,"graphAtt");
                     ...