Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

How can one use maple to linearized nonlinear ODE of this type Linearize.mw

with maple.

Best regards.

 

Hi,

 

  Suppose I have the following code

 

***

f[1]:=x1^2+x2;
f[2]:=-x2^3+x1-2;

solve({f[1], f[2]}, {x1, x2}):
a,b:= eval([x1, x2], %[1])[]:
a, b;

***

 

It could work and give me "1 and -1". 

 

I would like to have a general subroutine, solving f[1]...f[n]. I cannot input  "solve({f[1], f[2],....,f[n]}, {x1, x2}):" just like that. Is ther any way to use "solve" or similar module to handle arbitary number of variables?

I am trying to use solve to determine 5 unknowns from 12 equations each with seperate data points); however, maple requires you to have equal number equations and variables. Is there any way around this?

There are two ways of expressing the solution to a cubic equation, one of them uses cos and arccos [1]. How do I / is there a way to ask Maple to get this form?

More generally, can Maple be instructucted to solve equations using trig identities?

[1] http://en.wikibooks.org/wiki/Trigonometry/The_solution_of_cubic_equations

I am fitting polynomial differential equations to data, and I came across some troubling behavior. I have sets of coordinates in x, y, z, and values w at these coordinates. The x,y,z coordinates and the data at these coordinates are concatenations of three subsets of data each, and depending on the order of the concatenation, I get different polynomials from the stats fit command. The difference is not trivial -- it makes a huge differentce in the stability of the differential equations. I have written a demonstration worksheet to show this problem.

John Starrett

Solve IVP with complex coef. with compplex varables numerically..

the sys. is x'=-iDelta1x(t)+y(t)+epsilon

y'=-iDelta2y(t)+x(t)z(t)

z'=-2(x*(t)y(t)+x(t)y*(t)), where * means complex conjugate 

I solve it as:

epsilon:=5:Delta1:=4:Delta2:=4:assume(z(t),real):

var:={x_R(t),y_R(t),z_R(t),x_I(t),y_I(t),z_I(t)}:
dsys :={diff(x(t),t)=-I*Delta1*x(t)+y(t)+epsilon, diff(y(t),t)=-I*Delta2*y(t)+x(t)*z(t), diff(z(t),t)=-2*(conjugate(x(t))*y(t)+conjugate(y(t))*x(t))}:
functions := indets(dsys, anyfunc(identical(t))):
redefinitions := map(f -> f = cat(op(0, f), _R)(t) + I*cat(op(0,f), _I)(t), functions):
newsys := map(evalc @ Re, redefinitions) union map(evalc @ Im, redefinitions):

incs := {x_R(0)=0, x_I(0)=0, y_R(0)=0, y_I(0)=0,z_R(0)=-1/2, z_I(0)=0}:
dsol1 :=dsolve({newsys,incs},var,numeric, output=listprocedure, abserr=1e-9, relerr=1e-8,range=0..1):

but it seems there is not runing propebly

 

Hey

I have a excel spreadsheet with a number of data points which I want to plot in Maple.  I have imported the data and I did manage to produce a point plot, but how can I exclude some of the data points when I plot the point plot?
I would, for example, like to exclude the negative values and only plot the data which is greater than zero.

Thanks

Hi,

 

  I have the following code for using "PolynomialSystem" solve equations of polynomial

 

*********

with(SolveTools):

f:=PolynomialSystem({x+y-3, x^2+y^2-5}, {x, y}):

print(x,y);
print(f);
f[1];f[2];
a:=f[1][1];b:=f[1][2];
print(a,b);
c:=-evalf(a);d:=-evalf(b);

****************

 

The output is

***

x, y
{x = 2, y = 1}, {x = 1, y = 2}
{x = 2, y = 1}
{x = 1, y = 2}
x = 2
y = 1
x = 2, y = 1
-x = -2.
-y = -1.

***

From what I have seen, I cannot subtract the values of x and y as 2 and 1. Is there any way that I can get the values of solutions of variables, namely I can assign a variable "a" as 2, and the other variable "b" as 1?

 

Thank you very much!

 

 

 

 

How can one deduce values of f(eta), diff(f(eta),eta) and diff(f(eta),eta$2) from values of f1,f2,f3,f4,.. in the results of the finite difference here Trial.mw . I need the outcome for f and its 1st and 2nd derivatives as listed above. Please help.

How do I compute the Lyapunov exponents for Lorenz system

Is it possible to organize output of plots in Mathematica-like manner

1
in one list (placed in one cell)?

In Maple I use tables for these purposes:

1

Maybe, does there exist a more elegant solution?

hi...plese help me

:

how i can solve set of non linear differential equation with a unknown parameter with numeric role in maple software??

thank a lot ...

Hello, my problem is that I've changed in tools maple default mode to "Maple input" from "2-d" I've aplied it globbaly and it doesnt'work. Still when type Enter it switch automatically to 2-d Math, and I have to change it by Ctrl-M to 1-d. Any ideas how to solve it? I want 2-d just to output.


eqn1 := (3*y/(y^2+1)^(5/2)+(3*(x+y))/(1+(x+y)^2)^(5/2)+(3*(y+z))/(1+(y+z)^2)^(5/2)+(3*(x+y+z))/(1+(x+y+z)^2)^(5/2))*(-2*x^2/(x^2+1)^2+1/(x^2+1)-2*x*(x+y)/(1+(x+y)^2)^2+1/(1+(x+y)^2)-2*x*(x+z)/(1+(x+z)^2)^2+1/(1+(x+z)^2)-2*x*(x+y+z)/(1+(x+y+z)^2)^2+1/(1+(x+y+z)^2))+(-3*x/(x^2+1)^(5/2)-(3*(x+y))/(1+(x+y)^2)^(5/2)-(3*(x+z))/(1+(x+z)^2)^(5/2)-(3*(x+y+z))/(1+(x+y+z)^2)^(5/2))*(-2*y^2/(y^2+1)^2+1/(y^2+1)-2*y*(x+y)/(1+(x+y)^2)^2+1/(1+(x+y)^2)-2*y*(y+z)/(1+(y+z)^2)^2+1/(1+(y+z)^2)-2*y*(x+y+z)/(1+(x+y+z)^2)^2+1/(1+(x+y+z)^2)):

eqn2 := x/(x^2+1)+x/(1+(x+y)^2)+x/(1+(x+z)^2)+x/(1+(x+y+z)^2)-y/(y^2+1)-y/(1+(x+y)^2)-y/(1+(y+z)^2)-y/(1+(x+y+z)^2):

eqn3 := subs({x = (tan(alpha)-tan(beta)+tan(gamma))*(1/2), y = (tan(alpha)+tan(beta)-tan(gamma))*(1/2), z = (-tan(alpha)+tan(beta)+tan(gamma))*(1/2)}, eqn1):

eqn4 := subs({x = (tan(alpha)-tan(beta)+tan(gamma))*(1/2), y = (tan(alpha)+tan(beta)-tan(gamma))*(1/2), z = (-tan(alpha)+tan(beta)+tan(gamma))*(1/2)}, eqn2):

My question is how to solve eqn3 and eqn4 of tan(alpha)&&tan(beta).

I want to solve the equations eqn3 and eqn4 to solve the  tan(alpha)  and tan(beta),  give me a help .thanks a lot 

Download 0112.mw

I just started my Calculus I class and I have a project to prove some various theorems. Maple rocks for this. 

 

How can I insert the Q.E.D. symbol (looks like a little black square) at the end of my proofs?

First 1278 1279 1280 1281 1282 1283 1284 Last Page 1280 of 2224