I have a second order ODE, eq, that depends on a paramter p. Certain events take place when solving the ODE and I would like to study those events as a function of p. Here is a toy version of the command that yeilds the information I need.

>dsolve({eq, f(1)=1,D(f)(1)=0},f(t), numeric, events=[  [ [f(t)-Pi,t>0],halt   ]  ,

 [  [f(t) =-1, t>0],halt]  ])(15);

When f hits  -1 or Pi the solution stops and gives...

Hi

I'm trying to make a phase portrait for coupled oscillators. No problem for the following lines

Dear All:

  I met an strange problem when using fsolve() function in maple. the following is my system equations, ILr0,Vcr0,M,T1 is the values i need to solved by maple . 

   ILr2+ILr0=0,Vcr0+Vcr2=0,ILr1=ILm1,Iout*RL=Vout   (1)  where

   ILr2--> ILr0*f1(T1)+Vcr0*f2(T1)+M*f3(T1);

   Vcr2--> ILr0*g1(T1)+Vcr0*g2(T1)+M*g3(T1);

   ILr1,ILm1-->  ILr0*h1(T1)+Vcr0*h2(T1)+M*h3(T1);

Hello,

I've got a recurrence equation system like:

fn:={f(x)=sum(f(k)*f(x-k-1),k=0..x-1),f(0)=1}

appearently

rsolve(fn,f(k))

doesn't work.

Is there a way to solve this with maple?

Actually I also would like to solve this step by step, by hand. Could you refer me to something that could help me, solve my problem?

Hi

I am getting the floating point error during the job running. it happens instantly for the first fopen command and after some calculations it happens for the second fopen command also.

I would be appreciated if you could help  me in this regards,

Thanks.

here is the code:

> restart;
> w := 0.5695419e-1; T := evalf(2*Pi/w); Cl := 137;
> A := -(int(sin(w*t), t));
> P0 := -subs(t = t1, A)/Cl;
> mrdot0 := int(P0+A/Cl, t = t1 .. t2);

So I am fairly new to Maple and cannot figure this out.

I had to define and equation and then differentiate it and then evaluate this derivative at 2 and now I need to solve the derivative where the derivative=0 for x.  I can do this real easy on paper but I cannot figure out the right set of commands to do so in Maple.  Thank you

Is there a command the takes (A*B)^a and returns A^a*B^a ?

I tried expand and  simplify and also used assuming a>1 but without luck here...

I know this is very basic but I have monster product expressions that I need to be able to raise every term inside the parenthesis by the power so I can isolate some terms of interest.

In MapleSim, it is possible to constrain a rigid body in a multibody system to an arbitrary constraint? 

Similarly to how a mass at the end of a simple pendulum can only move along a circle defined by the length of the pendulum, I was hoping to constrain one to a shape like an ellipse or sine wave.  I could then give it initial velocities and see how it would move along the shape over time.

I'd like to be able to define this constraint with an equation. ...

Hello everybody !

Do you know any mean to find the number of iterations of the simplexe algorithm used in the Minimize pocedure?

Thank you very much for your help

Martin 

How do you plot the following parametric equation (equation of an ellipse) in the same graph

Hi there,

Unfortunately I don't have access to Maple because I am away from my office for the week, but I need a

system of linear equations solved. Can someone please solve it for me and post the solutions (if any)?

The system is:

a_1 x_1 + b_1 x_2 + c_1 x_3 + d_1 x_4 = e_1 + e_4

a_2 x_1 + b_2 x_2 + c_2 x_3 + d_2 x_4 = e_2

a_3 x_1 + b_3 x_2 + c_3 x_3 + d_3 x_4 = 0

a_4 x_1 + b_4 x_2 + c_4 x_3 + d_4 x_4 = e_1 + e_4

I am supposed to adjust the code my teacher gave us for a project. The only problem is I dont know anything about maple or much about what I am supposed to do. I tried copying and pasting the code but I couldn't.

The assignment says to adjust the code to incorporate variable window size. There is a place in the code where it says insert code here and the teacher said it shouldn't need more than like 8 lines of code.

If anybody thinks they they can help I will be more then happy to forward you the email.

I need the list of x and y separately of S[2] like the:

[.8098169753, 0.8098169753, .8098169753, .8098169753, .8098169753, .8098169753]

[1.339721710, 1.309892054, 1.282898312, 1.258277638, 1.235670705, 1.214793049]

restart; with(LinearAlgebra);
a := 2; b := .29; d := 1.85; for h from .5 by .1 to 1 do
eq1 := x*(-b*x^2-x+1);
eq2 := y*((a*x*x)/(b*y^2)-d-h*y);
S := solve({eq1, eq2}, {x, y});

working on my exam notes for tomorrow, I get this...

and only 1/3 or so of my document opens... anyone that can help me out here..?  

I use Maple 17... 

I want to choose four points A, B, C, D on the sphere  (x-2)^2 + (y-4)^2 +(z-6)^2 -81=0 from the list L so that there are not any the right triangle are formed from the points A, B, C, D. How to get it? I tried

> restart:

with(geom3d):

eqS:=Equation(sphere(S,(x-2)^2 + (y-4)^2 +(z-6)^2 -81=0,[x,y,z],'centername'=T)):

L:=[]: