Hey guys..

I need some help on determining the correlation coefficient between my data points for an Influenza epidemic and the SIR model, with my estimated parameters Beta and Alpha. The initial conditions for the equations are defined at the beginning of my worksheet:

SIR_model.mw

Hope someone can help me ;)

How to successfully deploy Matlab function [B, IX] = sort (A, ...)
IX is a permutation vector of the corresponding column of A

M := Matrix(6, 3, [2, 3, 3, 5, 7, 8, 12, 5, 9, -3, 4.1, 7, 7, 7, -3, 9, 3, 8]):
M;
L := <(seq(op(convert(sort(M[() .. (), j]), list)), j = 1 .. op(M)[2]))>;
L := ArrayTools[Alias](L, [op(M)[1 .. 2]], Fortran_order);
P := Matrix(op(M)[1 .. 2],  (i, j)-> ListTools[Search](L[i, j], M[() .. (), j]));
# from Matlab

I had MAPLE 12 on WINDOWS VISTA machine and the harddrive crashed.  Now I have a newer WINDOWS 7 machine & would lke to put my individually licensed MAPLE 12 on this machine.  I do not have a DVD because the original installation occurred via MAPLESOFT website if I recall correctly.

Can anyone direct me as to what I need to do?

THANKS

Hi,

I'm having some problems with the generation of the typical stiffness matrices seen in FEA formulations. I have got an elementary matrix that is looped through the entire stiffness matrix to generate it. As there is juxtaposition between the elementary matrices, there is some recursive assignment involved that I can't handle properly. Here is the code for a 6x6 elementary matrix with a 3x3 juxtaposition. n is the nombre of recursions needed to fill the stiffness matrix. M is the elementary matrix

On my keyboard, Cntrl = action for evaluating an expression inline does not work. The reason is that = is the upper level action (with the shift key) on a key with 0 at the lower level. Doing Cntrl Shift = does not work either. How do I get this done on my keyboard?

What is the number of all the solutions of the equation frac(x*floor(x)) = 1/2 belonging to RealRange(1,100)? How to count it with Maple?

Hi,

is the series expansion for the following expression bugged in maple or am I missing some crucial thing of series?

series((1/epsilon-log(4)),epsilon=0);

results in epsilon^-1 -2 -2*ln(2)

for xmaple 16.01. The console version does not add the wrong term "-2".

Apparently it works here: 

How is it possible to let multiple cores compute a loop faster?

For exaple:

for i = 1 to 1000 do

*computation(i)*

end do

The goal is to let a quad core cpu split the 1000 independant computations in 4 packages of 250 each to speed up the computation by the factor 4 so that core1 computes i from 1 to 250, core2 from 251 to 500 and so on, at the same time. Is this possible?

A triangle is given.The points

Hi, i have the following question: i've been asked to investigate the following family of functions:

f(x,y)=(be^a(x)^2)+y^2

the question is: How does the shape of each curve depend of a and b.

Thank you

Hello every one,

I want to plot two parameters against each other from a system of ode.

The system is

restart:

with(plots):with(LinearAlgebra):

Eq1:=diff(f(eta),eta$3)+f(eta)*diff(f(eta),eta$2)-diff(f(eta),eta$1)^2-M^2*diff(f(eta),eta$1)+M^2*b+b^2=0;

Eq2:=diff(g(eta),eta$3)+f(eta)*diff(g(eta),eta$2)+(Lambda-2*diff(f(eta),eta$1)-M^2)*diff(g(eta),eta$1)+diff(f(eta),eta$2)*g(eta)=0;

N:=10: M:=0.1:

bc:=f(0)=0,D(f)(0...

Dear friends,

I ran into another problem while using Maple to do residue calculus. The following call

residue(1/x^(1/3)/(x^2+2*x*cos(phi)+1), x=-exp(I*phi));

returns zero -- it does not recognize the residue. On the other hand, if I do a subsitution like this

subs(x=-exp(I*phi), 1/x^(1/3)/(x^2+2*x*cos(phi)+1));

followed by

convert(%, trig); simplify(%);

then I get a divide by zero error, which shows that the...

Hello all,

I recently helped someone with a contour integration problem, you may find my calculation at the following link. It occurred to me to test this integral with Maple and Mathematica.

For example,

int(1/x^(1/3)/(x^2+2*x*cos(Pi/7)+1), x=0..infinity);

produces an unusable output from Mathematica. Maple is a little better,...

I have a procedure which accepts arguments of a (created by myself) type which are records, like so:

Subs:=proc(eqn::seq(equation),elemt::Element,$);
...

Here "Element" is  type I have defined; it is a record with certain fields present.

My question: will this work with option remember? I have the procedure sort-of working without the remember option, but I need to make it work recursively (the Element record can have fields which contain names of other records of type Element...

How to find the maximum of sqrt((x-1)*(y-x))+sqrt((7-y)*(1-x))+sqrt((x-y)*(y-7)) for x=-2..3,
y=0..11 with Maple?
Here are  a few unsuccessful attempts.
> Optimization:-Maximize(sqrt((x-1)*(y-x))+sqrt((7-y)*(1-x))+sqrt((x-y)*(y-7)),
 x = -2 .. 3, y = 0 .. 11);
Error, (in Optimization:-NLPSolve) complex value encountered
> maximize(sqrt((x-1)*(y-x))+sqrt((7-y)*(1-x))+sqrt((x-y)*(y-7)), x = -2 .. 3, y = 0 .. 11);