Hi, I am using Maple 18 and struggling with plotting Newton's Method.

I am wanting use the function f(x)=x^3 +cx + 1 where c is a parameter and uses 100 parameter values between -2 and 0, with 100 iterations of each parameter.

Any help would be brilliant.

Thanks in advance,

Neil

deqn4 := {4*eta^4*(diff(f(eta), eta, eta, eta))*(diff(f(eta), eta))-12*eta^3*(diff(f(eta), eta, eta, eta))*f(eta)+12*eta^4*(diff(f(eta), eta, eta))*(diff(f(eta), eta))-36*eta^3*(diff(f(eta), eta, eta))*f(eta)-24*eta^3*(diff(f(eta), eta))^2+12*eta^2*(diff(f(eta), eta))*f(eta)-36*eta*f(eta)^2-4*eta^3*(diff(f(eta), eta, eta))*(diff(f(eta), eta))-12*eta^2*(diff(f(eta), eta, eta))*f(eta)+12*eta^3*f(eta)*(diff(f(eta), eta))+36*eta^2*f(eta)^2-4*eta^3*(diff(f(eta), eta, eta, eta, eta))-12*eta^3*(diff(f(eta), eta, eta, eta))+12*eta^2*(diff(f(eta), eta, eta))-12*eta*(diff(f(eta), eta))+24*f(eta) = 0};
/
| 4 / d / d / d \\\ / d \
< 4 eta |----- |----- |----- f(eta)||| |----- f(eta)|
| \ deta \ deta \ deta /// \ deta /
\

3 / d / d / d \\\
- 12 eta |----- |----- |----- f(eta)||| f(eta)
\ deta \ deta \ deta ///

4 / d / d \\ / d \
+ 12 eta |----- |----- f(eta)|| |----- f(eta)|
\ deta \ deta // \ deta /

3 / d / d \\
- 36 eta |----- |----- f(eta)|| f(eta)
\ deta \ deta //

2
3 / d \ 2 / d \
- 24 eta |----- f(eta)| + 12 eta |----- f(eta)| f(eta)
\ deta / \ deta /

2
- 36 eta f(eta)

3 / d / d \\ / d \
- 4 eta |----- |----- f(eta)|| |----- f(eta)|
\ deta \ deta // \ deta /

2 / d / d \\
- 12 eta |----- |----- f(eta)|| f(eta)
\ deta \ deta //

3 / d \ 2 2
+ 12 eta f(eta) |----- f(eta)| + 36 eta f(eta)
\ deta /

3 / d / d / d / d \\\\
- 4 eta |----- |----- |----- |----- f(eta)||||
\ deta \ deta \ deta \ deta ////

3 / d / d / d \\\
- 12 eta |----- |----- |----- f(eta)|||
\ deta \ deta \ deta ///

2 / d / d \\ / d \
+ 12 eta |----- |----- f(eta)|| - 12 eta |----- f(eta)|
\ deta \ deta // \ deta /

\
|
+ 24 f(eta) = 0 >
|
/
init4 := { {f(0) = 1, (D(f))(0) = 0, ((D@@2)(f))(0) = 0, ((D@@3)(f))(0.1e-2) = 0} };
{{f(0) = 1, D(f)(0) = 0, @@(D, 2)(f)(0) = 0,

@@(D, 3)(f)(0.001) = 0}}
dsol4 := dsolve(deqn4 union init4, numeric, method=rkf45);

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

Hi,

I used to use windows 32 bit, but I have 64 bit windows now. I installed the MapleSim 6.4 and Maple 18 and try to run the simulations that I created with the 32 bit windows. I have this error 'Unable to compile (rc=1), please try again, and if that fails verify your Windows compiller installation'. Could you please let me know what should I do to eliminate this problem?

Best

Onder

i want to plot a arc that starting angle and finishing angle are variable

but it didn't show the arc there is my code

Hello,

      I would like to solve a system of 9 nonlinear equations, with the constraints on all 9 variables to be that they are nonnegative. How can I do this?

My code is below - I am trying NLPSolve and have tried solve, but am getting stuck.

with(Optimization);

