Items tagged with system system Tagged Items Feed

Dear all;

I need your help to get all the solution of this nonlinear  system  with four parameters: r1, r2, q1, q2  assumed to be positives. Let the system:

the first equation is: r1*x*(1-x/q1)=x*y/(1+x)

the second equation: r2*y*(exp(-y)-q2)=-x*y/(1+x)

How can get the positive solution ( x, y) of the previous system.

Thank you for your help.

Hi

I need your help .. I solve a system and I get the error Warning, solutions may have been lost.

The maple code is attached.

testerror.mw

Thank you for your help.

 

I have a three paramter ode problem that involves three tanks with given initial concentrations.  Overtime the concentration equalizes but one of the steps is to determine all bifurcation values.  Not sure how to continue with this number of variables.

 This is our given system with initial conditions

sys_ode := diff(x(t),t) = (-r*x(t))/100+0+(r*z(t))/50, 
> diff(y(t),t) = (r*x(t))/100+(-r*y(t))/25+0,
> diff(z(t),t) = 0+(r*y(t))/25+(-r*z(t))/50;
> x0:=0; y0 := 200; z0:=0;

Hello people in mapleprimes,

I want to solve the next system of equation for B/A and C/A.

eq1:=A+B=F+G;
eq2:=k*(A-B)=kappa*(F-G);
eq3:=F*exp(I*kappa*a)+G*exp(-I*kappa*a)=C*exp(I*k*a);
eq4:=kappa*F*exp(I*kappa*a)-kappa*G*exp(-I*kappa*a)=k*C*exp(I*k*a);


But, though it is well-known, solve({eq1,eq2,eq3,eq4},{B/A,C/A})
does not work well, as the values I want to solve it for are
expressions: B/A and C/A not variables.

Then, you might thing the next works well.
eq:=subs({B=A/t,C=A/u},{eq1,eq2,eq3,eq4}):
solve(eq,{t,u});

But, this doesn't work well, with the answer was
only the ratio of t and u expressed as the following:

t = t, u = exp(I*k*a)*(exp(-I*kappa*a)*k^2-exp(I*kappa*a)*k^2-exp(-I*kappa*a)*kappa^2+exp(I*kappa*a)*kappa^2)*t/(4*kappa*k*exp(I*kappa*a)*exp(-I*kappa*a))

Isn't there nice way to solve the above system of equation, except that
sol1:=solve({eq3,eq4},{F,G});assign(sol1);
sol2:=solve({eq1,eq2},{A,B});assign(sol2);

Best wishes
taro

I have a nonlinear system with 4 equations and 4 unknowns. I am using fsolve. I know that there are multiple solutions for each variable but am only getting one. I need the others. what do I do??

This is my code:

R__1 := Matrix([[1, 0] , [0, 1] ]);

R__2 := Matrix([[1/2, sqrt(3)/2] , [-sqrt(3)/2, 1/2] ]);

R__3 := Matrix([[-1/2, sqrt(3)/2] , [-sqrt(3)/2, -1/2] ]);

R__4 := Matrix([[-1, 0] , [0, -1] ]);

R__5 := Matrix([[-1/2, -sqrt(3)/2] , [sqrt(3)/2, -1/2] ]);

 

d__1 := Vector( [ 0, 5.4] );

d__2 := Vector( [ 6.4, 4.539] );

d__3 := Vector( [ 11, 4.078] );

d__4 := Vector( [ 15.5, 2.079] );

d__5 := Vector( [ 19, 1.039] );

 

a := Vector( [ a__x, a__y] );

 

A__1:=R__1.a+d__1;

A__2:=R__2.a+d__2;

A__3:=R__3.a+d__3;

A__4:=R__4.a+d__4;

A__5:=R__5.a+d__5;

 

OO:=Vector([O__x,O__y]);

 

DA1:=A__2.A__2-A__1.A__1-2*(A__2-A__1).OO;

DA2:=A__3.A__3-A__1.A__1-2*(A__3-A__1).OO;

DA3:=A__4.A__4-A__1.A__1-2*(A__4-A__1).OO;

DA4:=A__5.A__5-A__1.A__1-2*(A__5-A__1).OO;

 

fsolve({DA1,DA2,DA3,DA4},{a__x,a__y,O__x,O__y});

Thanks for any tips you may be able to offer

 

Hello,

I need your help to display the matrix generated by a lenear system.

The linear system is okay, well written, but  there is an error in generate matrix, I get the wrong matrix.

Thank you very much for your help.

 

test.mw

Hi,

Thank you very much for your idea in previous discussion: this is the link.

