Hello,

I am tring to do if–then–else in the do-loop for the parameter A[1]. Example: when k=3 we want to use A[1]=B[3] insted of B[1] but it seems my if–then–else did not do what I want. Any ideas?

Thank you

This is maple code:

> restart;
> k := 0; B[1] := 0.001/365: B[2] :=0.002/365:B[3] :=0.005/365:

> ode := diff(U(t), t) = -(A[1]+A[2]*U(t))*U(t);
> ic[0] := U(365*k) = 1000;

I have a 6 x 6 matrix. I want to employ Cramer's rule to solve Gx=Θ. I need to sequently substitute column Vector Θ into the matrix columns. I can't get this to work. I have only being using Maple for a couple of weeks so still finding my way.

myf:=(mu,sigma)->int(exp(-(x-mu)^2/2/sigma^2)/sigma/sqrt(2*Pi)/(1+exp(-x)),x=-infinity..infinity);

plot3d('myf(mu,sigma)',mu=-5..5,sigma=0.5..10);

This plot was never created. I think it's taking too long to evaluate it when it's divegent.

Hi all,

For my research purpose, I need to integrate this function. I understand that there maybe values for which this integration is divegent. But is there a way to look into it? Say creat...

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.