Hi,

> seq(seq(p[j,c],j=2..3),c=1..3)=seq(seq(1/(1+exp(-(mu+cat(tau,j)+cat(eta,c)+mix[j,c]))),j=2..3),c=1..3);

Basically, I want to use these expressions in subs.

How do I tell Maple to make the substitution of "lhs=rhs"?

(in the form of {p[2,1]=something,p[3,1]=something,p[2,2]=something,p[2,2]=something,p[2,3]=something,p[3,3]=something,})

I can then use subs() to replace them in my desired expression.

Casper

Hi all,

I just learned about the with(Threads) package and I wanted to do a bit testing and see if it's 'faster' or more 'efficient'.

Test.pdf

This was done on a i7 940X (QuardCore) laptop.

I dont really see how it is better in a sense?

Or am I using Seq for the wrong purpose?

Casper

Hi all,

I have been trying to improve my code to speed up a bit.

Here is the toy example:

######################

restart:

myfunc:=proc(phi,p,CH)
 # dummy code
 # dummy code
 # dummy code
 # dummy code
end proc;

a:=myfunc(phi,p1,CH1)
b:=myfunc(phi,p2,CH1)
c:=myfunc(phi,p3,CH1)

ANS1:=w1*a+w2*b+(1-w1-w2)*b

a:=myfunc(phi,p1,CH2)
b:=myfunc(phi,p2,CH2)
c:=myfunc(phi,p3,CH2)

I have bracket matching turned on but donot get the rectangular yellow box around matching brackets as shown in the attached picture. Can anyone help me how to turn this feature on for both input and output?

Can somebody explain me this:

I created a null tetrad, but Maple's GRQuery returns false.  Here is the code:

Hello! I have the task to calculate the two angles of an 3d linkage. now in step 11 i have these two equations:

I have to solve these for phi, that is the angle I'm looking for. But the solution for phi should be indepented of psi. alpha and beta are variables. Everything else are known constant values.
What I tried: I solved Gl7 for psi (step 12...

Hi all,

#####################################
myp:=proc(a,b,c)

local ans;

ans:=zip(`+`,a,b);

ans:=ans*c;

return ans;

end proc;
myp([$1..4],[x,y,z,w],2);
map2(myp,[$1..4],{[x,y,z,w],[x2,y2,z2,w2],[x3,y3,z3,w3]},4);

map2(myp,[$1..4],{[x,y,z,w],[x2,y2,z2,w2],[x3,y3,z3,w3]},10);
map2(myp,[$1..4],{[x,y,z,w],[x2,y2,z2,w2],[x3,y3,z3,w3]},0.5);

#####################################

Is it possible to use something like

Hi all,

##################################################################

restart:

K:=50;C:=20;

st:=time():

    for j from 2 to K do
        for c to C do
            solve({log(cat(p,j,'C',c)/(1-cat(p,j,'C',c)))=mu+cat(tau,j)+cat(eta,c)+cat(mix,j,c)},cat(p,j,'C',c));
            assign(%);
  ...

Hi,

Here an example of integral that Maple 17 cannot calculate (mean of standard Gumbel distribution)

> G:=int(x*exp(-x)*exp(-exp(-x)), x = -infinity .. infinity);

The result is known as the Euler-Mascheroni constant as shown by a numerical computation:

 > evalf(G);

0.5772156649

>  identify(%, extension = [gamma]);

gamma

> evalf(gamma)

0.5772156649

Also, it cannot compute the variance

Good morning everbody,

I am currently trying to get string and numerical datas from a txt.file but I don't get how to do this.
I have done this in matlab but I would like to make it in maple.
I have read other post about the subject and try tried this command:

 fd := fopen(dataFile,READ,TEXT):
    theData:=readdata(fd,string, 3): # 3 columns
  fclose(fd):

  A:=convert(theData,Array);

but it didn't help me. 

Is there any way to detect programmatically that a WARNING has been issued? A method that converts WARNINGs to errors would be acceptable.

Hi all,

restart:
a:=<1,2;3,4>;
whattype(a);
b:=Matrix(2,2,[1,2,3,4]);
whattype(b);
evalb(a=b);

Is there any real difference between these two methods of creathing matrices?

How do I check if two matrices are equal?

Thanks,

Casper

I'd like to be able to plot a function that not only has a range for the dependent variable(s) but also for any "variable constants" that may be in the expression.

For example, if I wanted to plot the following expression :

                    y = A sin(wx + c)

and also to explore its behaviour over the following ranges for A, w and c :