MaplePrimes Questions

Pls what is the maple code for the fourth order Runge-kutta method for five systems of ODE of several variables?

I want to be able to make arrows on my axis to indicate positive direction, since it is a requirement at my school. I would like to do it with the "plots" function.

Here is what I wish I had in maple http://imgur.com/a/OXi5t

Hi Mapleprimes,,,


We see Maple file.

 

We are aware of oeis.org/A103435/


sequence data

0 2
1 4
2 16
3 48
4 160
5 512
6 1664
7 5376
8 17408
9 56320
10 182272
11 589824
12 1908736
13 6176768
14 19988480
15 64684032
16 209321984
b103435

 

Where F(m)=2*F(m-1)+4*F(m-2).


Goal - to increase refereed database.

 

M.mw

Regards,

Matthew Charles Anderson

Legislative Assistant

Oregon USA

Age 42.

5*x-7 = 3*x+2;
                             "(->)"
lhs(5*x-7 = 3*x+2) - rhs(5*x-7 = 3*x+2) = 0;
                          2 x - 9 = 0
                              [x]
isolate( 2*x-9 = 0, %ARG1 );
Error, (in isolate) 2*x-9 = 0 does not contain %ARG1
plz explain me what is this error?

If input a group of outputs which are binary numbers

can it simplify to give a logic that can output these outputs

If input ["1","1","2","2","2","77"]

output a graph such as

1 has 2 times

2 has 3 times

77 has 1 time

3.mw

There are 2 questions actually. The first as the title says is about taking the integral. I have 2 functions that were found numerically from the system of differential equations (see the file), and I need to take the integral of the expression that includes both of them. Maple gives me something like Int () = Int (), so it doesn't solve anything. Why can it be?

The second question is about varying the boundary conditions. If, for example, I have the system with the condition like R(x_0)=R_0, can I get the plot  of R(x,R_0)? In my case I need to vary conditions on R and mu (R(0)=R_0 and mu(0)=mu_0) and then get the plot of the integral in relation to R_0 and mu_0. Is it even possible?

Hello,

I have two simple module's:

1) ------------------------------------------------------------------------

Point := proc(xx::float,yy::float)
    return module()
        local x := xx,y := yy;
        export ShowPoint,GetX,GetY;
        
        ShowPoint := proc()
            printf("Point X,Y -> [%f,%f]",x,y);
        end proc;

        GetX := proc()
            return x;
        end proc;        
        GetY := proc()
            return y;
        end proc;    
    end module:
end proc:

2) ------------------------------------------------------------------------

PointMath := module()
    option package;  
    export `+`;
    
    `+` := proc(a::Point(float,float),b::Point(float,float))
        option overload;
        Point(a:-GetX()+b:-GetX(),a:-GetY()+b:-GetY());
    end proc;
end module:

------------------------------------------------------------------------

Next I use first module:

p1:=Point(1.2,1.4);

p2:=Point(1.0,2.0);

Finally I want to add above two points:

with(PointMath)

p2:=p1+p2

The results is:

p2:=p1+p2

Why is not called operator '+' and two points are not added?

Best,

Rariusz

 

 

 

Hi

How can I remove a subset containing a number from a partition?

eg S := [seq(2..5)]

P:=combinat:-setpartition(S, 2);

[[[2, 3], [4, 5]], [[2, 4], [3, 5]], [[2, 5], [3, 4]]]

I want to remove any partition involving 5

leaving

[[[2, 3]], [[2, 4]], [[3, 4]]]

i tried

remove(has,P,5)

 

 

Can someone answer why i get Float(undefined) while computing integral

Here is my example

NULL

restart

with(plots):

NULL

alpha0 := .3;

.3

(1)

R := .1;

.1

(2)

sigma := 636;

636

(3)

Kfc := 101;

101

(4)

Ksc := 9;

9

(5)

E := 2*10^5;

200000

(6)

`δth` := 2*10^(-7);

1/5000000

(7)

`δc` := 0.8e-1*10^(-3);

0.8000000000e-4

(8)

p := 400;

400

(9)

`lδz` := 12.1*10^(-3);

0.1210000000e-1

(10)

lKz := 20.3*10^(-3);

0.2030000000e-1

(11)

eta := 10^(-5);

1/100000

(12)

l0 := 0;

0

(13)

xi := p/sigma;

100/159

(14)

KImax := evalf(p*sqrt(Pi*l));

708.9815404*l^(1/2)

(15)

``

`Nδ` := evalf(E*sigma*(int((-xi^2+1)*(-KImax^2+Kfc^2)/((KImax^2-Ksc^2)*((-R^4+1)*(KImax^2+Ksc^2)+eta*E*sigma)), l = l0 .. `lδz`))/alpha0);

Float(undefined)

(16)

``

NULL

NULL

NULL

NULL

NULL

``

``

``

``

``

NULL

NULL


 

Download N.mw

Thanks in advance

"01000001" is A

how convert string of binary into a characters?

and

how to convert this string into 65 too?

Dear All,

restart;

with(QDifferenceEquations);

z := .8; q := .9;
                              0.8
                              0.9
QGAMMA(z, q);
Error, (in QDifferenceEquations:-QGAMMA) wrong type of arguments
How do we get the numerical result from Maple? It just only do algebraic calculation. 

Thanks

 

Primetable :=table[(2=b,3=c)]

so far, hard code in above code

is there a generic code for general

according alphabet order 5=d,7=e etc

I want to use maple2015 to perform a Monte-Carlo simulation of a one dimensional Harmonic Oscillator.

Consider a particle having displacements, x=na where n is an integer, i.e 0 plus minus one, etc.

The potential energy of the particle is given by U(x)= 1/2 kx^2 = 1/2 ka^2 n^2.

For this simulation will choose ka^2 = 0.01 k_B T where k_B is Boltzmann's constant, and T is the temprature in Kelvin.

So we should perform a Monte-Carlo simulation starting from x(0)=20a, and perfroming individual steps by selecting an attmept to move x=+-a ,with equal probability, and choose whether to make the actual step according to energy change expected in that step.

I need to draw particle's position and x(t), and particle's energy E(t).

After that I need to calculate and draw the autocorrelation functions \phi_x(t) and \phi_E(t).

I found some code that might help me, but I am nor sure how to implement the above two lines on maple; your help is appreciated.

http://pubs.acs.org/doi/suppl/10.1021/ed100414p/suppl_file/ed100414p_si_001.pdf

 

Op(1,ifactor(12))

is 2^2

it has power

how to get 2 only

First 985 986 987 988 989 990 991 Last Page 987 of 2434