This is the code  hw2_final.mw

Let me explain it.

I am sure that the mistakes must be about the expresstion of the I1(t) and I2(t). Actually if you delete I1(t) and I2(t) , the whole code works and get the picture at the bottom. 

What I want is to put the expresstion of the I1(t) and I2(t) into 'sol:=dsolve...' and 'plots...' to get the picture of I1(t) and I2(t) with respect to t. Before the t* which subject to Phi(t*)=0 (The blue line in the picture at the bottom is Phi) I want I1(t) and after t* I want I2(t).

I1(t) = (int(sqrt(2*(H(t)+omega*cos(q(t)))), q(t) = q(t)-2*Pi .. q(t), numeric))/Pi.    what I want of this experesstion is to get  'int(sqrt(2*(H(t)+omega*cos(q(t)))' from  'q(t)-2*Pi' to 'q(t)' by numeric method.This q(t) is the solution of the ODE sys.

For example(the number I used is not true,just for example) , at the point t=20, q(t)=30-2*Pi.

so I1(t)= (int(sqrt(2*(H(t)+omega*cos(x))), x = 30-2*Pi .. 30, numeric))/Pi.The I2(t) I want is similar to I1(t).

How can I solve it?

Hello, 

I am pretty new to Maple and Im trying to find the parameters of an equation using some kind of fit routine but I can only find such a routine to fit an expression to data and not vice versa.

my equation is as follows

epsilon(E):= a0 + a1*ln(E) + a2*ln(E)^2 + a3*ln(E)^3

I do have data for E which I imagine I need. 

E:= {121, 244, 344, 411, 444, 778, 867, 964, 1085, 1112, 1212, 1299, 1408} all in keV :)

Any suggestions/help would be much appreciated, although I am new to the program so go easy on me :) 

Cheers

Integral.mw

Hi all

Please help me to solve this Integral...

There have come unwanted lines and marks . I donot know how to remove them. Using doc.block, remove block seems to be little tough to incorporate! Please enlighten me. Modified doc. is most welcome. Thanks. Ramakrishnan V 

Download GausianFinal15Nov2015.mwGausianFinal15Nov2015.mw

There are the complexes (C), quaternions (H or Q), octionions (O), sedenions (S) and the pathions (P). I have found the multiplication tables of them, although according to signs (+ or -) there differents at pathions. The important question is that How can I multiply two bases, i_n and i_m of higher dimensions, like in the routions or in the voudions?

Should I xor the indexes of the bases? Like this way: i_1 * i_2 = i_(1^2) = i_3

What is about the signs?

hello dear
how can i solve this three couplde pde?please help me

I want to solve an ODE from Game Theory, the Cournot competition.

It says

p(q1+r2(q1))+p'(q1+r2(q1))*r2(q1)-c2'(r2(q1))=0

 where, I think,

' means diff(,q1),

c2(q2)=c*q2 for a fixed c in [0,1]

and

p(q)=max(0,1-q).

So c2,p and r2 are functions.r2 goes from [0,inf) to [0,inf).

I look for r2, which should be r2(q1)=(1-q1-c)/2 when correctly solved.

However, the command dsolve says Error in dsolve (divison by 0).

 What is wrong? How do I obtain the solution for r2 in Maple?

How many numbers exist which are: less than 1000, multiple of 18 and also multiple of 10? Use Maple to find the answer 

Hello,

I tried to solve below equation, but it gives me zeros result. Please help me to find their inverse laplace. 

 It will be clearer if was pasted on Maple:

restart

Ps := [P[0], P[1], P[2], P[3], P[4]]:

eqs := [P[0](s) =~ (P[1](s)*mu[1]+P[2](s)*mu[2]+P[3](s)*mu[3]+P[4](s)*mu[4])/(s+lambda[1]+lambda[2]+lambda[3]+lambda[4]), P[1](s) = lambda[1]*P[0](s)/(s+mu[1]), P[2](s) = lambda[2]*P[0](s)/(s+mu[2]), P[3](s) = lambda[3]*P[0](s)/(s+mu[3]), P[4](s) = lambda[4]*P[0](s)/(s+mu[4])];

