Hi, I am trying to find adomian's polynomial of exp(y), but after execution it shows DF(e^y), D^2F(e^y) as it is. why it does'nt show its derivative?plz help meNewADP.2.mw

Hi, I am working on an assignment and I am having difficulties with this Eurler method since I am working with vector.

The assignment:

If the aircraft's engine transmits power P to the air, the force of the propellar, Fp, satisfies P=Fpv. Since the force Fp acts in the direction of motion, it can be written in vector form as: (what I have so far)

with(linalg);

v:=Vector(2,1,[vx,vy]); velocity vector

nu:=v->Norm(Vector(v),2); (length of v)

Fd:=-1/2*C*rho*A*nu*v # air resistance

Fl:=1/2*Cl*rho*A*nu*v # liftforce in y direction here v=[-vy,vx]

Fp:=(P/nu^2)*v^2 # power from aircraft

All this adds up to:

F_T:=v->Vector(evalm(Fp/nu(v)^2*v+Fd*nu(v)*v+Fl*nu(v)*zip((x,y)->x*y,v([2,1]),Vector([-1,1]))+Vector(2,1,[0,-m*g]))); # I have tried it and it works

so up til here I am good. It is this next part that I can't seam to get.

Write a routine which uses Euler's method to solve the equation of motion numerically for velocity v(t):

m*dv/dt=F_T(v)

Test the program by trying a plane with intitial velocity [v0,0] such that it maintains this speed when the engines power is P=P0

We have been given an Esolve we can use:

Esolve:=proc(f::procedure,h,x0,y0,N::integer) # Calculates y(x0+n*h) for n=1..N given y0=y(x0) and f(x,y)

local n,y,x;

y[0]:=y0; # Start Value
x:=x0;

for n from 1 to N do
y[n]:= y[n-1]+h*f(y[n-1],x); # Euler's formula
x:=x+h; # Next x
end do;

y;

end proc;

If anyone could help me understand how it works and how to get this to work with vectors I would be very gratefull.

Vahid Borji

How can I plot the p stability region in H₁and H₂ region I already gotten the polynomial 

Hi every body.

I want to solve a differential equation:

ode=D(r(θ),θ).....

I have the initial conditions:

e=10^(-20)

r(e)=A where A is parameter.
r'(e)=0

sol:=dsolve({ode,r(e)=A,D(r)(e)=0},parameter=[A],numeric)

and:

All in the file.

#@#

Ode-Problem.mw

hi every one

please see attached file blew and help me for dsolve set equations.i dont know aboue errors with in

thanks

equation.mw

When solving ODE with dsolve, I don't want to use "_C1" as the name of constant, I want to specify by myself, how can I do?

Take for instance the signal [1,0,-3,2,1,0,1,2] (considered in Wavelets for Kids, Vidakivic & Mueler, AMS, 1991).

I want to anlalyse it down to the level zero with Haar wavelets. It seems that MAPLE's commnad only applies once.

Hi !

I am trying to solve a pde with initial/boundary conditions, in the numeric mode. It works very well when I provide 3 numerical conditions : 1 initial condition (uniform distribution at first), and two boundary conditions.

Now I want to switch the initial condition to a function of r. It's a polynome I obtained after interpolation of another result. I tested it and the function interp20(r) works. But the pdsolve doesn't seem to evaluate the function, when it comes to start the procedure and pdsolve doesn't return a module as it should, although it doesn't send any message error.

You can see the code following, with the error messages I get.

Is it strictly impossible to use a function as an initial condition ?
Is it just a problem I could solve by converting the function or its result to something else ?
Is float the right type of input ? For example, if I write a:=283.15, is it a float ?
Have you already had similar problems ? How did you solve them ? Where could I find working code examples on this ?
Could I solve this problem with non-uniform initial condition with another Maple function ?

Then you very much for your help !

Hello 

Does any one have the package DESOLV not deslove , I need it with its help file.

Carl Love 9260

hi.i have problem for dsolve equation.please help me again

thanks alotequation.mw

Hi everybody.

In the following attached file, I have created 3 matrices with 13*13 and 13*6 dimensions and their components are numeric. I closed Maple and ran it again. Surprisingly, when I try to browse the matrices again, I encounter with this problem " empty RTable structure" while I imported the components, manually.

Thanks in advance for your answers.

Ho.mw

I want to solve these two differential equations. I have the initial conditions:
x(0)=0
y(0)=0
x'(0)=5.7
t'(0)=8.1
What am I doing wrong?

How can I do a color inversion on a jpg image in Maple?

Hello 

I have maple 13, and I installed two new packages, the problem that I don't know where should I saved them and how to make them active to start using there routines ??? 

So pleas any one know how to add packages to maple 13, help me    :)