Items tagged with system system Tagged Items Feed

Dear Maple users

In order to solve a problem in Biology, I have come across four equations with three unknowns, making it overdetermined. The equation system does not have any specific solutions, but since the equations can be considered containing errors, I am looking after the best possible solution - which will minimize the coefficients in the original equation system. My question is twofold:

1. Does there exist a specific Mathematical Theory to handle non-linear equation systems of this type?

2. Which command in Maple can solve it?

 

 

Regards,

Erik

 

Can Maple simplify these DE's by eliminating the d/dt VL(t) by taking the derrivative of the bottom equation and substituting in the first one? 

I solve a linear system of equations which is rank deficient. Naturally, when Maple solves it symbolically, it chooses some of its variables to use them as a basis to express the solution. 

In a specific problem I'm solving, the basis chosen by Maple is -very- smart, showing a good exploitation of the problem structure. 

I'm curious as to what kind of factorization is used by default, or if there's a lot of by hand "black magic" involved, what are its general characteristics. 

 

Best regards

Claudio

how to calculate the polynomial map for a system of  polynomials

assume system of polynomial is in terms of a,b,c

how to find polynomial map

(r - something in terms of a,b,c)

(u - something in terms of a,b,c)

(v - something in terms of a,b,c)

 

A lot of my life is at the moment spent using solve to solve systems of equations, and then trying to weed through the solutions maple gives to find the ones I am interested in. Specifically i'd like to have a program that can weed through the solutions and eliminate those that include equalities of the  form p[i]=-p[j] or p[i]=0  where i and j are integers (or equalities of that form with the letter q replacing p). Specifically i don't want to exclude equalities of the form p[i]=-p[j]*something+something else-another thing.... as they can be useful (or equalities of that form with the letter q replacing p).

Here is a (simple) example of the kind of equations I am likely to be solving and their output from solve:
A := solve([p[1]*p[2]*p[3] = q[1]*q[2]*q[3], p[1]+p[3] = q[1]+q[3], p[2]^2+p[3]^2 = q[2]^2+q[3]^2])

I have some code which gets rid of solutions where one variable is set to 0 

with(ArrayTools);
GetRidOfDumbSolutions := proc (sols)
local Nsols, Npars, GoodSol, GoodSols, GoodSolsCounter, i, j;
Nsols := numelems(sols); Npars := numelems(sols[1]);
GoodSols := []; GoodSolsCounter := 0;
for i to Nsols do
GoodSol := 1;
for j to Npars do
if IsZero(rhs(sols[i, j]))
then GoodSol := 0
end if
end do;
if GoodSol = 1 then
GoodSols := Concatenate(1, GoodSols, sols[i])
end if
end do;
GoodSols
end proc

but i can't see how (in maple) to detect an expression of the form p[i]=-p[j] especiall if that is being written in 2-d math. (i don't quite understand the different maths environments or how to convert from one to another or to string)

Dear Maple T.A. users

I have just begun using Maple T.A. I have access to a number of questions, some of which involves placing points in a coordinate system. It works for most students, but for a few, including myself, it doesn't work. I am not able to place those points in the coordinate system at hand when leftclicking. What can be the reason for this issue?

Erik

 of equations from a set of solution after solve ?

what are these ways?

 

[a = s/RootOf(_Z^2-s^2+s), b = -RootOf(_Z^2-s^2+s)/s, c = RootOf(_Z^2-s^2+s)]

Dear all
Please guide me how to convert system of expressions into system of equations, so as solve them using "solve command".

The following expressions are just coefficients extracted from certain equation.

16*a[2]^4*delta[1]^2-48*a[2]^3*a[3]*delta[1]*delta[2]-2*a[2]*b[2]*delta[1]^2, 48*a[2]^2*a[3]^2*delta[1]*delta[2]-80*a[3]^4*delta[1]*delta[2]-4*a[3]*b[3]*delta[1]*delta[2], 64*a[2]^3*a[3]*delta[1]*delta[2]-64*a[2]*a[3]^3*delta[1]*delta[2]-2*a[2]*b[3]*delta[1]*delta[2]-2*a[3]*b[2]*delta[1]*delta[2]

