MaplePrimes Questions

Hi
I am confused by solve and SolveEqutions and their options!
Solving same equation results different solutions.
Which  solution is correct and why?
How i can be sure one of the solutions is correct about the other equations?
I attached my code.

Thanks.

Test2_MaplePrimes.mw

Its difficult for me to explain what I'm looking for so I will try with an example:

Example 1

Is there anyway any function that will say for example the sum of squares of integers for a given n say n=4, 1^2 + 2^2 + 3^2 + 4^2 = 30 that will give the equation ((n^3)/3) + ((n^2)/2) + (n/6) that gives the sum for a given n?

It would be determined from the series 1^2 + 2^2 + 3^2 + 4^2 + ... or the equation for the elements 1^2, 2^2, 3^2, 4^2 the equation that gives the sum; in this case the well known ((n^3)/3) + ((n^2)/2) + (n/6) .

 

So far I've got FindSequenceFunction[{1, 4, 9, 16}, n] but this only finds the individual elements not an equation for the sum of them. Of course I already know the equation for producing the individual elements.

How can plot  f := proc (x, y) options operator, arrow; sin(x)*cos(y) end proc  as surface and contour (contour is project on surface)

Hi

I want to solve quadratic eqution involving more than 2 parameters...want to analize unique soltuions and real roots also want to plot the real and unique region on graphs ....thanx

the equation is -delta*(Q*S*alpha*b-a*alpha^2+M*c+b*delta) where alpha is the varible and rest are paramters 

 

Hei!

In my Maple project I want to put a simple sin graph on top of an image. But how do I create a background image with a graph on top?


 

restart

PI := proc (p, e, q) options operator, arrow; p*S(p, e, q)+v*(q-S(p, e, q))-w*q-g(e)+(w-c)*q end proc

proc (p, e, q) options operator, arrow; p*S(p, e, q)+v*(q-S(p, e, q))-w*q-g(e)+(w-c)*q end proc

(1)

S := proc (p, e, q) options operator, arrow; q-(int(F(x), x = 0 .. q)) end proc

proc (p, e, q) options operator, arrow; q-(int(F(x), x = 0 .. q)) end proc

(2)

NULL

PI(p, e, q)

p*(q-(int(F(x), x = 0 .. q)))+v*(int(F(x), x = 0 .. q))-w*q-g(e)+(w-c)*q

(3)

g := proc (e) options operator, arrow; (1/2)*mu*e^2 end proc

proc (e) options operator, arrow; (1/2)*mu*e^2 end proc

(4)

``

F := proc (x) options operator, arrow; int(f(x), x = 0 .. q) end proc

proc (x) options operator, arrow; int(f(x), x = 0 .. q) end proc

(5)

P := proc (p, e, q) options operator, arrow; simplify(eval(PI(p, e, q), [alpha = 50, w = 10, mu = 10, c = 5, v = 1, f(x) = 1/2])) end proc

proc (p, e, q) options operator, arrow; simplify(eval(PI(p, e, q), [alpha = 50, w = 10, mu = 10, c = 5, v = 1, f(x) = 1/2])) end proc

(6)

Diff_p := diff(P(p, e, q), p)

q-(1/2)*q^2

(7)

Diff_e := diff(P(p, e, q), e)

-10*e

(8)

Diff_z := diff(P(p, e, q), q)

-p*q+p+q-5

(9)

``

``


 

Download Profit_code.mw

 

 

 

why i can not evaluate 29 polynomial. maple try to evaluate last 7hr, how many time required too solve it?

 

How to change the numbers shown in the diagram shapes from y=0.5- to rigth form y=-0.5a.mw

Download a.mw


 

 

 

 

 

 

How to calculate the following for a=1,b=3

where  is the Wigner 3j symbol.

 

 

 

Top part of worksheet gives me a and b; lower part does not

Please tell me where I go wrong. Thanks

Test1.mw


 

restart

p := (alpha+c*(beta-k^2/mu)+z-Lambda(z))/(2*beta-k^2/mu)

(alpha+c*(beta-k^2/mu)+z-Lambda(z))/(2*beta-k^2/mu)

(1)

x := simplify((p-c)/(p-v))

