Hello,

So I want to create an image of the cusp catastrophe that looks like this 

but instead I have been getting this image where the discontinuity is plotted out to infinity

I have tried to split up the surfaces on either side of this but I haven't been able to display both on the same plot whilst using implicitplot3d. 
This is the line of code for the image, there't not much too it
"
with(plots);

implicitplot3d(x^4+x^2*z+x*y, x = -7 .. 7, y = -50 .. 50, z = -10 .. 5, axes = boxed, style = surfacecontour, grid = [60, 60, 60], orientation = [180, 8, 175], 'transparency' = .1);
"
If anyone could help or point me in the right direction that would be great.
Cheers,
Darryl A.M.

using the code generator assistant I entered the following function

p := proc (z::(float[8]))

local a::integer, accm::(float[8]), k::integer, k1::(float[8]), c;
c := Array([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], order = C_order, datatype = float[8]);
k1 := 1;
c[1] := evalf(sqrt(2*Pi));
a := 12;
for k to a-1 do c[k+1] := evalf(exp(a-k)*(a-k)^(k-1/2)/k1); k1 := -k1*k end do;
accm := c[1];
for k to a-1 do accm := accm+evalf(c[k+1]/(z+k)) end do;
accm := accm*evalf(exp(-z-a)*(z+a)^(z+1/2));
return accm/z
end proc

the code-generated julia code follows

function input(z)
c = [0,0,0,0,0,0,0,0,0,0,0,0]
k1 = 1
c[0] = (sqrt(2 * pi))
a = 12
for k = 1:a - 1
c[k] = (exp(a - k) * (a - k) ^ (k - 1//2) / k1)
k1 = -k1 * k
accm = c[0]
for k = 1:a - 1
accm = accm + (c[k] / (z + k))
accm = accm * (exp(-z - a) * (z + a) ^ (z + 1//2))
return(accm / z)
end

two things are wrong

1: no end after loop end

2: array index starts at 0, it should be 1 and of course the array references should reflect that

btw, it would be nice to be able to enter code tags like [code] code here [/code]

Hello everyone.  Let n be a positive number; I seek how to write the matrix Z\in M_{2^n,n} that is defined as follows:

for every i,j, Z_{ij}\in \{-1,1\}.

When I consider the 2^n rows of Z, I want to find ALL the possible sequences of length n with the entries +-1. A toy example, when n=2, is Z=\begin{pmatrix}1&1\\1&-1\\-1&1\\-1&-1\end{pmatrix}.

Thanks in advance;

The method of solving underdetermined systems of equations, and universal method for calculating link mechanisms. It is based on the Draghilev’s method for solving systems of nonlinear equations. 
When calculating link mechanisms we can use geometrical relationships to produce their mathematical models without specifying the “input link”. The new method allows us to specify the “input link”, any link of mechanism.

Example.
Three-bar mechanism.  The system of equations linkages in this mechanism is as follows:

f1 := x1^2+(x2+1)^2+(x3-.5)^2-R^2;
f2 := x1-.5*x2+.5*x3;
f3 := (x1-x4)^2+(x2-x5)^2+(x3-x6)^2-19;
f4 := sin(x4)-x5;
f5 := sin(2*x4)-x6;

Coordinates green point x'i', i = 1..3, the coordinates of red point x'i', i = 4..6.
Set of x0'i', i = 1..6 searched arbitrarily, is the solution of the system of equations and is the initial point for the solution of the ODE system. The solution of ODE system is the solution of system of equations linkages for concrete assembly linkage.
Two texts of the program for one mechanism. In one case, the “input link” is the red-green, other case the “input link” is the green-blue.
After the calculation trajectories of points, we can always find the values of other variables, for example, the angles.
Animation displays the kinematics of the mechanism.
MECAN_3_GR_P_bar.mw 
MECAN_3_Red_P_bar.mw

(if to use another color instead of color = "Niagara Dark Orchid", the version of Maple <17)

Method_Mechan_PDF.pdf

Write a Maple procedure which takes any cubic polynomial [−3, 3] p −→ R and returns the following:

How i can plot of p and its tangent lines illustrating the nature of each critical point.

Demonstrate this in the instance that p(x) = 1/3 x^ 3 − 6/10 x^ 2 − 187/100 x + 1.

Any ideas how i can solve this?

I wonder if there is any way to sort expressions like

 to

also:

how to sort  to

Hello,

I'm sorry to bother you but I have a problem with the numeric resolution of a system of 3 differential equations. The system is as follows  : sysdif :=

As you can see the system is composed of 3 differential equations, and I enter initial conditions in the object "sysd". Then I try a numeric resolution by executing the following command (I give a value to parameters before)  :

Then Maple's answer is : Error, (in dsolve/numeric/process_input) missing differential equations and initial or boundary conditions in the first argument: sysdif.

I can't see where I'm wrong, does anyone notice something that could explain this error message ? There's no help page about this error so I ask the question here.

Thank you very much for your time if you answer this,

Louis

I noticed that Maple 2016 did not add entry in the start menu->All programs, after installing it on window 7, as all the other Maple releases did.  I looked everywhere and do not see it. This is very strange. Only time I had problem like this, where I install Maple but it does not show up in the start menu.

I did not have problem installing 2016, and I can use it fine.  But each time I need to start it now, I have to do start menu->Search and type Maple 2016, to find Maple 2016.

Is there a way to make it show in the start menu? Should I uninstall it and install again? Any one else had this problem on windows?

Using Maple 2016, windows 7, 64 bit. I see 2 problems allready. One is that long display do not wrap around as it does with Maple 2015. Second problem, I see strange characters inside the numbers displayed on the screen.  Here is 2 screens shots, same computation, one from Maple 2015 and one from Maple 2016. Both on same PC. Windows 7.

Notice also the result for 2016 do not wrap. I had to scroll to the right to see the full result.

Here is the code in plain text. Can someone verify if they get the same problem on Maple 2016?

V:=(m,n)->binomial(m-floor(1/2*(n+1)),m-n)+binomial(m-floor(1/2*(n+2)),m-n);
r:=V(2000,500):
r:=r/10^6:
r:=r/(60*60*24*365);

In this course you will learn automatically using Maple course Statics applied to civil engineering especially noting the use of components properly. Let us see the use of Maple to Engineering.

Static_for_Engineering.mw

(in spanish)

Atte.

Lenin Araujo Castillo

Ambassador of Maple

@Markiyan Hirnyk 

Thanks for helping.

(1) But how to find period using code?

(2) how did you determine n=29 ?

(3) I don't need sequence of numbers, but would like to extract the coefficients of exp(..) terms. For example in the follwoing expression, how to extract the coefficients: (note that the expression is missing one more term and i don't know how to fix?)

Hi,

restart;
.123456789*.123456789;
                                               0.01524157875

i need Maple to respond   0.015241578750190521 by adjusting Digits number dynamically. I don't want to adjust  the Digits Number to a fixed value.

Thanks!

OnesM:=Matrix(`%id`=119376536)

Anyone can solve this??

Many thanks!

How can I have a sequence as follow by using "seq" command?

x1,x2,...x200.

Please note that I need "xi" not "x_i".

Hey,

is there a way to substitute funktion combinations when they have the same arguments? I want to substitute abs(exp)*abs(1,exp) with exp. Algsubs works as long as I know the excat expressions. However, I want to do this substitution for any exp. Simplify doen't work here ether.

I also tried to write my own simplification rule with no success:

Download simplify_abs.mw

Thanks in advance!