Maple 18 Questions and Posts

These are Posts and Questions associated with the product, Maple 18

I encountered this bizarre inconsistency issue that Maple18 generates different outputs when executing the same command:

test_res2:= factor( simplify( expand( value( subs( Perturbation_Sol, EQ_PX2_order_7 ) ) ) ) )

'EQ_PX2_order_7' is a rational expression in sin(i0), cos(i0), sin(uL), and cos(uL) with rational coefficient terms. It also has inert differentiation terms Diff( * , uL ).

'Perturbation_Sol' is a set of 171 elements in the form of 'parameter_name = expression'.

My goal is to check if substituting 'Perturbation_Sol' into 'EQ_PX2_order_7' yields 0. Since 'EQ_PX2_order_7' has inert differentiation terms, I've applied 'value' after using 'subs'. Then I apply 'expand', 'simplify', and 'factor' to reduce the result to the simplest form.

However, Maple18 generates different outputs when I just execute this repeatedly. Please see the worksheet "test.mw" for details. Any insight will be greatly appreciated! Also, I wonder if the same issue would happen when the worksheet is executed with newer versions of Maple.

EQN_SOL_test1.mla

test.mw

 

 

I encountered the problem with .m files originally. But MaplePrimes doesn't allow uploading .m files, so I had to save the expressions into the file "EQN_SOL_test1.mla", which is included in this question. Below we load the expressions from the .mla file first, and then save them into a .m file in order to recreate the problem that I encountered.

restart;

>

 

read "EQN_SOL_test1.mla":

# Load 'EQ_PX2_order_7' and 'Perturbation_Sol'

 

save

EQ_PX2_order_7,
Perturbation_Sol,

"EQN_SOL_test1_m.m";

# Save the expressions into a .m file

 

Now we demonstrate the inconsistency problem with .m files. Notice that Maple generates 3 possible outputs:

test_res2 := 0

test_res2 := -(1/4)*rho0^2*a0^2*Be^2*cos(uL)*J2re*R_earth^2*(5*cos(i0)^2*cos(uL)^2-7*cos(i0)^2-5*cos(uL)^2+4)/sha

 

test_res2 := -(1/8)*rho0^2*a0^2*Be^2*cos(uL)*J2re*R_earth^2*(5*cos(i0)^2*cos(uL)^2-7*cos(i0)^2-5*cos(uL)^2+4)/sha

 

The last 2 outputs cannot be reduced to 0 since 5*cos(i0)^2*cos(uL)^2-7*cos(i0)^2-5*cos(uL)^2+4 is nonzero as shown below.

 

 

plot3d( 5*cos(i0)^2*cos(uL)^2-7*cos(i0)^2-5*cos(uL)^2+4 , uL=0..2*Pi, i0=0..2*Pi );

 
 

restart;

 

read "EQN_SOL_test1_m.m":

 

length( EQ_PX2_order_7 );

939346

(1)

length( Perturbation_Sol );

2082306

(2)

numelems( Perturbation_Sol );

171

(3)

Perturbation_Sol[1..5];

# Just to give an example of what the elements in 'Perturbation_Sol' look like

{PX1[1] = 0, PX1[2] = 0, PX1[3] = -(1/4)*rho0*a0*Be, PX1[4] = (1/2)*rho0*a0*Be*WEra*cos(i0)-(3/16)*R_earth^2*a0*rho0*(3*cos(i0)^2-1)*J2re*Be/sha+(1/4)*Be*a0*rho0*X10[3]/sha, PX1[5] = (1/4)*rho0*a0*X10[4]*Be/sha-(1/256)*R_earth^4*a0*rho0*(163*cos(i0)^4-110*cos(i0)^2+19)*J2re^2*Be/sha^2+(3/16)*R_earth^2*a0*rho0*(3*cos(i0)^2-1)*J2re*Be*X10[3]/sha^2+(3/8)*cos(i0)*R_earth^2*WEra*a0*rho0*(3*cos(i0)^2-1)*J2re*Be/sha-(1/48)*Be^3*a0^3*rho0^3*s1/sha^2-(1/8)*Be*a0*rho0*X10[3]^2/sha^2-(1/2)*cos(i0)*WEra*a0*rho0*Be*X10[3]/sha-(1/16)*rho0*a0*(3*cos(i0)^2+1)*Be*WEra^2-(1/32)*Be^2*J2re*R_earth^2*a0^2*rho0^2*sin(i0)^2*sin(2*uL)/sha^2}

(4)

 

 

for j from 1 to 50 do
    test_res2:= factor( simplify( expand( value( subs( Perturbation_Sol, EQ_PX2_order_7 ) ) ) ) );
end do;

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

 

0

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

 

0

 

0

 

0

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

 

0

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

 

0

 

0

 

0

 

0

 

0

 

0

 

0

 

0

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

 

0

 

