how to write a procedure to ceiling a number such as 3.5 to 4

HI everyone, im currently using maple 15 and im trying to animate at least the mass of my pendulum. My problem is that the  coordinate of each mass would be the solution of two numerical ODE and I can't animate this. I am able to plot the solution of my ODE (wich correpond to the two angle in the pendulum), but when it comes to animating... its an another piece of work. So here is what I tried: the code would correspond to the animation two masses... no rope yet.

I have a set of multi variate polynomial  equations. I need to setup a matrix [G][x]=[Θ]  Θ is a function of c₁ and c₂ up to 4th order terms

 My problem x=[c₁⁵, c₁⁴c_{2 ,} c₁³c₂² , c₁²c₂³ , c₁c₂⁴ , c₂⁵ ] these 6 5th order terms...

Hi,

I'm trying to plot a few phase portraits for this four equation system but maple wont let me...

> restart;

> s := 10; d[T] := 0.2e-1; k1 := 2.4*10^(-5); k2 := 3*10^(-3);
> N := 5; delt := .24; c := 2.4; T0 := 1; 

with(DEtools); 

with(plots);

[INSERT YOUR CODE HERE BY COPY-PASTING IN 1-D FORMAT]

Can someone please help me? Thanks

Hi everyone, im trying to solve the following system of ode (motion of a double pendulum) where im trying to plot numerical solution for T1(t) and T2(t). Been searching for a while... found nothing :) Here is my code , i also loaded the library plots, DEtools. . Im using maple 15

> Systeme := [2*(diff(T1(t), x, x))+cos(T1(t)-T2(t))*(diff(T2(t), x, x))+sin(T1(t)-T2(t))*(diff(T2(t), x))^2+19.6*sin(T1(t)) = 0, diff(T2(t), x, x...

Hi everyone,

I want to create a Gaussian PDF so I need to calculate Determinant(sigma) with sigma the covariance matrix of a gaussian variable.

If we call this variable alpha (which is a 12 dimension vector and represents the noise in a discrete dynamical equation), then sigma_ij=ExpectedValue(alpha_i*alpha_j)-ExpectedValue(alpha_i)*ExpectedValue(alpha_j)

and this is zero most of the time! So the covariance matrix is singular and the determinant is zero. 

Hello,

Could you please help me with the following system of equations:

> restart; with(student);
> Z[0] = 377;
> solve({Z = 1/(j*w*C[1]+1/(j*w*L[1])+1/R)+1/(j*w*C[2]+1/(j*w*L[2])), C[1] = w[c]*(w[0]^2-w[c]^2)*(w[2]^2-w[1]^2)/(2*Z[0]*(w[1]^2-w[c]^2)*(w[0]^2-w[1]^2)), C[2] = w[c]*(w[0]^2-w[c]^2)*(w[2]^2-w[1]^2)/(2*Z[0]*(w[1]^2-w[c]^2)*(w[2]^2-w[c]^2)*(w[2]^2-w[0]^2)), L[1] = 1/(w[1]^2*C[1]), L[2] = 1/(w[2]^2*C[2]), R[1] = subs(w = w[1], 2*Z[0]*(1-S[21...

I'm trying to write a program that solves sudoku's using a Groebner basis. I introduced 81 variables x₁ to x₈₁, this is a linearisation of the sudoku board.

The space of valid sudokus is defined by:

for i=1,…,81 : F_i=(x_i−1)(x_i−2)⋯(x_i−9) This represents the fact that all squares have integer values between 1 and 9.

for all x_i and x_j which...

circle_line_intersec.mw

i'm trying to write some code that will output an angle phi when a null event is reached.

see worksheet

Hello,

is there a way to activate the "undo" command after executing the entire worksheet? 

Before executing every change can be undone, but since I press the execute button 

often, I loose my chance to undo.


thanks

Does anyone know the release date of Maple 16?

I want to assign the output of a function that has been passed to timelimit. Currently I have to pass, say, X in by reference (as a single element list) and do  timelimit( 0.5, F(X,args) ); where F := proc( X::list ); X[1] := f(args); end proc; --- which is not really desireable. 

More generally I want my proceudre to skip loop iterations that are taking too long. 

Hi Everyone:

As far as I know maple has no way to define key bindings for it's actions. 

I would however like to be able to navigate around the worksheet using hjkl keys - namely alt+j would be equivalent to down_arrow, and so on. I find that I spend more time reaching for the arrow keys or the mouse than writing useful code (I might be exaggerating slightly here ;) ). I can implement this behavior easily using my window manager (FVWM) - meaning the window...

Dear all,

I am trying to simulate the penetration of the wavefunction into a potential barrier consisting of two barrier steps of different heights, V1 and V2.

My code looks like that:

> restart;
> m := 0.91e-30;h := 0.1055e-33; V1 := 10*0.16e-18; V2 := .8*V1;a := 0.10e-9;En := .6*V1; w := En/h;
                      ...

Hi experts,

I would like to ask what is the defult method's name for the numerical integration?

Thanks,