16*a[2]^4*delta[1]^2-48*a[2]^3*a[3]*delta[1]*delta[2]-2*a[2]*b[2]*delta[1]^2, 48*a[2]^2*a[3]^2*delta[1]*delta[2]-80*a[3]^4*delta[1]*delta[2]-4*a[3]*b[3]*delta[1]*delta[2], 64*a[2]^3*a[3]*delta[1]*delta[2]-64*a[2]*a[3]^3*delta[1]*delta[2]-2*a[2]*b[3]*delta[1]*delta[2]-2*a[3]*b[2]*delta[1]*delta[2]

(1)

It possible for me to write (1) in the following form

for EQ in 16*a[2]^4*delta[1]^2-48*a[2]^3*a[3]*delta[1]*delta[2]-2*a[2]*b[2]*delta[1]^2, 48*a[2]^2*a[3]^2*delta[1]*delta[2]-80*a[3]^4*delta[1]*delta[2]-4*a[3]*b[3]*delta[1]*delta[2], 64*a[2]^3*a[3]*delta[1]*delta[2]-64*a[2]*a[3]^3*delta[1]*delta[2]-2*a[2]*b[3]*delta[1]*delta[2]-2*a[3]*b[2]*delta[1]*delta[2] do EQ = 0 end do

16*a[2]^4*delta[1]^2-48*a[2]^3*a[3]*delta[1]*delta[2]-2*a[2]*b[2]*delta[1]^2 = 0

 

48*a[2]^2*a[3]^2*delta[1]*delta[2]-80*a[3]^4*delta[1]*delta[2]-4*a[3]*b[3]*delta[1]*delta[2] = 0

 

64*a[2]^3*a[3]*delta[1]*delta[2]-64*a[2]*a[3]^3*delta[1]*delta[2]-2*a[2]*b[3]*delta[1]*delta[2]-2*a[3]*b[2]*delta[1]*delta[2] = 0

(2)

But I want to write (1) in the following form

16*a[2]^4*delta[1]^2-48*a[2]^3*a[3]*delta[1]*delta[2]-2*a[2]*b[2]*delta[1]^2 = 0, 48*a[2]^2*a[3]^2*delta[1]*delta[2]-80*a[3]^4*delta[1]*delta[2]-4*a[3]*b[3]*delta[1]*delta[2] = 0, 64*a[2]^3*a[3]*delta[1]*delta[2]-64*a[2]*a[3]^3*delta[1]*delta[2]-2*a[2]*b[3]*delta[1]*delta[2]-2*a[3]*b[2]*delta[1]*delta[2] = 0

16*a[2]^4*delta[1]^2-48*a[2]^3*a[3]*delta[1]*delta[2]-2*a[2]*b[2]*delta[1]^2 = 0, 48*a[2]^2*a[3]^2*delta[1]*delta[2]-80*a[3]^4*delta[1]*delta[2]-4*a[3]*b[3]*delta[1]*delta[2] = 0, 64*a[2]^3*a[3]*delta[1]*delta[2]-64*a[2]*a[3]^3*delta[1]*delta[2]-2*a[2]*b[3]*delta[1]*delta[2]-2*a[3]*b[2]*delta[1]*delta[2] = 0

(3)

``

 

Download Maple_Query.mw

Hello guys and gals!

I'm not strong enough with maple to get what the result I want.

It seems that it's because I'm asking for two lenths, and not a lenth and an angle, but I have no Idea how to tackle it diferently.

If you know a trick, please share it!

 

Here's an image:

http://imgur.com/xavAUoB

And here's the maple file attached (I think)

 complex_problem_from_the_internet.mw

Thanks,

Happy new year!

How to find the value of alpha and beta?

restart

u := proc (x) options operator, arrow; sin(x)+cos(x)-4*x+(alpha-4)*x end proc; v := proc (x) options operator, arrow; sin(x)-cos(x)+beta end proc; `assuming`([solve({alpha = int(u(t)+v(t), t = 0 .. Pi), beta = int(u(t)-v(t), t = 0 .. Pi)}, {alpha, beta})], [alpha <> 0, beta <> 0])

{alpha = 4+beta*Pi+(1/2)*alpha*Pi^2-4*Pi^2, beta = -beta*Pi+(1/2)*alpha*Pi^2-4*Pi^2}

(1)

``

thank you for helping

Download Qu_3.mw

Hi. I want to solve this system of equations by varying the value of n. I managed to aolve and plot for n=1, 1.1 and 1.2 but it happens to be a problem when I let n=1.3.

>restart;

>Digits := 15;

>with(plots):n:=1.3: mu(eta):=(diff(U(eta),eta)^(2)+diff(V(eta),eta)^(2))^((n-1)/(2)):

