If i have to solve the eq. i have to guess what U2 is  until i get close to U1 that the error is ok. can i do this in maple. or du i have to make program.

Does the installation of Maple toolbox for matlab removes MuPad in matlab? 

Now when I type in mupad in matlab, it says undefined function.

Is there a way to keep both?

Hello Everyone,

I would like  to write a code to calculate the riemann sum approximation of the curve sin(sqrt(x^2+y^2)+1 

calculate the actual volume using integration

use the ranges x=[-2pi to 2pi]=y, and 20 subdivisions.

also display the curve and the parallelepiped approximations on the same plot.

 my code give me  an error..please someone can advise me which part of my code is not correct..

Thanks very much in advance!

> g := proc (x, y) options operator, arrow; sin(sqrt(x^2+y^2))+1 end proc;

 Z := g(x, y);

 X := Array(0 .. n); Y := Array(0 .. n); Z := Array(0 .. n);

c := -2*Pi; d := 2*Pi;

a := -2*Pi; b := 2*Pi;

n := 15;

s := 0;

X[0] := a; Y[0] := c; Z := g[X[0], Y[0]]; s := 0;

dX := (b-a)/N;dY := (d-c)/N;

    for i from 0 to N do      for j from o to N do     s:=s+Z[i]*dX*dY     X[i]:=X[i-1]+dX: Y[i]:=Y[i-1]+dY: Z[i]:=g([Xi],[Yi])  end do;

This the error message i am getting (Error, unterminated loop)

> 

Hi,

Recently a a simple problem which i can not handle by myself, made me confused.

I have simple code of maple which is not stable at all. Everytime I run the code, the final result which is the determinant of a matrix, changes and I can not see the problem with the code. In fact i noticed that problem occures when the matrix is being build by culculating the coefficients of some constant values.  I have attached the code. Could you see what is wrong here?

Thanks by the way.

Download Code.mw

code:

restart;
convert([92*x/790+18*x*(1-y)/1000=0.125,ln(46.59/x)=1/1.5*(ln(0.553/y)+2.5*ln((1-y)/(1-0.553)))],rational);
fsolve(%);
plots:-implicitplot(%%,x=0..2,y=0..1,numpoints=1000,thickness=2,gridlines,color=[blue,red]);

question:

there is something missing with the red line. I can't find the intersect point on the figure.

Is there something wrong with implicitplot? 

I purchased the student edition today, and after install and activation, it says the product is only available in 20 days. If I click the activate button, I repeat the successful activation process, but after another restart of Maple, the prompt comes out again, Eventually I clicked on OK, and I found that the license for me really expires in 20 days. 

I think it could be that my proof for student is not processed yet (no email), or a problem.

I wonder if there is anyone else purchased the student edition and encounter the same situation, and what is following next.

Hi everyone,

i'm trying to find out the euler angles from a rotation matrix.

I have a matrix that contains the result of multpilying the 3 axis rotation R(z,c)*R(y,b)*R(x,a) without knowing the values of the angles a,c,b (that's what I want to find out), so there are sinus and cosinus everywhere. There is another matrix containing the expected values that each equation in the first matrix will match.

My problem is that I eventually will change the order of the multiplication of the axis (i.e. R(x,a)*R(z,c)*R(y,b)) and I'm try to make maple compute this for me.

I defined R(x,a), R(y,b) and R(z,c) as follows:

Rx := Matrix (3,3, [1,0,0,0,cos(a),-sin(a),0,sin(a),cos(a)]);    
Ry := Matrix (3,3, [sin(b), 0, cos(b), 0, 1, 0, -sin(b), 0, cos(b)]);
Rz := Matrix (3,3, [cos(c), -sin(c), 0, sin(c), cos(c), 0, 0, 0, 1]);
and
RT := Multiply(Multiply(Rx,Rz),Ry);

Since here is alright. Now I want to match RT to the solution matrix.

g1 := RT[1,1] = mat[1,1];
g2 := RT[1,2] = mat[1,2];
g3 := RT[1,3] = mat[1,3];
g4 := RT[2,1] = mat[2,1];
g5 := RT[2,2] = mat[2,2];
g6 := RT[2,3] = mat[2,3];
g7 := RT[3,1] = mat[3,1];
g8 := RT[3,2] = mat[3,2];
g9 := RT[3,3] = mat[3,3];

soll := fsolve({g1,g2,g3,g4,g5,g6,g7,g8,g9},{a,b,c});

At this point I'm getting an error. I know that there are more equations than variables but the system is solvable anyway.
What maple trick can I do to solve this system or to find the good 3 equations?

Extra question:

Is there any method to match RT to mat without all those gX equations?

Thank you everyone.

Hi,

The default way a prompt shows on my maple is

[> blah blah blah

but i see in many worksheets including maple help pages that it shows as

> blah blah blah

My Question is how do i remove that left square bracket "[" from appearing ?

   Hi, there

 I work with maple 13 and made a file that its contet is plotting curves and its size is 

151 MB.Now when I want to open it the file cannot be opened completely only about a quarter of its contents rises.

Thanks for your help

REGARDS

Yegan

Hi,

I tried to write a procedure which performs a differentiation on list elements:

list1:= [x^3 , 2x , 3]
list2:= [x^5, 2x^4 , 2x]

dList:=proc(a::list, b::list)::list;
    local res::list, i::integer, n::integer; 
    description "add list elements and diff";
    n:=min(nops(a), nops(b)) ;
res:=[seq(0, i=1..n)];
for i from 1 by 1 to n do
     res[i]:=diff(a[i]+b[i],x)
end do;
res;
end proc:


dList(list1,list2)

However the Line:  res[i]:=diff(a[i]+b[i],x)

causes the Error: Error, illegal use of an object as a name;

Can someone give me a hint how to fix this?

Thank you very much.

I hope you be health and happy :-)

Please send me a pdf of Maple 17, I have many problem about Maple 17, and can't ask whole of them here!!!

restart:

Xsubs:=X=x*exp(epsilon*alpha[1]);

Ysubs:=Y=y*exp(epsilon*alpha[2]);

psubs := P = p*exp(epsilon*alpha[3])

where P(X,Y) and p(x,y). epsilon, alpha[i], i=1,2,3 are parameters.

how to changeover from p(x,y) to P(X,Y) in the following equation?

Eq:= diff(p,y$2)+diff(p,x)=0;

Cheers!

Hello,

 I need help writing a maple code to calculate the normal to the curve X^2-2X at x=1.5

and also display the curve, the tangent, and the normal on the same plot in the range x=[-10,10].

I manage to answer the question, but i am not quite certain if i have written a proper code..please if someone can review my solution and advise me i would really appreciate it....

I also would like to also make visisble the point of tangency at x=1.5 and give description of the figure but i am not sure how to do so...

how can i cut and paste the exact code with the images....it seesm that it doesnt allow me to do direct cut and paste so i changed it to 1-D and paste it into this page..

Thank you for your help.. 

restart;

    f := proc (x) options operator, arrow; x^2-2*x end proc;

   (D(f))(x);

 #Then, the slope of the line tangent to the graph of f(x) at the point with x-coordinate, x=1.5 is,

    D(f)*at*x = 1.5; 

   evalf((D(f))(1.5));

#The y-coordinate of the point is 

  f(x)*at*x = 1.5;

   evalf(f(1.5));

   showtangent(f(x), x = 1.5, color = [black, purple]);

   plot([f(x), (D(f))(1.5)*(x-1.5)+f(1.5)], x = -10 .. 10, y = -10 .. 10, color = [black, purple]);

  plot([x^2-2*x, x-1.5-.75, -x+.75], x = -10 .. 10, y = -10 .. 10, color = [black, purple, pink]);

hi, is there a way to change color of the page in Maple 18? in fact I am preparing my lectures using slideshow option and want to change the color instead of the default white.

Hi guys,

       I have the following Lagrange function. 

       (Note :w is the generalized co-cordinate.) 

       L := (diff(w(r, t), t, r))*v(r, t)

       ( I know

        If L =diff(w(t),t)*v(t)..... I can write EulerLagrange(L,t,w(t) but I dont know how to handle above function)

        Thanks in advance!!!