hi, I want to solve linear equations with the assumption that Maple round every step of the evaluation accurate to 2 decimal places. for example: > solve({1.11 x - 2.13 y + 1.14 z = 1.23, 0.32 x + 1.44 y + 3.02 z = 4.10, 1.27 x - 1.14 y + 2.54 z = 5.67}, [x, y, z]) x=19.59 y=7.38 z=-4.24 and not [[x = 19.09759807, y = 7.184806478, z = -4.091838645]] thx for your help

Hello I have tried solving quite simple equations, but maple seems to give me a complex solution or complicated one. Ex.: ode := diff(h(t), t) = -.168*sqrt(h(t)) ics := h(0) = 225 dsolve({ode, ics}) h(t) = RootOf(250*sqrt(_Z)+21*t-250*sqrt(225)) I know it is a correct solution but it seems a bit wrong anyway, with both sqrt(_Z) and Rootof

I am working on a Maximum Likelihood problem and would like to know whether the LinearAlgebra package can be helpful. Up to now I have been using Maple only to check my work. However, someone pointed out that Maple's capabilities reach much further: it also can determine whether the system of equations is consistent and, if so, Maple can solve the equations. Unfortunately, I am only trained as a psychologist and I am still at the beginning of learning Maple's capabilities. Here are the three equations: f = (Y'b - 1t'b)(b'b)^-1 equation#1

I just wnt to know that in maple programing if i have a set of system of equations in which subscripted variables are used then how we find the solution of that system for example
we have set of equations as
eq_set:={-v2,1-a2=0 ,3v0,3+2v2,1=0 ,
-b2+3v3,0-2v1,2=0 , v1,2=0}
here we want to find the values of subscripted variables v2,1 ,v0,3 v3,0 , v1,2.Kindly help me.

Hello, I'd like to use the BVP solver to find the numeric solution to a large set of differential equations. The thing is, I have about 60 equations differing only by constants in the derivitave. Instead of typing in a set 60 equations with different constants, is there an easy way I can create an array of equations and use it with dsolve? Thanks!

I am using the solve command in a proc to solve some simple equations. Is there an easy way of determining if the solution is an imaginary number? (basically I only want real solutions... sometimes the solutions will be in the form of a=b+1, i.e. in terms of another variable, and sometimes it will just be some number... once getting the solutions using something like ans:=solve(eqn,{vars}), how can I determine what type it is? JM

Running into this problem, im trying to graph concentric circles (think planets) with a set radius and a set amount of time used to complete one orbit. I am using parametric equations with t being a value of time in terms of one complete orbit. However, when I try to graph these circles using plots[multiple] I get error codes. Heres my text~ We have our 5 planets, A - E orbiting around the sun affixed at the origin. We can easily graph these together to give us an idea of the orbit and time required to complete one. > plots[multiple](plot, > [0.71 cos(2.985074627 Pi t), 0.71 sin(2.985074627 Pi t), t = 0 .. 67],

March 08 2006
nmcg 0
i have a program to write and i dont know how to do it. if anybody could hwlp me cheers. Q implement a computer program which fora given positive integer n, a given square matrix A of order n and given vector b which n components, decides whether A is a band matrix with p=2=q and in the positive case solves the system of equations Ax=b by simplification of Gauss Elimination for a band matrices with p=2=q

Hello, I am trying to animate a solution to a set of coupled second order differential equations but I am getting what looks like the concatenation of the velocity with the graph I am plotting (note the bump to the left of the solitary wave). I am also looking to make it a bit more precise so I do not see all the distortion as the animation runs, I believe this my be due to the approximation of the numeric solution. Any help would be greatly appreciated.
Sincerely,
M. Hamilton

Hello, I am trying to analyze the numeric solution to the following equations - now what I want to do is a little different, I want to plot u_(n)(some number) as a function of n. This would show me the progression through the particles/distance rather than with respect to time. Here is the code:
> with(DEtools):
> n:=5: #(n can be as large as 500..)
>
> sys:=[seq(diff(u||i(t),t$2)=exp(u||(i+1)(t)-u||i(t))-exp(u||(i)(t)-u||(i-1)(t)),i=2..n-1)]:
> eqn1:=diff(u||1(t),t$2)=exp(u||(2)(t)-u||1(t))-exp(u||(1)(t)-u||(n)(t)):
> eqnn:=diff(u||n(t),t$2)=exp(u||(1)(t)-u||n(t))-exp(u||(n)(t)-u||(n-1)(t)):

Lately I have been experimenting with

structured Gaussian elimination. This is a technique for reducing large sparse systems of linear equations to much smaller dense ones, which can then be solved using a modular method. Needless to say, I had to generate some large sparse linear systems.
I wanted the equations to be written as polynomials, because that is the natural sparse representation in Maple and it makes programming structured Gaussian elimination easier (you can use has and indets, for example). So I tried my favorite randpoly command. This was me trying to generate one linear equation:

Maplets for Calculus is a collection of maplets designed to help students practice their calculus
problem-solving skills and to assist instructors in providing effective classroom demonstrations (including 2- and 3-D visualization -- even animation). The maplets cover all major topics in single-variable calculus - limits, derivatives, integrals, differential equations, sequences, series, and polar coordinates. Some of the maplets help to build intuition and some provide practice with routine computational techniques.
An individual license for Maplets for Calculus is available through the

MapleConnect program at

. Lab/Classroom bundles and site licenses are also available. The complete list of maplets and sample videos may be seen at

.

Hi guys ;)
I have a question. Is it possible to apply a PID controller with the following parameters(P=6,I=4,D=0.02)to a Non Linear System described with the differentiate
equation like this
mk*diff(x(t),t,t)+kfv*diff(x(t),t)-(ki^2*u(t)^2*kc)/(mk*(x(t)-x10)^2)-mk*g=0
or in the second form, as two differentiate equations like
dx1 := (x1,x2) ->x2;
dx2 := (x1,x2) ->((u*ki)^2*kc)/(mk*(x1-x10)^2)-g*mk/mk-kfv*x2/mk;
without linearisation in operating point, where u - is input variable, x - position,
dx/dt - velocity, d^2x/dt - acceleration ?
Thanks for your help.

Hi everybody, I'm new to this site, I did a search and looked through the topics but did not find anything regarding a problem I'm having with Maple 10.02 :
EDIT BELOW: Problem now solved!
While running a "for" loop containing some very very large equations the evaluation would stop and a window would say "Connection to kernel was lost." I finally determined that the mserver.exe process had quit, causing this message. I have configured my firewall to permit all communications that maple and mserver want to do.
The worksheet is a translation from an older version (Maple V R4 ...), the worksheet and the loop ran fine in Maple VR4, but I get the above mentioned problem when running the worksheet in Maple 10. I've translated the code into the updated language, but I still get the error. The loop in question is:

Following are coupled PDEs governing the system.
c1 ∂_{x1}S11 + c2 ∂_{x1}S22 + c3 ∂_{x2}S12 = 0 --- (1)
c2 ∂_{x2}S11 + c1 ∂_{x2}S22 + c3 ∂_{x1}S12 = 0 --- (2)
∂_{x2 x2}S11 + ∂_{x1 x1}S22 - 2 ∂_{x1 x2}S12 = 0 --- (3)
where c1, c2, and c3 are constants. S11, S22, and S33 are 2-dimensional field.
And boundary conditions are appropriately defined.
In fact, Eq.(1) and (2) are the equilibrium equations and Eq.(3) is the compatibility
equation of 2D static strain-stress problem.
I don't have any experiences on constructing finite difference equations of coupled