>Eqn1 := 2*U(eta)+(1-n)*eta*(diff(U(eta), eta))/(n+1)+diff(W(eta), eta) = 0;

>Eqn2 := U(eta)^2-(V(eta)+1)^2+(W(eta)+(1-n)*eta*U(eta)/(n+1))*(diff(U(eta), eta))-mu(eta)*(diff(U(eta), eta, eta))-(diff(U(eta), eta))*(diff(mu(eta), eta)) = 0;

>Eqn3 := 2*U(eta)*(V(eta)+1)+(W(eta)+(1-n)*eta*U(eta)/(n+1))*(diff(V(eta), eta))-mu(eta)*(diff(V(eta), eta, eta))-(diff(V(eta), eta))*(diff(mu(eta), eta)) = 0;

>bcs1 := U(0) = 0, V(0) = 0, W(0) = 0;

>bcs2 := U(20) = 0, V(20) = -1;

>R1 := dsolve({Eqn1, Eqn2, Eqn3, bcs1, bcs2}, {U(eta), V(eta), W(eta)}, initmesh = 30000, output = listprocedure, numeric);


Error, (in dsolve/numeric/bvp) precision is insufficient for required absolute error, suggest increasing Digits to approximately 23 for this problem

>for l from 0 by 2 to 20 do R1(l) end do;
>plot1 := odeplot(R1, [eta, U(eta)], 0 .. 20, numpoints = 2000, color = red);


Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

 

I have tried increasing the Digits as suggested to 23, 25, 30, 31 up untill 500 yet still same error occur suggesting to increase the Digits. Is there any other way to solve this kind of error? Can someone help me? Thank you in advance.

 

Hi. I want to solve a system of equations. But I got this type of error. 

>restart;

>Digits := 15;
>with(plots):n:=0.7:Pr=1: mu(eta):=(diff(U(eta),eta)^(2)+diff(V(eta),eta)^(2))^((n-1)/(2)):
>Eqn1 := 2*U(eta)+(1-n)*eta*(diff(U(eta), eta))/(n+1)+diff(W(eta), eta) = 0:
>Eqn2 := U(eta)^2-(V(eta)+1)^2+(W(eta)+(1-n)*eta*U(eta)/(n+1))*(diff(U(eta), eta))-mu(eta)*(diff(U(eta), eta, eta))-(diff(U(eta), eta))*(diff(mu(eta), eta)) = 0:
>Eqn3 := 2*U(eta)*(V(eta)+1)+(W(eta)+(1-n)*eta*U(eta)/(n+1))*(diff(V(eta), eta))-mu(eta)*(diff(V(eta), eta, eta))-(diff(V(eta), eta))*(diff(mu(eta), eta)) = 0:
>Eqn4 := (W(eta)+(1-n)*eta*U(eta)/(n+1))*(diff(theta(eta), eta))-(mu(eta)*(diff(theta(eta), eta, eta))+(diff(mu(eta), eta))*(diff(theta(eta), eta)))/Pr = 0:
>bcs1 := U(0) = 0, V(0) = 0, W(0) = 0, theta(0) = 1:
>bcs2 := U(20) = 0, V(20) = -1, theta(20) = 0:
>R1 := dsolve({Eqn1, Eqn2, Eqn3, Eqn4, bcs1, bcs2}, {U(eta), V(eta), W(eta), theta(eta)}, initmesh = 20000, output = listprocedure, numeric);

Error, (in dsolve/numeric/bvp/convertsys) too few boundary conditions: expected 8, got 7

>for l from 0 by 2 to 20 do R1(l) end do;
>plot1 := odeplot(R1, [eta, theta(eta)], 0 .. 20, numpoints = 2000, color = red);

 

What is the problem actually because based on the paper that I refer to, there is only 7 bc. 

Can anyone help me?

Thankyou in advance.

 

Consider a taper steel plate of uniform thickness t := 25mm as shown in the Fig. In addition to its self weight, the plate is subjected to a point load P := 100N at its mid point. Find the global force vector [F] , global stiffness matrix [K] , displacement in each element (1 and 2) , stresses in each element  (1 and 2) and reaction force at the support.Take E := 2*10^5N/mm2; rho := 8.2*10^(-5)kg/m3;

restart

t__1 := 150:

t__3 := 75:

w := 25:

l := 600:

t__2 := (t__1-t__3)/l*((1/2)*l)+t__3 = 225/2

