The Maple dsolver is very powerful, but everything has advantages and disadvantages. I was recently asked the following question.
Let us consider the system of ODEs
>restart; sys := [diff(y(x), x) = -(4*cos(x)*y(x)+z(x)*cos(x)^2+3*z(x))/(sin(x)*(cos(x)^2-9)),
>diff(z(x), x) = -(y(x)*cos(x)^2+3*y(x)+4*z(x)*cos(x))/(sin(x)*(cos(x)^2-9))]:
The functions
>y1 := C[1]*(cos(x)+1)^(1/2)/(cos(x)+3)^(1/2)+C[2]*(1-cos(x))^(1/2)/(3-cos(x))^(1/2):
>z1 := -C[1...

I've used the "unapply" command to calculate the several separarte transformation matrices (working a classical lamination theory problem in mechanics of composite laminates).  I'm now having a problem using these values to perform further calculations.  I've uploaded my file.  Basically, I;ve got to calculate some matrices and use those to calcuate other matrices.  Would anyone mind taking a look at it and letting me know if you have any suggestions...

I have a Matrix 'A' in maple with symbolic entries and I need to use that Matrix in matlab M file. The symbols with in 'A' are already initialized in M file.

How is this possible?

For Matrices 'A' of small dimension, I can use the codegeneration package to convert 'A' to Matlab Matrix and then copy paste.

But for matrices of 1000x1000, lets say, does not seem as a option.

Suppose one had 2 levels of evaluation and you only wanted to go back 1 level to uncover it's original components.  That would be no problem if nothing was cleared, but very unfortunate if it did. 

example,

a:=5*x^2-3*y^3:

b:=3*x^2+10*y^3:

c:=a+b
               

hi

for ex. i put the function X^2+2 on y. y:=x->x^2+2

how can i display what function is behind y later? of course i always could scroll up to the point where i defined it, but is there a shorter way? with variables, you enter them and see their value. it something similar possible with functions?

i now that i can see the function inmy example with y(x), but what if you have a function with 5 more?

thanks for help

I came across this post http://www.stiglerdiet.com/2012/02/10/mathematica-a-love-story/

and in the commented section near the end someone mentions this can not be elegantly done in Maple.  I would like to disprove this comment.  He says "There’s NO WAY any numeric language like R or Maple can try to do this."

I don't believe it, and I welcome anyone to disprove his...

`if`(p::`+`, nops, 1)(p);

'if'(p::'w', nops,1)(p);

w not exist

 I want to calculate successive roots of the following complex eq.

how can I set interval (?) for fsolve in the following code:

I'm pretty new to Maple 16.

I need to calculate a 6x6 transformation matrix for 8 different angles (8 different plys...classical lamination theory problem in mechanics of composites class), and I don't want to cut and paste the matrix 8 times changing the angle.  I want to put the 8 angles in a column vector (I think that's right...maybe it needs to be an array) and then have the transformation matrix call the angles one at a time to calculate a transformation matrix for each angle ...

Optimal_Portfolio_Ch.mw

I was using maple as a documentation software and was preparing some notes for my exam. I had a sudden computer shutdown and the file didn't saved properly.

The file now opens but most of the content is gone. Even the one which was properly saved.

I see Joe Riel 

Hello,

I need to solve bond orders in Maple the formula is

Bond order = [number of bonding electrons – number of anti-bonding electrons]/2

I need to implement it in Maple, I dk how

Thanks

when calculate the intersection of I2 and J2, i exactly follow this picture
https://skydrive.live.com/redir?resid=E0ED7271C68BE47C!311

however, it is G list which is not the stated in this picture

restart;
with(Groebner):
I2 := {x^2+y^3-1, x-y*x+3};
J2 := {y*x^2-1};
G1 := Basis({expand(w*I2[1]),expand(w*I2[2]), expand((1-w)*J2[1])}, 'tord', deglex(x,y));
G2 := Basis({expand(w*I2[1]),expand(w*I2[2]), expand((1-w)*J2[1])}, 'tord', deglex(w,x,y));

Dear All

How do I solve the followinh set of non-linear equations for the unknowns as given below:
    

    Zth1=R1(1-e^(-t1/tau1))+R2(1-e^(-t1/tau2))+R3(1-e^(-t1/tau3))+R4(1-e^(-t1/tau4))+R5(1-e^(-t1/tau5))+R6(1-e^(-t1/tau6))

    Zth2=R1(1-e^(-t2/tau1))+R2(1-e^(-t2/tau2))+R3(1-e^(-t2/tau3))+R4(1-e^(-t2/tau4))+R5(1-e^(-t2/tau5))+R6(1-e^(-t2/tau6))

    Zth3=R1(1-e^(-t3/tau1...

I'm trying to graph a contour map of this function: e^(y/x)

>with(plots)

>contourplot(e^(y/x),x=-3..3,y=-0.5..0.5)

I thought the contour map for this function would looks more like several lines intersects at the origin. But when I graph it on Maple, it doesn't looks like the graph I'm looking for.

Did I put in a wrong range?

Hi everyone, I'm new to the forum and also new to Maple.

I am trying to do the following:

"Write a procedure, marks, to take a student’s mark and output the appropriate grade, based on the following table:"