Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

can Maple solve bessel function by itself?

Or it seems i need to work out the solving process on myown?

my function is just like this, a simple one

please let me know if yo have any idea about it.

best regards

memdream

 

hello

as described in the title

what i want to define are these 

 

but i get this

as the result of expression of

I wonder what I should do to get this?

i.e define them at the same time.

if anybody knows, please help me.

best regards,

Memdream.

in my solution is 0.0000 . how can show it 0 and how can i show the solutions for example 10^-30 to zero(round to zero)

 

1.mw

 

Dear Maple users,

Is it possible to draw a portion of the figure on the figure it self (Inset) In maple? Below is an example pic:

 

Thanks,

 

how can I print my equation or configure maple so that they don't try to autofit themselves in the window? and what would the best way to export them to MATLAB be? I know there's the codegeneration command.

 

For example, I noticed that the output of some commands like StateSpace prints the equations in a way that if they are too long, they won't try tu fit, but rather it lets you scroll to the right:

 

I have the following inequalty

GAMMA*sigma^2 < 1+2*sigma[H]^2*p[th]-2*K*(1+sigma[H]^2*p[th])/N, assuming all of the variables are positive.

How could I make a 3d plot with three axes for Г, K and N for example?

i solve a dynamic problem but i get wrong answer. F (load) after 0.02 then to zero but d (displaciment) not go to zero 

1.mw  

New_Microsoft_Excel_Worksheet1.xlsx

A. how to find xx1,xx2,xx3,yy1,yy2,yy3 that
Determinant(Matrix([[xx1,yy1,1],[xx2,yy2,1],[xx3,yy3,1]])) =(1/2)*aa*d*s*u+(1/2)*aa*d*s*a*t+(1/2)*d*v*u*t+(1/4)*d*v*a*t^2;
 
B. how to find x1,x2,x3,x4,y1,y2,y3,y4 that expand(
(x2 - x1)*(y4 - y3) - (y2 - y1)*(x4 - x3)) = (1/2)*d*s*aa*v+(1/2)*d*aa*v*u*t+(1/4)*d*aa*v*a*t^2+(1/2)*aa*d*s*u+(1/2)*aa*d*s*a*t+(1/2)*d*u^2*t+(3/4)*d*u*a*t^2+(1/4)*d*a^2*t^3;

I have a set of 15 equations in state space form, and I'm just trying to program a way to substitute all variables representing the time derivates as for example x__s_dot (t) as diff(x__s (t), t)  and diff( x__s_dot (t),t) as diff(x__s (t), t,t) so I can manipulate my equations easier.

The last half of the list xx has the variables that I want to substitute, the x,y and u (with respective diff(x__s (t), t) ) variables are always constant, but in the case of zeta, eta and u I might have indexes from 1 to n , in this case just 1 and 2.

I tried a bunch of methods including creating vectors with the seq, op commands and substituting them into my original equations with subs, algsubs but none of the things I programmed worked the way I wanted, I ended up with a very inefficient way , doing it manually, but the point is that I wanted to automate this bit. I've been doing maple for about 2 weeks and I still struggle a lot with it, I would greatly appreciate any guidence/advice.


My original equations are actually much bigger and my original code too, I tried to shorten it a bit just for making my problem simpler to understand but the equations still look really big here, If there is another way to format this question please do advice so.

EDIT : I think my last code was hard to understand because my equations were too big and it looked all messy and horrible, I edited to show just one vector of variables, I think that if it works here, it should work on the rest of my code:


 

NULL

interface(rtablesize=50):

xx := [x__s(t), y__s(t), u__s(t), x__r(t), y__r(t), u__r(t), x__c(t), y__c(t), u__c(t), zeta__1(t), eta__1(t), u__1(t), zeta__2(t), eta__2(t), u__2(t), x__s_dot(t), y__s_dot(t), u__s_dot(t), x__r_dot(t), y__r_dot(t), u__r_dot(t), x__c_dot(t), y__c_dot(t), u__c_dot(t), zeta__1_dot(t), eta__1_dot(t), u__1_dot(t), zeta__2_dot(t), eta__2_dot(t), u__2_dot(t)]:

nq := 15;            
xxUNDOT := [ seq ( xx[i], i = 1..nq),
             seq ( diff(xx[i],t), i = 1..nq)];
 

15

 

[x__s(t), y__s(t), u__s(t), x__r(t), y__r(t), u__r(t), x__c(t), y__c(t), u__c(t), zeta__1(t), eta__1(t), u__1(t), zeta__2(t), eta__2(t), u__2(t), diff(x__s(t), t), diff(y__s(t), t), diff(u__s(t), t), diff(x__r(t), t), diff(y__r(t), t), diff(u__r(t), t), diff(x__c(t), t), diff(y__c(t), t), diff(u__c(t), t), diff(zeta__1(t), t), diff(eta__1(t), t), diff(u__1(t), t), diff(zeta__2(t), t), diff(eta__2(t), t), diff(u__2(t), t)]

(1)

DOTL := [ seq ( xx[i], i = 1..nops(xx)) ]:
DOTR := [ seq ( xxUNDOT[i], i = 1..nops(xx))]:
Vector[column](DOTL),Vector[column](DOTR);

for i from nq+1 to nops(xx) do
newxx    := subs[inplace][eval](
  op(i,DOTL) = op(i,DOTR) ,xx):
end do:

Vector(30, {(1) = x__s(t), (2) = y__s(t), (3) = u__s(t), (4) = x__r(t), (5) = y__r(t), (6) = u__r(t), (7) = x__c(t), (8) = y__c(t), (9) = u__c(t), (10) = `&zeta;__1`(t), (11) = `&eta;__1`(t), (12) = u__1(t), (13) = `&zeta;__2`(t), (14) = `&eta;__2`(t), (15) = u__2(t), (16) = x__s_dot(t), (17) = y__s_dot(t), (18) = u__s_dot(t), (19) = x__r_dot(t), (20) = y__r_dot(t), (21) = u__r_dot(t), (22) = x__c_dot(t), (23) = y__c_dot(t), (24) = u__c_dot(t), (25) = `&zeta;__1_dot`(t), (26) = `&eta;__1_dot`(t), (27) = u__1_dot(t), (28) = `&zeta;__2_dot`(t), (29) = `&eta;__2_dot`(t), (30) = u__2_dot(t)}), Vector(30, {(1) = x__s(t), (2) = y__s(t), (3) = u__s(t), (4) = x__r(t), (5) = y__r(t), (6) = u__r(t), (7) = x__c(t), (8) = y__c(t), (9) = u__c(t), (10) = `&zeta;__1`(t), (11) = `&eta;__1`(t), (12) = u__1(t), (13) = `&zeta;__2`(t), (14) = `&eta;__2`(t), (15) = u__2(t), (16) = diff(x__s(t), t), (17) = diff(y__s(t), t), (18) = diff(u__s(t), t), (19) = diff(x__r(t), t), (20) = diff(y__r(t), t), (21) = diff(u__r(t), t), (22) = diff(x__c(t), t), (23) = diff(y__c(t), t), (24) = diff(u__c(t), t), (25) = diff(`&zeta;__1`(t), t), (26) = diff(`&eta;__1`(t), t), (27) = diff(u__1(t), t), (28) = diff(`&zeta;__2`(t), t), (29) = diff(`&eta;__2`(t), t), (30) = diff(u__2(t), t)})

(2)

Vector[column](newxx);

Vector(30, {(1) = x__s(t), (2) = y__s(t), (3) = u__s(t), (4) = x__r(t), (5) = y__r(t), (6) = u__r(t), (7) = x__c(t), (8) = y__c(t), (9) = u__c(t), (10) = `&zeta;__1`(t), (11) = `&eta;__1`(t), (12) = u__1(t), (13) = `&zeta;__2`(t), (14) = `&eta;__2`(t), (15) = u__2(t), (16) = x__s_dot(t), (17) = y__s_dot(t), (18) = u__s_dot(t), (19) = x__r_dot(t), (20) = y__r_dot(t), (21) = u__r_dot(t), (22) = x__c_dot(t), (23) = y__c_dot(t), (24) = u__c_dot(t), (25) = `&zeta;__1_dot`(t), (26) = `&eta;__1_dot`(t), (27) = u__1_dot(t), (28) = `&zeta;__2_dot`(t), (29) = `&eta;__2_dot`(t), (30) = diff(u__2(t), t)})

(3)

``


 

Download subs.mw

I tried to make a for loop that substitutes each variable one by one, but it only seems to work for the last term, I don't know why, please help or suggest another method? I also have diff(x__s,t,t) terms in my equations, not just diff(x__s,t), but I guess only by specifying the first derivate it should work?

 

Thanks.

 

how can i solve coupled differential of equation with rkf or another method

maple.mw

A1 , A2, A3 ??

v=u+at                      (1)
s=u*t+1/2*a*t^2        (2)

below 3 equations, can substitute  (1)  into it to form (2)
s=1/2*(u+v)*t       (3)
v^2=u^2+2*a*s    (4)
s=v*t-1/2*a*t^2    (5)

can these 5 equations be considered as a solution set of solve function?

or

is only first 2 equations be a solution set?

if so, number of equations less than 5 variables, is there something missing?
 

i want to dsolve frome this equation . but i cant

1.mw


 

op(2, -(1/4)*z^2-(1/18*I)*z^3+(1/96)*z^4+(1/600*I)*z^5-(1/4320)*z^6-(1/35280*I)*z^7+(1/322560)*z^8+(1/3265920*I)*z^9-(1/36288000)*z^10-(1/439084800*I)*z^11+O(z^12))

-((1/18)*I)*z^3

(1)

op(2, series(Ei(I*z)-gamma-ln(z)-I*csgn(I*z)*Pi-I*z, z, 12))

2

(2)

series(Ei(I*z)-gamma-ln(z)-I*csgn(I*z)*Pi-I*z, z, 12)

series(-(1/4)*z^2-((1/18)*I)*z^3+(1/96)*z^4+((1/600)*I)*z^5-(1/4320)*z^6-((1/35280)*I)*z^7+(1/322560)*z^8+((1/3265920)*I)*z^9-(1/36288000)*z^10-((1/439084800)*I)*z^11+O(z^12),z,12)

(3)

is(series(Ei(I*z)-gamma-ln(z)-I*csgn(I*z)*Pi-I*z, z, 12) = -(1/4)*z^2-(1/18*I)*z^3+(1/96)*z^4+(1/600*I)*z^5-(1/4320)*z^6-(1/35280*I)*z^7+(1/322560)*z^8+(1/3265920*I)*z^9-(1/36288000)*z^10-(1/439084800*I)*z^11+O(z^12))

false

(4)

``


 

Download op_error.mw

Hello,

Is there a Maple command for copying a text to the clipboard ? (Apart from using the contextual menu or the Edit-Copy menu option)

For example:

ee := x^2+y^2+z^2

CopyToClipboard( convert(ee,string) )

Thank you,

César Lozada

 

 

 

4 5 6 7 8 9 10 Last Page 6 of 1286