A__1 := t__1*w = 3750``

A__2 := t__2*w = 5625/2``

A__3 := t__3*w = 1875``

Revised areas:

A__1e := (A__1+A__2)*(1/2) = 13125/4``

A__2e := (A__2+A__3)*(1/2) = 9375/4``

  E := 2*10^11:m2; F__1 := R__1:is support reaction N; F__2 := 100:N;``

rho__1 := 82*10^(-6) = 41/500000  N/mm2

rho__2 := 82*10^(-6) = 41/500000 N/mm2

l := 600:``

Number of elements,

n__e := 2:

l__e := 300 = 300````

q__0 := 100:N/m ; l := 1: m; n__e := 4:  elementsl  l__e := l/n__e: m;

We shall consider a two element system as shown in the Fig.
For element 1 Stiffness matrix K is

                                           Vector[row](2, {(1) = 1, (2) = 2})
K__1 := A__1e*E/l__e.(Matrix(2, 2, {(1, 1) = 1, (1, 2) = -1, (2, 1) = -1, (2, 2) = 1})) = Matrix([[2625000000000000, -2625000000000000], [-2625000000000000, 2625000000000000]])  Vector(2, {(1) = 1, (2) = 2})

For element 2 Stiffness matrix K is

                                         Vector[row](2, {(1) = 2, (2) = 3})
K__2 := A__2e*E/l__e.(Matrix(2, 2, {(1, 1) = 1, (1, 2) = -1, (2, 1) = -1, (2, 2) = 1})) = Matrix([[1875000000000000, -1875000000000000], [-1875000000000000, 1875000000000000]])  Vector(2, {(1) = 2, (2) = 3})

Global stiffness matrix obtained by adding all the elemental stiffness matrices and given b

           Vector[row](3, {(1) = 0, (2) = 0, (3) = 0})

K__g := Matrix(3, 3, {(1, 1) = K__1[1, 1], (1, 2) = K__1[1, 2], (1, 3) = 0, (2, 1) = K__1[2, 1], (2, 2) = K__1[1, 2]+K__2[1, 1], (2, 3) = K__2[1, 2], (3, 1) = 0, (3, 2) = K__2[2, 1], (3, 3) = K__2[2, 2]}) = Matrix([[K__1[1, 1], K__1[1, 2], 0], [K__1[2, 1], K__1[1, 2]+K__2[1, 1], K__2[1, 2]], [0, K__2[2, 1], K__2[2, 2]]])  Vector(3, {(1) = 0, (2) = 0, (3) = 0})

For element 1 Load matrix F is

  F__1e := (1/2)*`&rho;__1`*A__1e*l__e*(Vector(2, {(1) = 1, (2) = 1})) = Vector[column]([[861/25600], [861/25600]]) Vector(2, {(1) = 1, (2) = 2})

``

For element 2 Load matrix F isNULL

F__2e := (1/2)*A__2e*l__e*`&rho;__2`*(Vector(2, {(1) = 1, (2) = 1})) = Vector[column]([[123/5120], [123/5120]]) 

``

 

Download wrong_answers.mwwrong_answers.mwwrong_answers.mw

Ramakrishnan V

rukmini_ramki@hotmail.com

Math_File.docxMath_File.docx 

I have been using MapleSoft since 2013 to solve mathematical problem. However, I could not solve the following set of differential equations.  Would you please, solve the problem and return the file how you solve it? Thanks in advance.

Hello,

i need help to translate a system which is given below to a for loop.

Other wise i am writing it with myself. 

instead of doing it like that 

with(LinearAlgebra);

sys := [galerkin_funcs[1], galerkin_funcs[2], galerkin_funcs[3], galerkin_funcs[4], galerkin_funcs[5], galerkin_funcs[6], galerkin_funcs[7], galerkin_funcs[8], galerkin_funcs[9], galerkin_funcs[10]];

var := [w[1], w[2], w[3], w[4], w[5], w[6]];

Kmat, Fmat := GenerateMatrix(sys, var);

i want to do it like that.

N:=10

for i to N+4 do

sys(1,i) := galerkin_funcs[i] 

end do

for i to N do

var(1,i) := w[i] 

end do

After that i will generate matrix with this comman Kmat, Fmat := GenerateMatrix(sys, var);

But this for loop i wrote is not doing the i want to do.

Thanks for your help.

 

 

1 2 3 4 5 6 7 Last Page 1 of 17