restart; eq1 := 531062-S/(70*365)-(.187*(1/365))*(H+C+C1+C2)*S/N = 0;eq2 := (4/365*(T+C))*S/N-(.187*(1/365))*(H+C+C1+C2)*T/N-(1/(70*365)+1/(5*365))*T = 0; eq3 := (.187*(1/365))*(H+C+C1+C2)*S/N-(4/365)(T+C)*H/N-(1/(70*365)+1/(4*365))*H = 0; eq4 := (.187*(1/365))*(H+C+C1+C2)*T/N+(4/365*(T+C))*H/N-(1/(70*365)+3/(8*365)+.2*(1/365)+.1)*C = 0; eq5 := .1*C-(1/(70*365)+1/(4*365)+1/60+.5)*C1 = 0; eq6 := (1/60)*C1-(1/(70*365)+1/(4*365)+1/210+.5)*C2 = 0; eq7 := .5*C1-(1/(70*365)+1/60+0.1e-2)*CT1 = 0; eq8 := .5*C2-(1/(70*365)+1/210+(1/9)*(0.1e-2*7))*CT2+(1/60)*CT1 = 0; eq9 := N-S-T-H-C-C1-C2-CT1-CT2 = 0; soln := NLPSolve({eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, {C, C1, C2, CT1, CT2, H, N, S, T}, assume = nonnegative);

The problem is The square root of 16-x^2 over the interval [0,-4]  0 being the upper bound, -4 being the lower bound.  I have solved 3/4s of this problem but I don't understand what they mean by "Solve the definite integral exactly by geometry". 

I Want to plot only the intersection between this Surfaces 

z = sqrt(x^2+y^2)

y >= abs(x)

x^2+y^2+z^2 >= 4

x^2+y^2+z^2 <= 9

Anyone knows how ?

Thanks!

Greetings,

Maple 15 allows the following syntax

omega := sqrt(w0^2*(1+((z-zf)/z0)^2))

But it does not allow this one 

omega^2 := sqrt(w0^2*(1+((z-zf)/z0)^2))

Why is this so for functions and variables? Is there any way around this, I am really bugged by this issue. The Maple Math told me that the syntax was invalid so I had type the commands directly.

Hi,

I've encountered a problem with Maple.

Is it possible to create a procedure which both displays some plots and also returns a value?

When I write a simple example:

It plots a circle as expected.

But when I add return:

It returns value 2, but no plot is displayed.

Is it possible to do both, display and return a value?

I used to write two copies of procedure, one for plotting and one for returning value.

But thats not "pretty" solution.

Thanks for help.

I have a formula, which is called the cube law. It requires input between 0 <= x <= 1. When creating  a procedure I run into all kinds of problems and get nothing but error codes. Trying to use assume works, it puts a tilde mark on the variable x, but doesn't restrict the input. Here is the formula the way it is writtin in a book on calculus:

f(x) := x^3/x^3+(1-x)^3 (0 <= x <= 1).

Trying to restrict the input to be restricted as required by the relational statement has been proven fruitless.

Any help would be greatly appreciated, thanks,

Rolf.

Hi there,

I used the implicitplot command to plot an implicit equation as follows:

with(plots):

implicitplot(x^3+y^3-3xy=0,x=-3..3,y=-3..3, scaling=constrained)

but maple13 takes error and specifies {xy} as an extra unknown.

Your help will be appreciated.

Regards

Can maple solve maximization problem like

q := proc (a, b, c) options operator, arrow; .2*b+.1*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.2 b + 0.1 c
w := proc (a, b, c) options operator, arrow; .7*a+.1*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.7 ab + 0.1 c
e := proc (a, b, c) options operator, arrow; .7*a+.2*b*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.7 a + 0.2 b c

with(Optimization)

Maxmize(int(min(100+(.7*a+q)*(1/2), a), q)+int(min(100+(.2*b+w)*(1/2), b), w)+int(min(100+(.1*c+e)*(1/2), c), e)-a-b-c-ab-ac-bc)

Error, (in Optimization:-NLPSolve) cannot convert procedures to piecewise

Thanks alot if you can help me.Urgent! Really appreciate.

Hi,

Anyone knows how to define a ring of polynomials with variables x_1,...,x_n such that x_i*x_j = -x_j*x_i if i \neq j?

I tried using the Physics package and AntiCommutator but the problem is that in that case the variables anticommute with themselves so I have x_i^2=0.

If there is no direct way to this, I guess I could define a procedure that would look at the monomials of a polynomial and order them in lexicographic order and each time it switches two variables with different index it would multiply the monomial by -1.

I also don't know how to do that since I don't know how to look at a specific term in a monomial. For example, if my monomial is x_ix_jx_k, is there a way to find the first two variables and then switch them if LexOrder(variable1,variable2)= LexOrder(x_i,x_j)= false, i.e. if i>j?

Thanks,

M-A

Hi there,

I used the following plot command in maple13.

with(plots):

implicitplot(abs(x)+abs(y)=2, x=-2..2, y=-2..2, scaling=constrained)

But the output(square) is not complete at vertices.In fact at the vertex the graph is like a trapezoid.

Any help will be appreciated.

Yegan

For different reasons I need to ocasionally export a number of Maple worksheet in a folder to pdf files. Is there a way to automate this? I would want that the worksheet is opened, output removed, then executed and eventually exported to pdf. It can take quite a while to do this manually for about 50 worksheets.