Int not working with GlobalSolve why? It shows that bounds should be specified to all variable. but i have given bounds to all variables already. Thank you

How do I solve an overdetermined system of algebraic equations in Maple? solve command returns trivial solution for variables which are not actually trivial when I solve them by hand.

Hi all

When I solve using maple the first-order differential equation: diff(y,x)=0  for x in the closed interval [0,1] we obtain a constant function as a solution

but one can define the piecewise constant function see please the attached code

diff_piecewise.mw

why when we differentiate the piecewise function gives undefined derivative at point zero and a half.

Whats is the relationship between this example and Existence and uniqueness theorem for fist order ode

Many thanks

Hello everybody,

While i was trying to work on a physical math problem, a system of 4 integral equations is obtained. The right hand sides of these equations are known functions of r. The left hand sides contain double integrals with respect to lambda and t. i believe that an analytical determination of the 4 unknown functions f_1(t), f_2(t), f_3(t), and f_4(t) is far from being trivial, thus recourse to a numerical technique is necessary and indispensable.

i tried to express the unknown functions as series expansions in t and solve the resulting linear system of equations for the expansion coefficients, but unfortunately the coefficients are very large and the solution is strongly dependent on the number of coefficients. i was wondering whether someone here has some experience with such integral problems and is willing to assist and help. Any hint is highly appreciated.

i attach a Maple script including the equations.

Thank you,

>>>>>> Question.mw

I am trying to use Maple to solve a set of 5 equations, but cannot get a solution. Or there is no solution??

Any help? (Yes, the L function is a likelihood function and I am doing MLE for 5 variables..)

Download PlayGround.mw

Please, is it possible to enlarge the vertices of a graph (DrawGraph in GraphTheory)? Thanks

Hi there, 

I am trying do to a solve procedure, so you can solve an equation using the following argument. 

Solution(eqn, value, var) which uses the embedded solve command. 

When is an integer it works fine with this code here. 

MySolution:=proc(xl::algebraic,is::integer, x::algebraic)
return solve(xl=is,x);
    end proc; 

However the problem arises when I will like to make it to recognize is a an decimal number aka float. 

MySolution:=proc(xl::algebraic,is::integer, x::algebraic)
    return solve(xl=is,x);
    if type(is,float) then
return solve(xl=is,x);
    end proc; 

Any idears on what I am doing wrong? 

Many thanks in advance and have a nice day.

I would like to find a fixed point of f^4 in tems of a and b. I define function as

I calculate f(f(f(f(x,y))))  and Iet f(f(f(f(x,y)))) = (x,y), then use the solve command as:

solve({b^4*(a*(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)^2-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)*b^2*(a*x-x^2-x*y)*x*y)*(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)*(a*x-x^2-x*y)*x*y = y, a*(a*(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)^2-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)*b^2*(a*x-x^2-x*y)*x*y)-(a*(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)^2-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)*b^2*(a*x-x^2-x*y)*x*y)^2-(a*(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)^2-(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)*b^2*(a*x-x^2-x*y)*x*y)*b^3*(a*(a*x-x^2-x*y)-(a*x-x^2-x*y)^2-(a*x-x^2-x*y)*b*x*y)*(a*x-x^2-x*y)*x*y = x}, {x, y})

My computer was freezing. How can I get my result. Thank you

How would I solve a linear algebra differential equation with initial conditions.

For example, what if I had:

x'(t) = [1 2 ; 3 4] x(t) such that the ics: (0,1)

I try to write a Maple Code. But I can' t finish it. 

This is very important for me. Could you help me? download_Maple Code.mw

I am using foslve to get the numerical solution of a very complex equition. When I rerun the code with different Digits settings, I got different solutions. I guess it is not surpprising since it is a numerical solution, but how can I trust/choose which one to use?? Is there anyway to unify the solutions regardless the Digits settings.

This attached are when Digits:=30 and Digits:=20 respectively.. Note, the K is a very large nubmer.. the X is a beta distribution..

Hi guys, 

I have tried to create a loop to solve a set of two equations, but can't seem to get it working. My initial equations are given by;

nstar := (F, L, sigma) -> ceil((ln(k*F) - ln(c(L, sigma)*B))/ln(Phi(L, sigma)))

and 

i := (F, L, sigma) -> r*(1 - (G(L, sigma)*Phi(L, sigma))^nstar(F, L, sigma)*B/F)/(1 - G(L, sigma)^nstar(F, L, sigma))

in which both are based on further rather simple equations. To these I am trying to apply the proc function where I am trying to find which i makes borth the equations above work :

i := proc(F,L,sigma)  

local k :=0.01 ;  

local eps := 0.01 ;  

do while(eps>0.001)  

nstar:= (6)

i := (7)

eps:= i -k:

k=i:

end do;

k;

end proc;

Error, Got internal error in Typesetting:-Parse : "'_Inert_DELAYLESSTHAN' is not a valid inert form"
 

But as you see I am here getting a error which I have not managed to fix. Can anyone see where I might have gone wrong? Could this be done by solve or fsolve? If yes, then how (have tried it as well without succeding)?

help me! 

I have a problem with the system, looking forward to everyone's help!

I have a lot (2000+) of irregular 3D data in an nx3-Matrix "M" (in the added example created via RandomMatrix, in reality actual experimental data), representing 3d coordinates. I manage to create a 3d plot of M via surfdata.

I also have a 3d-function "UC", I want to plot this function when it results into "1". I managed to plot this via implicitplot3d.

So far, so good, see the Maple file. (In reality, the UC function should be approximating the data in M, which obviously isn't the case here because of the randomisation :-)

But I actually want a 2D version of these plots:

• with contour lines for equal values of the third column in M (the vertical axis) for the surfdata.
• with contour lines for equal values of "n" in the UC function.
• a combination of the two above in the same plot (M and UC together, as in the last 3D plot in my example, but in a 2d version). Of course, then there should be some way to make a distinction between the contour lines of M and UC (by using colours or appropriate symbols? Or one of the plots with clour bands, the other with contour lines?)

Is this at all possible? I've been trying by adding "dimension=2", "contour=10", ....: all efforts without succes....

Any help appreciated!

Stef

Download test.mw

expand((x-c)^2+(y-d)^2-R^2) = 0; algsubs(-R^2+c^2+d^2 = f, %); P := proc (x, y) options operator, arrow; -2*x*c-2*y*d+x^2+y^2+f = 0 end proc; 2 2 P := (x, y) -> -2 x c - 2 y d + x + y + f = 0 P(a*cos(theta), b*sin(theta)); G := unapply(%, theta); #usage des formules d'Euler simplify(expand(4*(exp(I*theta))^2*subs(cos(theta) = (exp(I*theta)+exp(-I*theta))*(1/2), sin(theta) = (exp(I*theta)-exp(-I*theta))/(2*I), G(theta)))); poly := sort(subs(exp(I*theta) = X, exp((2*I)*theta) = X^2, exp((3*I)*theta) = X^3, exp((4*I)*theta) = X^4, %)); coeff(lhs(poly), X^4)/tcoeff(lhs(poly)); # exp(I*theta1),exp(I*theta2),exp(I*theta3),exp(I*theta4) sont les racines de ce polynôme unitaire : exp(I*theta1)*exp(I*theta2)*exp(I*theta3)*exp(I*theta4) =1 exp(I*(θ1+θ2+θ3+θ4)=1 d'où θ1+θ2+θ3+θ4 ≡ 2*Pi