Ls := solve(eqs, Ps(s))[];

P(t)=~inttrans[invlaplace]~(rhs~(Ls), s, t);

Thank you

I'm trying to run c code generated from a simple MapleSim example (Spin Up Beam). To makes things even simpler, I set the input real expression to just time, and changed the input function in my c code to:

static void inpfn(double T, double *U)
{
U[0] = T;
}

I set all the state variables to 0 for the initial conditions. I ran the code and my output doesn't resemble the MapleSim output at all, and I'm getting errors. So my question is about the state variables, which I initially set to 0:

* State variable(s):
* x[ 0] = `Main.SpinUpBeam.FB3.uf1`(t)
* x[ 1] = diff(`Main.SpinUpBeam.FB3.uf1`(t),t)
* x[ 2] = `Main.SpinUpBeam.FB3.uf2`(t)
* x[ 3] = diff(`Main.SpinUpBeam.FB3.uf2`(t),t)
* x[ 4] = `Main.SpinUpBeam.FB3.vf1`(t)
* x[ 5] = diff(`Main.SpinUpBeam.FB3.vf1`(t),t)
* x[ 6] = `Main.SpinUpBeam.FB3.vf2`(t)
* x[ 7] = diff(`Main.SpinUpBeam.FB3.vf2`(t),t)

What do these represent for this flexible beam (FB3)? I.e. what is uf1, uf2, vf1, vf2? Should they be set to 0, or no?

Thanks.

I have several *.mw files which I have built seperately. Now, for readability reasons, I would like to create a Master document, which will ues the first files, and then generate a PDF document.
I tried the read command, but it does not recognize the files.

Any help is appreciated.

Wonderful.mw

Hi all

I don't know what is its meaning???

The plot and the solved Zero-roots aren't coincident with each other...

Zero-roots: 7.532332868 Wheras from the plot we noticed that there's a root near 14...

And why this Code can't find other Zero-roots?

Thanks a lot

please I want a clear explanation for the command (contourplot) in the maple?

contourplot

Hi Mr Preben Alsholm .

i before asked you about a problem.now i have another question...

for determine unknown parameter oemega(or eigenvalue) according to line ''newsys2 := subs(omega^2 = omega2*10^5,newsys)'' and also line ''approxsoln = [omega2 = 1, f(x) = x^2*(1-x)^2]''

power (5) in term'' 10^5''  in  (omega2*10^5,newsys) and initial guess (1)  (omega2 = 1) are very effective on final result eigenvalues .sometimes gain different eigenvalues ,which it is not impossible recognize that

which one is correct and convege.however i need Positive

 eigenvalue that is minimum between them,If and only if , it converge after
some iteration in two section of maple file which is attached as
Eigenvalue.mw

  According to my code converge occuerd until  eigenvalue gained in first section (before line `11` := ((1+6*alpha2)*(1/12))*(int(fy11^2, x = 0 .. 1))...........) are equal with those obtain in second section .(please see below file for example for_example.mw)

it is necessary mention that between section 1 and 2 is relations.for example amont of first obtained eigenvalue which obtain from first section must be repleaced in line ''approxsoln = [omega2 = 0.661514014001420, h(theta) = theta^2*(1-theta)^2]...........''in second section,

and this procedure should be continued between section 1 and 2 until convergence with desired accuracy occurred.Another relate is that,the first ODE system  can be solved using the first set of boundary data to obtain first estimate for  second section (fy11). At the next step by repeating in the same manner, 

this time by obtained function (g) at the end of maple code (Solution of the second ODE system  with second set of boundary conditions leads to the first estimate of function g),that at this stage, the first iteration is completed.

Now by replacing omega2 which is determined in this section in to the first section ,fy11 is updated and gained Further.

Next, the updated function , by continuing the iterative procedure.

in matlab bvp4c rule is used for this purpose.is this impossible in maple software ?if not please help me for solve and gain correct omega..

thanks alot

When I right click on a plot and add a label name for the vertical axis, it is shown in horizontal way next to it, what should I do if I want to have it in a vertical display along the axis written upward?