http://www.mapleprimes.com/questions/202744-Calculs-Using-Maple

I asked how can I plot the phase portrait of this system

Sys1 := {diff(r(t),t) = r(t)^2*sin(theta(t)), diff(theta(t),t) = -r(t)^2*(-2*cos(theta(t))^2+1)};

I get this answer:

 

1) For me it's not necessary to give an initial condition to plot the pharse portrait.

2) What is the line in this phase portrait.

3) Here, I try a second method using Matlab code, I get:

 

and this is the code:

[r, theta] = meshgrid(0:pi/4:2*pi, 0:pi/4:2*pi);
rdot =r.^2.*sin(theta) ;
thetadot = r.^2.*(1-2*sin(theta).^2);
quiver(r,theta,rdot,thetadot)
xlabel('r')
ylabel('theta')

the arrow in the phase portrait are not the same.

Can someone give me more clear information about this problem.

Many thinks.

 

I am numerically solving a nonlinear system of nine equations. How long can I expect it to take?

I have run it for 30 minutes and it has not solved yet.

Here is the system of equations:

0=Lambda-mu.*S-beta.*(H+C+C1+C2).*(S./N)-tau.*(T+C).*(S./N);

0=tau.*(T+C).*(S./N)-beta.*(H+C+C1+C2).*(T./N)-(mu+mu_T).*T;

0=beta.*(H+C+C1+C2).*(S./N)-tau.*(T+C).*(H./N)-(mu+mu_A).*H;

0=beta.*(H+C+C1+C2).*(T./N)+tau.*(T+C).*(H./N)-(mu+psi.*mu_A+mu_T+lambda_T).*C;

0=lambda_T.*C-(mu+mu_A+rho_1+eta_1).*C1;

0=rho_1.*C1-(mu+mu_A+rho_2+eta_2).*C2;

0=eta_1.*C1-(mu+rho_1+gamma).*CT1;

0=eta_2.*C2-(mu+rho_2+gamma.*(rho_1)./(rho_1+rho_2)).*CT2+(rho_1).*CT1;

0 = N - S - T - H - C - C1 - C2 - CT1 - CT2;

and I have numeric values for Lambda, beta, tau, mu, mu_T, mu_A, rho_1, rho_2, psi, gamma. The only parameters left are eta_1, eta_2.

Thank you.

Having solution of an inequations system, is there a way/function/algorithm to find a particular numeric solution (as simplex[minimize] can do) ?

ex:

Q := {1 < x - y, x + y < 1};

R := solve(Q);

      { x < 1 - y, y < 0, y + 1 < x }

manually it's easy to find some numeric solutions:


      y = -1, x = 1
      y = -2, x = 0

but I need an automatic way.

Thank you for your help
s.py

 

a+a*c

a+a*b

a+b+c

assume system of polynomials are above, how to test whether inverse exist

and find inverse of them

 

I have a set "EQ" containing N linear equations in N unknowns. The only symbolic variables in each "EQ[i]" are the unknowns. I want to write a procedure that derives the matrices "A" and "b" where A.x=b is the same linear system stored in "EQ". In other words, I want to write the linear system in the matrix form.

Can anyone guide me through writting such a procedure?
As an example do it with the system written in the following file.

Note: This procedure will be used for large linear systems (e.g. 2000 Equations, 2000 Unknowns) so it will be important that the procedure uses the least operations required.

LinearSystem.mw

Thanks in advance.

I have a large system of linear algebraic equations that I want to solve (2005 Unknowns, 2005 Equations). I was wondering that what are the proper commands to use in maple for solving the system as fast as possible. Take a look at the files in the download link if you want to see the system of linear algebraic equations.

http://pc.cd/h79

Please provide me any suggesitons that you may think will be helpful like using other sofwares that are good in doing this work such as MATLAB or something else.


Thanks in Advance




I'm trying to solve this system of ODEs by Laplace transform. 

> de1 := d^2*y(t)/dt^2 = y(t)+3*x(t)

> de2 := d^2*x(t)/dt^2 = 4*y(t)-4*exp(t)

with initial conditions 

> ICs := y(0) = 2, (D(y))(0) = 3, x(0) = 1, (D(x))(0) = 2

 

Using 

> deqns := de1, de2

and

> var := y(t), x(t)

 

I need to solve it for both y(t) and x(t), I have tried this by:

> dsolve({ICs, deqns}, var, method = laplace)

And

> dsolve({ICs, deqns}, y(t), method = laplace)

> dsolve({ICs, deqns}, x(t), method = laplace)

 

However I get this error message:

Error, (in dsolve/process_input) invalid initial condition

 

Any help is appreciated

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