Hi all I am truing to write a maple program for Newton's Method for finding a zero of a function.  I got pretty far but I am trying to get it to list all the values not just one.  Any help would be great:

Here is  what I got:

0.2000000000

0.2016393443

0.2016396757

What I want is if I type in 2 (like the last line I showed) I want the output to show for n=0, n=1, =2.  Instead it is only showing the output just for 2.  Any help would be great thanks, Matt.  

Hi,

I have a code compute some function : 

         alpha1:(n, m,1) -> (n + 1) (int(K(|m h - y|), y = n h .. (n + 1) h))

              int(K(|m h - y|) y, y = n h .. (n + 1) h)
            - -----------------------------------------
                                  h                    
         alpha2:  (n, m,2) -> -n (int(K(|m h - y|), y = n h .. (n + 1) h))

                int(K(|m h - y|) y, y = n h .. (n + 1) h)
              + -----------------------------------------
                                    h                    
and  I have a Matrix "MatA" .

My aim, when I give the value of the Kernel K used in alpha1, and alpha2, like K(x,y)=ln|x-y| , I want a numerci Matrix.

How can I do it.

Many thinks.

Fred.mw

Hi,

m := Matrix(8, 2, [1, 2, 1, 4, 4,1,4,7,1, 2,4,1,1,4,4,7]);
plot(m)

I would like to add a mame of each point in the graph as: alpha:=(1,2), beta:=(1,4) , gamma:=(4,1),eta:=(4,7)

and the symbol of each point is a bleu sphere.

and how to remove axis borders, tick marks.

Many thinks

Several years ago, I used to plot, in Maple 7 I think, 3D scalar functions by using procedures to create a 3D mesh and populate the data points before I could use plot3d.

I wonder if there is a more convenient (least coding) way to do it today? Consider for example f(x,y,z) = x^2 + y^2+z^2

A way to visualize a number of concentric isosurfaces of f is to loop with (is there a tag to write this in code block?): 

for i to 10 do iso[i] := implicitplot3d( f = i, x = -10 .. 10, y = -10 .. 10, z = -10 .. 10) end do

display(seq(iso[j], j = 1 .. 10))

Of course, visualizing the output is another issue.  I wish for an app to explore 3D data like Paraview. It does not have to be as sophisticated, but to display standard elements like isosurfaces and arbitrary cutting plane views would suffice.

If you know a package/app/procedure in Maple of this nature, please share it here. Thank you

I am trying to solve an equation by continuation method (one vairable and one parameter) containing more than 1000 terms and also hight exponents like 800. I want to get omega' values for different 'p' values.The maple file  is appended below. the link is maple_equation.mw

Many thanks in advance.

How do I find the constants in a solution containing whittaker functions?

the boundary conditions are:

c(x,0)=1

c(0,y)=1

[c'(1,y),x]=0

is there a built in command to get rid of terms of some order or greater in a given variable?

say i have:

myeq:=diff(mu[m](t), t, t) = (-2*d^2*phi[0]^2-2*phi[0]^2)*mu[m](t)/(2*C_J*phi[0]^2*L)+(J*L*sin(mu[m](t)/phi[0])*d^2+J*sin(mu[m](t)/phi[0])*L)*mu[p](t)^2/(2*C_J*phi[0]^2*L)+J*cos(mu[m](t)/phi[0])*d^3*mu[p](t)/(phi[0]*C_J)+(I_b*L*d*phi[0]^2-2*J*L*sin(mu[m](t)/phi[0])*phi[0]^2-d^2*Phi[xx]*phi[0]^2-phi[0]^2*Phi[xx])/(2*C_J*phi[0]^2*L);

and i want to get rid of all terms of O(d^2).

i know i can do this

lhs(myeq)=algsubs(d^2=0, rhs(myeq));

or

lhs(myeq)=convert(series(rhs(myeq), d=0, 2), polynom);

but I wonder if there is a built-in command. 

thanks

Sometimes Maple can be endlessly frustrating!!

A recent question of mine on how to distribute integrals into a sum was answered:

http://www.mapleprimes.com/questions/201500-Integral-Of-A-Sum--Sum-Of-Integrals-In-Maple

Now I need the same thing, except (as part of a solution) Maple chose to use Intat rather than Int.

So I modify the rule I have to look like this:

ToSumInt:=Intat(A::algebraic*Sum(v::algebraic,r::{name,equation})*B::algebraic,u::{name,equation})=Sum(Intat(A*v*B,u),r);

.. and for the heck of it I cannot get Maple to accept this, no matter what type I make u to be. The error message always looks the same:

Error, invalid input: Intat expects its 2nd argument, dx_at_t, to be of type symbol = algebraic, but received u::{name, equation}

and you may replace the u::{name,equation} with any other u::type you care to think of.

Any help would be appreciated.

M.D.

Edit: Also, useInt will not replace the Intat (in the expression I want to apply the rule to) with Int.

Sorry if this has been already posted.

When print() is invoked from a proc into a module, non-English characters are not properly displayed with Maple 18.

It works ok if it is invoked from within the workbook.

Example:
print("Están en perspectiva")

Put this sentence in a proc into a module and the character "á" wont be displayed

Output: "Est�n en perspectiva"

Any hint about how to treat this issue?

Thank you very much.

César Lozada 

To those who read this post.

I have a question about the property of local variables in a procedure.

Are names at the local variables in procedures, strings not names? 

If so, why can strings be evalated to have Integer, which is numbers for counting for continuing repetition? 

This is the code which searches for the maximum value. 

**************************************************

Max2 := proc()

local i, m :: string;
if nargs = 0 then
return -infinity
end if;
m:= args[1];
for i from 2 to nargs do
if args[i] > m then
m := args[i]
end if;
end do;
m;
end proc:
Max2(1,7,5,4,6,10,35,63.5,-10,5);

 ************************************************

type(a,name);

              this answer was true.

type(1, string)

              this answer was false, of course.

I limited the values of local variables to be string.

But, to the sequence 1,7,5,4,6,10,35,63.5,-10,5, any error messages wouldn't appear.

I hope any of you will give me an answer.

Thank you in advance.

taro yamada

So im trying to write a maple script that computes the Jordan form of a given (3x3)- matrix
A. If {a,b,c} is a basis with respect to which A is in Jordan form, then I'm trying to make it
plot the three lines spanned by a, b and c, in the standard coordinate system. I was hinted to use plot3d here.

sidenote: I know how to compute the jordan matrix of A, such by find the eigen vectors and generalised eigen vectors and putting them in as columns in a 3x3 matrix say S,   where S is invertible    then  (S^-1)*(A)*(S) = (J).

Thanks in advance. <3

Hi:

i follow the code newton raphson for solve system nonlinear ODE in maple,every body have code for it?

How would you do a quick sort or bubble sort?

When creating questions with Question Designer, how can I include notes that are visible to instructors only? I don't want students to see certain info. Is it in the "information fields"?