mu*(beta*c+Lambda(z)-alpha-z)/(Lambda(z)*mu+((-c+2*v)*beta-z-alpha)*mu+k^2*(c-v))

(2)

``

xx := p-c-(p-v)*x

(alpha+c*(beta-k^2/mu)+z-Lambda(z))/(2*beta-k^2/mu)-c-((alpha+c*(beta-k^2/mu)+z-Lambda(z))/(2*beta-k^2/mu)-v)*mu*(beta*c+Lambda(z)-alpha-z)/(Lambda(z)*mu+((-c+2*v)*beta-z-alpha)*mu+k^2*(c-v))

(3)

Lambda := proc (z) options operator, arrow; int((z-u)*phi(u), u = 0 .. z) end proc

proc (z) options operator, arrow; int((z-u)*phi(u), u = 0 .. z) end proc

(4)

``yy := eval(xx, [alpha = 50, mu = 10, c = 5, v = 1, beta = 2, k = 2, phi(u) = 1/2])

100/9+(5/18)*z-(5/72)*z^2-10*(136/9+(5/18)*z-(5/72)*z^2)*(-40+(1/4)*z^2-z)/((5/2)*z^2-544-10*z)

(5)

``

zvalue := solve(yy, z)

z

(6)

yy

100/9+(5/18)*z-(5/72)*z^2-10*(136/9+(5/18)*z-(5/72)*z^2)*(-40+(1/4)*z^2-z)/((5/2)*z^2-544-10*z)

(7)

``


 

Download Dummy_Code.mw

Hello,

Please , how can we plot two figures with different X-axis and the same Y-axis ?

Thank you

 


 

restart

``

Lambda := proc (z) options operator, arrow; int((z-u)*phi(u), u = 0 .. z) end proc

proc (z) options operator, arrow; int((z-u)*phi(u), u = 0 .. z) end proc

(1)

k := diff(Lambda(z), z)

int(phi(u), u = 0 .. z)

(2)

k

int(phi(u), u = 0 .. z)

(3)

m := subs(phi(z) = 1/2, k)

int(phi(u), u = 0 .. z)

(4)

``

``


 

Download dummy2.mw

Hi,

 

I'm trying to learn the dsolve command on Maple. I have a second order differential equation and two initial conditions to solve it.

I'm following maple help's procedure to solve it. But in the end i don't get my constants solved (which i do if i use Matlab i.e.)

What am i doing wrong?

restart;
ode := diff(Y(y), y, y) = 3*(1-(y/b)^2)*q/(4*k*b);
                                    /     2\  
                                    |    y |  
                                  3 |1 - --| q
                                    |     2|  
                  d  / d      \     \    b /  
                 --- |--- Y(y)| = ------------
                  dy \ dy     /      4 k b    
dsolve(ode);
                         2       4                
                    3 q y     q y                 
             Y(y) = ------ - ------- + _C1 y + _C2
                    8 b k        3                
                             16 b  k              
ics := Y(b) = Ts-Tb; (D(Y))(-b) = 0;
                         Y(b) = Ts - Tb
                          D(Y)(-b) = 0
dsolve({ics, ode});
                      2       4          
                 3 q y     q y           
          Y(y) = ------ - ------- + _C1 y
                 8 b k        3          
                          16 b  k        

               16 _C1 b k + 16 Tb k - 16 Ts k + 5 b q
             - --------------------------------------
                                16 k                 

or a screenshot:

 

 

 

 

I have a package that uses fsolve, a Maple function known (to me) to often fail when UseHardwareFloats is set to true. Since I set UseHardwareFloats:=true in my .mapleinit (for I want hardware performance when using floats) I set UseHardwareFloats:=deduced; in the ModuleLoad routine of this package. ModuleLoad executes as evidenced by a print statement in the routine that does in fact print, and within ModuleLoad, UseHardwareFloats is set to deduced. But at the worksheet level, UseHardwareFloats remains set to true.
 The only way I can set it is to set it explicitly in the sheet that calls the package. I tried various things, like putting it into the body of the module, using :-UseHardwareFloats, and so on. Nothing seems to work. While not a fatal issue it is bothersome & has had me run around with code not working several times.

Any hint would be appreciated.

Mac Dude

 

First 741 742 743 744 745 746 747 Last Page 743 of 2426