-(1/4)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

0

(5)

 

 

It seems that with .mla files the problem occurs in a different way! With "EQN_SOL_test1.mla", the outputs for all 50 iterations stay the same as

test_res2 := -(1/4)*rho0^2*a0^2*Be^2*cos(uL)*J2re*R_earth^2*(5*cos(i0)^2*cos(uL)^2-7*cos(i0)^2-5*cos(uL)^2+4)/sha               (A)

 

but they may all change to the following different result after retarting many times:

test_res2 := -(1/8)*rho0^2*a0^2*Be^2*cos(uL)*J2re*R_earth^2*(5*cos(i0)^2*cos(uL)^2-7*cos(i0)^2-5*cos(uL)^2+4)/sha               (B)

 

In particular, after a large number of test runs (i.e., open the file "test.mw", execute the worksheet, close the file, and repeat), the result (B) has only occured twice. The second appearance is saved here for you to view. Once you re-execute this worksheet, most likely all outputs below will change back to (A), and (B) will only reappear after a large number of reruns.

 

restart;

 

read "EQN_SOL_test1.mla":

# Load 'EQ_PX2_order_7' and 'Perturbation_Sol'

 

 

for j from 1 to 50 do
    test_res2:= factor( simplify( expand( value( subs( Perturbation_Sol, EQ_PX2_order_7 ) ) ) ) );
end do;

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

 

-(1/8)*rho0^2*a0^2*Be^2*cos(uL)*R_earth^2*J2re*(5*cos(uL)^2*cos(i0)^2-5*cos(uL)^2-7*cos(i0)^2+4)/sha

(6)

 

Download test.mw

scmch.mw

I can't get a graph. Is this code is correct.Please help.

HTR.mw

In above problem, Additionally How to  plot  heat transfer rate  Q versus L^2  for distinct porosity parmeters(Sh) , using  heat transfer rate formula, Q = (q*L)/(k*A*T[b])=theta'(1).

using  [Sh = 0.1, L^2 = 0.1, Nr =0 .1, Ha =0 .1, Pe = 0.1],  [Sh = 0.3, L^2 = 0.3, Nr = 0.1, Ha = 0.1, Pe =0 .1],   [Sh = 0.5, L^2 =0 .5, Nr =0 .1, Ha = 0.1, Pe =0 .1].

How to solve and plot a ODE system in RK method.
eq1 := diff(f(x), x, x, x)-(1/2)*Sc*sin(alpha)*g(x)*(diff(g(x), x, x))+(1/2)*x*cos(alpha)*(diff(f(x), x, x))+(1/2)*sin(alpha)*f(x)*(diff(f(x), x, x)) = 0; eq2 := (diff(g(x), x, x, x))/Pm+(1/2)*x*cos(alpha)*(diff(g(x), x, x))+sin(alpha)*f(x)*(diff(g(x), x, x))-sin(alpha)*(diff(f(x), x, x))*g(x) = 0; eq3 := (diff(theta(x), x, x))/Pr+(1/2)*x*cos(alpha)*(diff(theta(x), x))+(1/2)*x*(diff(f(x), x))*(diff(theta(x), x))+sin(alpha)*(x*(diff(f(x), x))-f(x))*(diff(theta(x), x))-Nb*(diff(s(x), x))*(diff(theta(x), x))-Nt*(diff(theta(x), x))^2+(1/4)*Sc*Br*sin(alpha)^2*(diff(f(x), x))^2*(x*(diff(g(x), x))-g(x))+(diff(g(x), x))^2*(x*(diff(f(x), x))-f(x)) = 0; eq4 := diff(s(x), x, x)+S*((1/2)*cos(alpha)*x*(diff(s(x), x))+(1/2)*sin(alpha)*f(x)*(diff(s(x), x)))+Nt*(diff(theta(x), x, x))/Nb = 0

ics := f(0) = 0, (D(f))(0) = 1, g(0) = 0, (D(g))(0) = 1, theta(0) = 1, s(0) = 1; bcs := (D(f))(100) = 0, (D(g))(100) = 0, theta(100) = 0, s(100) = 0

alpha = - 30 degree, Sc = 1.0, Pm = .1, Pr = 6.2, Nb = .1, Nt = .1, Br = .5, S = 1

Good evening,could you please help me to find the values of Phi[2],Phi[3],Phi[4],....

The right hand side of the for loop is the equation i need.Is this maple code correct.

rps1.mw

Certainly a standard question.

I have an integer n*n matrix A (the entries are explicitly integers; there is no variable -type x- in the matrix). I want the Smith normal form of A, that is A=UDV where U,V are integer matrices with determinant +-1 and D is a diagonal matrix with -eventually- some zero and positive integers d_i s.t. d_i divides d_{i+1}.

