Items tagged with solve solve Tagged Items Feed

Hi, I have 10 equation system and 10 unknown variables. I however, want to reduce the equations to 2 with two unknowns. I'm wondering how this could be done in maple. The variables are Y, q, yd, y*, yx, H,  pd, w, P and Pv. I intend to solve the equations for Y and w. 

Thanks in advance for your help. The maple file has also been attached.

 

Regards

Small_Open_Economy.mw

 

hello everyone. I have an undergradute project i'm currently working on and I'm stuck where I have to use the Differential Transforms Method to solve a problem with boundary conditions at infinity


restart;

Digits := 5;

F[0] := 0; F[1] := 0; F[2] := (1/2)*A; T[0] := 1; T[1] := B; M := 2; S := 1;

for k from 0 to 10 do F[k+3] := (2*(sum((r+1)*F[r+1]*(k+1-r)*F[k+1-r], r = 0 .. k))-T[k]-3*(sum((k+1-r)*(k+2-r)*F[r]*F[k+2-r], r = 0 .. k))-M*(k+1)*F[k+1])*factorial(k)/factorial(k+3);

T[k+2] := (-3*(sum((k+1-r)*F[r]*T[k+1-r], r = 0 .. k))-S*T[k])*factorial(k)/factorial(k+2)

end do; f := 0; t := 0;

for k from 0 to 10 do

f := f+F[k]*x^k;

t := t+T[k]*x^k end do;

print(f);
print(t);

but the problem is that i cant seem to evaluate

or higer diagonal pade-approximant. any help will be greatly appreciated. thank you.

How can I solve this equation 18 *9^(x^2 + 2* x) + 768* 4^((x + 3)* (x - 1)) - 5 *6 ^((x + 1)^2)?

I tried

restart:

A:=18 *9^(x^2 + 2* x) + 768* 4^((x + 3)* (x - 1)) - 5 *6 ^((x + 1)^2);
solve(A=0);

I see that, the equation has three solutions: x = -2, x = -1 and x = 0. I check

f:=x->18 *9^(x^2 + 2* x) + 768* 4^((x + 3)* (x - 1)) - 5 *6 ^((x + 1)^2);

f(-2);

f(-1);

f(0);

Another question, Maple can not solve inequality 

18 *9^(x^2 + 2* x) + 768* 4^((x + 3)* (x - 1)) - 5 *6 ^((x + 1)^2) > = 0.

PS. We can easy to solve the above inequality with Mathematica

Reduce[18 9^(x^2 + 2 x ) + 768 4^((x + 3) (x - 1)) - 5 6 ^((x + 1)^2) >= 0 , x, Reals]

I got x <= -2 || x == -1 || x >= 0

 

 

hi guys ,

 

i have two set equations and i want to solve them. eq.mw

 

 

thanks guys

hi everyone,
I wanna to solve this equation using maple 12 but I receive this warning message

cos(x)*cosh(x)=1

cos(x)*cosh(x)=-1

tan(x)=tanh(x)

thanks in advance

Hello, dear experts.
I have a question...
solve the system of differential equations,where one of the initial conditions need to be chosen so thatcondition is metat the end of integration.
The task is not difficult, but I'm having trouble with the syntax.

1.I can't "pull"the desired function from the solution and find its value at a certain point.
I try to do so:
r_ravn:=s->subs(F,r(s));
evalf(r_ravn(s_end));
evalf(r_ravn(0));
but there is no result

2.In this case,instead of"for"it is better to use a while loop, but again the problem arises 1.
Tell me, please,how to implemen my program.

 

restart:
R:=0.3:
theta_min:=Pi/6:
theta_max:=Pi/2:
betta_max:=evalf(Pi/180*80);
p:=2*10^5:

theta0:=s->Pi/3/s_end*s+Pi/6:
r0:=s->R*sin(theta0(s)):
s_end:=evalf(R*(theta_max-theta_min)):

sol1:=solve({sin(betta_max)=c/r0(0)},{c});
const1:=0.1477211630;

betta0:=s->arcsin(const1/r0(s)):
betta:=s->arcsin(r(s)/r0(s)*sin(betta0(s))):
A:=s->cos(betta(s))/cos(betta0(s)):
T1:=s->rT1(s)/r(s):
T2:=s->T1(s)*tan(betta(s))^2:

step:=0.001:
delta:=0.001:
for i from 1 to 3000 do
r_min:=0.3-step:
rT1_n:=p*Pi*r_min^2/2/Pi/sin(theta_min):

sys := diff(rT1(s), s)-A(s)*T2(s)*cos(theta(s)),diff(theta(s), s)-A(s)/T1(s)*(p-T2(s)*sin(theta(s)/r(s))),diff(r(s),s)-A(s)*cos(theta(s)),diff(z(s),s)-A(s)*sin(theta(s));
fcns := {rT1(s),theta(s),r(s), z(s)};
F := dsolve({sys,rT1(0)=rT1_n, theta(0)=theta_min,r(0) = r_min, z(0) = 0}, fcns, numeric,output=listprocedure):
r_ravn:=s->subs(F,r(s)):
if abs(evalf(r_ravn(s_end))-R)=delta then break:
print(r_min):
end if:
end do:

r_ravn:=s->subs(F,r(s));
evalf(r_ravn(s_end));
evalf(r_ravn(0));
plot([r_ravn(s),r(s)],s=0..s_end);

Hi .please help me for solve this nonlinear equations , that attch below

thanks a lots.....

 

 

dsys3 := {8*(diff(f2(x), x, x, x, x))+9*(diff(f2(x), x, x))+10*f2(x)+11*(diff(f1(x), x, x, x))+12*(diff(f1(x), x))+13*(diff(f3(x), x, x))+14*f3(x)+f3(x)*f3(x)+(diff(f3(x), x))*(diff(f3(x), x))+(diff(f3(x), x, x))*f3(x) = 0, 16*(diff(f3(x), x, x, x, x))+18*(diff(f3(x), x, x))+19*(diff(f3(x), x, x))+22*(diff(f1(x), x))+23*(diff(f1(x), x))+24*(diff(f2(x), x, x))+25*f2(x)+26*f2(x)+27*f3(x)+29*f3(x) = 0, 2*(diff(f1(x), x, x))+3*(diff(f2(x), x, x, x))+4*(diff(f2(x), x))+6*(diff(f3(x), x))+7*f1(x)+(diff(f3(x), x, x))*(diff(f3(x), x))+(diff(f3(x), x))*f3(x)+(diff(f3(x), x))*f3(x) = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0, f3(0) = 0, f3(1) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(1) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(1) = 0}; dsol5 := dsolve(dsys3, 'maxmesh' = 500, abserr = .1, numeric, range = 0 .. 1, output = listprocedure)12.mw

I have a set of equations icsi and a set of variables v. How do i assign each solution to the variable before solving the next equation?

I tried- for i from 1 to 3 do solvedi := solve(icsi, vi);assign(solvedi); end. This didn't work for me.

I can't solve this problem, can someone explain me whats worng and how to fix it, here is what i've done:

with(plots);
with(DEtools);
sist := diff(x1(t), t) = -4*x1(t)+x2(t), diff(x2(t), t) = -x1(t)+7*x2(t);
    d                             d                          
   --- x1(t) = -4 x1(t) + x2(t), --- x2(t) = -x1(t) + 7 x2(t)
    dt                            dt                         


condi := [x1(0) = 0, x2(0) = 3], [x1(0) = 2, x2(0) = 1], [x1(0) = 7, x2(0) = 9];
        [x1(0) = 0, x2(0) = 3], [x1(0) = 2, x2(0) = 1],

          [x1(0) = 7, x2(0) = 9]
for i to 3 do sol[i] := dsolve({sist, condi[i]}, {x1(t), x2(t)}, numeric); gra := odeplot(sol[i], [x1(t), t], t = -1 .. 1); display(gra) end do;

when solving a system of equations, I want to get rid of all the absolute functions.

 

for example, |y-2|=x,I don't want maple solve this equation directly...because maple may has difficulties when dealing with absolute values. Instead, I want to transform this equation by squaring the both sides at the same time which end up this equation: (y-2)^2=x^2.

 

The example I provide is kind of simple...what if there are multiple absolute term in the equations? Is there a general way to get what I need? Or Is it practical to use Maple to achive such thing?

 

thanks in advance.

solve 30a+75b+110c+85d+255e+160f+15g+12h+120i=8000 over the positive integers

solve 30a+75b+110c+85d+255e+160f+15g+12h+120i=800000 over the positive integers

restart;
a=1-(1+x)*exp(-x);
SS:=solve(%,x,allsolutions);
evalf(eval(SS,[_Z2=0,a=0.3]));evalf(eval(SS,[_Z2=-1,a=0.3]));

I know there are two branches for the lambertw function meet my demand.

I wonder whether there is a method to find the _Z2 value with SS::real, for the above commands, they are 0&-1.

thank you for any advice.

Hi, I'm trying to solve a system of equation and I keep getting this error. Could anyone help me figure out what I'm doing wrong?

My problem is:

> alpha := .3; G := 3.5; L := 6; f := 1.1;

for i to 50 do

I0 := x(z)+y[i](z); ICon := x(0) = 1, y[i](0) = 0;

for j to 50 do

i <> j;

d1 := diff(x(z), z) = -G*x(z)*y[i](z)/IC-alpha*x(z);

d2 := diff(y[i](z), z) = G*y[i](z)*y[j](z)/IC-alpha*y[i](z);

dsys := {d1, d2};

F := dsolve({ICon, op(dsys)}, [x(z), y[i](z)], numeric);

end do;

end do;
Error, (in dsolve/numeric/process_input) unknown y[2] present in ODE system is not a specified dependent variable or evaluatable procedure


 

Hello! I'm really craving for a piece of advice from you pros.
For simplicity's sake let's start with this made-up system.

where

and

Here is solver options with pair of events to get idea about what I want.

Plot of solution



For my purposes it's essential to restrict variables m and D(fi1)(t), D(fi2)(t). While it's very convenient to use piecewise() operator in first case with m, I simply can't find working method for second.

I've tried various events for dsolve, but they trigger only once and I failed to provide continuous restriction on the derivatives in question.

That's what I need in the end:


If there's a way to obtain this result, please, don't make secrets of it :) Thank you in advance!

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