I'd like to bring up all the Maple 2016 applications from the application center.

How do I do that?

Dear all,

I would like to ask you to help me with the following problem that I got error.

Thanks in advance

Hello people in mapleprimes,
I want to collect
aa:=sigma*omega*alpha = mu(H, H)*s*omega*gamma*L+omega^(-sigma+1)*tau*mu(F, F)*s*L-omega^(-sigma+1)*tau*mu(F, F)*s*L*gamma;
with omega^(-sigma+1);

I know that the above expression could be collected with tau, which
produces the same expression as the one collected with omega^(-sigma).
And, in collecting with omega^(-sigma+1), I know that freeze-thaw could be used.

But, I want to know another way to do the same thing, using frontend or 
subsindets. Is this possible?

Thanks in advance.

taro
 

To summarize,

dsolve numeric parametric form gives wrong answers when fdiff is used to calculate the Hessian. Doing dsolve twice gives the correct answer. Hope I am not making syntax or programming errors.

Also, fdiff is not compatible with vector form.

Download Bugindsolveorfdiff.mws

hi

please help me for fsolve algebric equations...

thanks

AGM.mw
 

Download AGM.mw

I consider  100 .100 real matrices A,B=Matrix(100,100,(i,j)->rand()) (with 12 significant digits).  In general, ConditionNumber(A) is <10^5; also I choose Digits:=17.Theoretically, the complexity of the calculations of Determinant(A), CharacteristicPolynomial(A,x), A.B and MatrixInverse(A) are similar (~n^3). Yet, the times of these calculations are respectively: 0"13, 0"67, 0"60 and, what surprises me, 75" (moreover, I don't display any result).

My question: concerning the calculation of the inverse, where does this factor 100 come from ? Would Matlab  be 100 times faster ? I do not see why this would be the case; in particular, the standard methods for the calculation of the inverse are  easily programmable.

Thanks in advance.

it would be great if you could explain me why maple cant numerically solve this set of 3 ecuatios for h1, h2, k being my unknows terms, it works with an exponential function. there is a way to solve it?

thanks..
 

Download heat_ecuation.mw

I have code 
with(plots);

implicitplot3d((r*sin(theta)*cos(phi))^2+(r*sin(theta)*sin(phi))^2+(r*cos(theta))^2=36,r=-3..3, phi=0..2*Pi,theta=0..Pi,coords=spherical);

How can I fix the code to display the picture

How to plot f(x)=x in the interval -2<x<2 with a period 4?

The plot should look like this

Hello,
How do i generate 4d plot of f(x,y,z)=1 is it possible made it in 3d?

So I am trying to find the structure of a material from x-ray diffraction data. To accomplish this I have to solve a system of nonlinear equations. Since there is some unaccuracy in the measurements I wanted to use a numerical tool, so I am trying to solve using fsolve. To check that I am doing this correctly I am trying to solve a set of similar equations where I know the solution. My code is shown below.

eqs := {(2^2/a^2+2^2/c^2-2*(2*2)*cos(x)/(a*c))/sin(x)^2 = 1/2.2393^2, (2^2/a^2+4^2/c^2+4*(2*2)*cos(x)/(a*c))/sin(x)^2 = 1/1.5968^2, (1^2/a^2+sin(x)^2/b^2+2^2/c^2+(2*2)*cos(x)/(a*c))/sin(x)^2 = 1/2.7896^2, 1/sin(x)^2*(2^2*sin(x)^2/b^2) = 1/2.8650^2}

fsolve(eqs)

For this system I know that a = 5.44 b = 5.73 c = 7.89 x=pi/2, but fsolve only returns the input.

Am I doing something wrong or should I use a different function? I am quite new to Maple, so have some patience with me.

Thanks in advance.

Hello! Prompt please as in Maple can determine the distribution function of the resulting histogram distribution? I know about cdf function, but how to act in relation to the histogram do not know.

Histogram:

restart;
with(stats);
with(stats[statplots]);
data2 := [30, 30.5, 31, 31.5, 32, 32.5, 32.6, 33, 33.1, 33.3, 33.6, 34, 35, 36];
histogram(data2, area = count);

In other words, I need smoothing the histogram, get graph and get on it to obtain an analytical expression.

Aslam o Alikum!

Hope everyone going fine with sound health I want to make the graph of following system of ODEs

Eq1 := diff(F(eta), eta, eta, eta, eta)-M*(eta*(diff(F(eta), eta, eta, eta))+3*(diff(F(eta), eta, eta))+(diff(F(eta), eta))*(diff(F(eta), eta, eta))-F(eta)*(diff(F(eta), eta, eta, eta)))-Ha^2*(diff(F(eta), eta, eta));
Eq2 := diff(G(eta), eta, eta)+Pr*M*(F(eta)*(diff(G(eta), eta))-eta*(diff(G(eta), eta)))+Pr*Ec*(diff(F(eta), eta, eta))^2+Nb*(diff(H(eta), eta))+diff(G(eta), eta)+Nt*(diff(G(eta), eta))^2;
Eq3 := diff(H(eta), eta, eta)+M*Sc*(F(eta)*(diff(H(eta), eta))-eta*(diff(H(eta), eta)))+Nt*(diff(F(eta), eta, eta))/Nb;
 

using the BC's

IC1 := F(0) = 0, ((D@@2)(F))(0) = 0, (D(G))(0) = 0, (D(H))(0) = 0;
     F(0) = 0, @@(D, 2)(F)(0) = 0, D(G)(0) = 0, D(H)(0) = 0
IC2 := F(1) = 1, (D(F))(1) = 0, G(1) = 1, H(1) = 1;
 

for Ec := .1; Nt := .1; Nb := .1; Sc := .5; Pr := 10; M := .5; and different values of Ha like Ha:=0, Ha:=2, Ha:=4, Ha:=6 and Ha:=8 in one coordinate. 

I am waiting your positive response.

Hello all, first time posting on here so if I make any mistakes please let me know. So I have created a function called "g" and plotted it on a graph. After looking at the two graphs I have plotted (one with discont=false and one with discont=true), I can see 3 discontinuities are present. However, attempting to solve for these discontinuities by setting the denominator of "g" = 0 only returns 1 of the points. I'm not sure why this is happening and how would I go by solving for the other 2 points? Thank you.

Dear all,

My script has multiple animations, in which the range of the variable is larger then I wish to see. Therefore, I have used the view=[x₁..x₂,y₁..y₂] command, which has different values for every plot. Plotted seperately, this gives me exactly what I want. However, if I combine the plots using display(A,B,C) they loose their individual views. Instead, maple chooses the view of plot A, and plots B an C in the range defined by view A. 

Is there a way to combine plots and still maintain individual views, within the largest view?

Kind regards,
Bastiaan Overdorp