"SmithForm()" doesn't work directly (I get rational -non integer- matrices). Maybe it is necessary to declare the matrix A as 'Matrix(integer)' ...
Thank you in advance for your help.

How do  I solve system of differential equations in finite difference method or finite element method?

eq1 := (diff(f(x), x, x, x))*(a*beta*f(x)^2-1)+(diff(f(x), x))^2-2*a*beta*f(x)*(diff(f(x), x))*(diff(f(x), x, x))+(diff(f(x), x))*(M+k[1])-(diff(f(x), x, x))*f(x)-(alpha*theta(x)+delta*phi(x))/rho = 0;

eq2 := -(diff(theta(x), x, x))*K[SB]*(Df-(Rd+k[hnf]/k[bf])/Pr)+N[t]*K[SB]*(diff(theta(x), x))^2-N[b]*(diff(theta(x), x))*(diff(phi(x), x))-(diff(f(x), x))*(diff(theta(x), x))-lambda*theta(x)-mu*Ec*(M*(diff(f(x), x))^2+(diff(f(x), x, x))^2) = 0;

eq3 := diff(phi(x), x, x)+Le*Sr*(diff(theta(x), x, x))+Le*f(x)*(diff(phi(x), x)) = 0;

ics := f(0) = 0, (D(f))(0) = 0, theta(0) = 1, phi(0) = 1;

bcs := (D(f))(100) = 0, theta(100) = 0, phi(100) = 0;


Parameters1 := rho = 2063.905, k[hnf] = .29942, k[bf] = .2520, mu = .38694, a = .1, beta = 5, k[1] = 2.0, M = 10, alpha = 20, delta = 20, K[SB] = .5, Df = 3, Pr = 1.2, Rd = 5, N[t] = 1.2, N[b] = 1.0, lambda = 1.5, Ec = 5, Le = .1, Sr = .1;

 

Hello everyone,

I'm trying to learn how to use Shoot Library 9.

Unfortunately, I'm not doing very well. I'm getting an error that I don't understand. I don't know where it comes from.

Please help me solve this odes system using this library.

I attached my Maple worksheet file.

ShootLib_Test.mw

Please help to plot this equation 

fwf-v.mw

How to plot this equation

 y(x):=

where,

A := Matrix([[1, -1, 1, -1], [1, 1, -1, -1], [-1, 1, 1, -1], [1, 1, 1, 1]]);
B := Matrix([[1], [0], [1], [0]])

Could you please help me to solve this error.

my code is here.

CR-C.mw

How to conver a patial differetial equation to ordinary differential equation with or without dchange?
 

restart

declare(u(x, y, t), v(x, y, t), T(x, y, t), C(x, y, t), eta(x, y, t), psi(x, y, t), f(eta), theta(eta), phi(eta));

declare(u(x, y, t), v(x, y, t), T(x, y, t), C(x, y, t), eta(x, y, t), psi(x, y, t), f(eta), theta(eta), phi(eta))

(1)

eta := proc (x, y, t) options operator, arrow; y/(nu*t+nu*x/U[w])^(1/2) end proc:

eq1 := diff(T(x, y, t), t)+u*(diff(T(x, y, t), x))+v*(diff(T(x, y, t), y))-sigma*(diff(T(x, y, t), y, y))-epsilon*D[B]*(diff(T(x, y, t), y))*(diff(C(x, y, t), y)) = 0

diff(T(x, y, t), t)+U[w]*(D(f))(y/(nu*t+nu*x/U[w])^(1/2))*(diff(T(x, y, t), x))+(-(1/2)*f(y/(nu*t+nu*x/U[w])^(1/2))*nu/(nu*t+nu*x/U[w])^(1/2)+(1/2)*(D(f))(y/(nu*t+nu*x/U[w])^(1/2))*y*nu/(nu*t+nu*x/U[w]))*(diff(T(x, y, t), y))-sigma*(diff(diff(T(x, y, t), y), y))-epsilon*D[B]*(diff(T(x, y, t), y))*(diff(C(x, y, t), y)) = 0

(2)

``


 

Download pde_to_ode.mw

Can you change f(eta) to upflow curve and theta(eta) to downflow curve.

In my Problem,Boundary Conditions are

theta(infinity) = 0, (D(f))(infinity) = 1 , (Take, eta =infinity)

Flows will be correct for what value is taken for infinity .

I take  eta = 5. and also tried changing ranges  but could't find it.Please Help to fix the curve.

my code is,

SM.mw

How to find series values.I got this error.Please Help.

Maple code for the problem is

TFBE.mw

hallo every body 

Please, I have a small problem in this program in Maple 18 at the end of the programme when i want to calculat the intgral given in variable (f2i)  i see this message and I can't calculate the integral 

Please what is the solution to calculate this integral?

The program in Maple 18 is below 

Thanks a lot 

hopf12.mw

1 2 3 4 5 6 7 Last Page 1 of 88