Maple Questions and Posts Maple Questions and Posts Feed

These are Posts and Questions associated with the product, Maple

How to solve the system
{sqrt((x-1)^2+(y-5)^2)+(1/2)*abs(x+y) = 3*sqrt(2), sqrt(abs(x+2)) = 2-y}
over the reals symbolically? Of course, with Maple. Mathematica does the job.

I want to solve maximize of equation,but the maximize failed to solve it,who can help me.thanks.

c[1] := (1/8)*w*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+(1/(z^2+1))^(3/2)+1}+(1/8)*{x/((x+y+z)^2+1)+x/((x+y)^2+1)+x/((x+z)^2+1)+x/(x^2+1)}:

c[2] := (1/8)*w*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+[1/(z^2+1)]^(3/2)+1}+(1/8)*{y/((x+y+z)^2+1)+y/((x+y)^2+1)+y/((y+z)^2+1)+y/(y^2+1)}:

t[1] := diff(c[1], x);

(1/8)*w*{-(3/2)*(1/((x+y+z)^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/((x+y)^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/((x+z)^2+1))^(1/2)*(2*x+2*z)/((x+z)^2+1)^2-3*(1/(x^2+1))^(1/2)*x/(x^2+1)^2}+(1/8)*{1/((x+y+z)^2+1)-x*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-x*(2*x+2*y)/((x+y)^2+1)^2+1/((x+z)^2+1)-x*(2*x+2*z)/((x+z)^2+1)^2+1/(x^2+1)-2*x^2/(x^2+1)^2}

(1)

t[2] := diff(c[2], y);

(1/8)*w*{-(3/2)*(1/((x+y+z)^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/((x+y)^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/((y+z)^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2}+(1/8)*{1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2}

(2)

eliminate({t[1], t[2]}, w);

[{w = -{1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2}/{-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2}}, {{1/((x+y+z)^2+1)-x*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-x*(2*x+2*y)/((x+y)^2+1)^2+1/((x+z)^2+1)-x*(2*x+2*z)/((x+z)^2+1)^2+1/(x^2+1)-2*x^2/(x^2+1)^2}*{-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2}-{1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2}*{-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(x^2+2*x*z+z^2+1))^(1/2)*(2*x+2*z)/((x+z)^2+1)^2-3*(1/(x^2+1))^(1/2)*x/(x^2+1)^2}}]

(3)

w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*sqrt(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*sqrt(1/(x^2+2*x*y+y^2+1))*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*sqrt(1/(y^2+2*y*z+z^2+1))*(2*y+2*z)/((y+z)^2+1)^2-3*sqrt(1/(y^2+1))*y/(y^2+1)^2);

w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2)

(4)

sub(w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2), c[1]);

sub(w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2), (1/8)*w*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+(1/(z^2+1))^(3/2)+1}+(1/8)*{x/((x+y+z)^2+1)+x/((x+y)^2+1)+x/((x+z)^2+1)+x/(x^2+1)})

(5)

subs(w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2), c[2]);

-(1/8)*(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+[1/(z^2+1)]^(3/2)+1}/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2)+(1/8)*{y/((x+y+z)^2+1)+y/((x+y)^2+1)+y/((y+z)^2+1)+y/(y^2+1)}

(6)

"#"Iwant to maximize the equation (5)and (6),under the conditon of x,y,z are negative or positive at the same time.

 

NULL

 

Download maximize.mw

Does anyone know how to incorporate the tetrad with the directional derivative? I tried using the SumOverIndices, but get crazy results. I know Maple can find the answer easily because I have done the same thing by hand. What am I missing?

The directional derivative should take the form f,1 = eaμ df/dxμ . The answer is Y,1 = dY/dζ – Ybar dY/du.  I obviously do not get this result.

 


restart; with(Physics); with(Tetrads)

0, "%1 is not a command in the %2 package", Tetrads, Physics

(1)

`#msup(mi("ds",mathcolor = "#af00af"),mn("2",mathcolor = "#af00af"))` := Physics:-`*`(Physics:-`*`(2, dzeta), dzetabar)+Physics:-`*`(Physics:-`*`(2, du), dv)+Physics:-`*`(Physics:-`*`(2, H(zetabar, zeta, v, u)), (du+Physics:-`*`(Ybar(zetabar, zeta, v, u), dzeta)+Physics:-`*`(Y(zetabar, zeta, v, u), dzetabar)-Physics:-`*`(Physics:-`*`(Y(zetabar, zeta, v, u), Ybar(zetabar, zeta, v, u)), dv))^2)

2*dzeta*dzetabar+2*du*dv+2*H(zetabar, zeta, v, u)*(du+Ybar(zetabar, zeta, v, u)*dzeta+Y(zetabar, zeta, v, u)*dzetabar-Y(zetabar, zeta, v, u)*Ybar(zetabar, zeta, v, u)*dv)^2

(2)

X = [zetabar, zeta, v, u]

X = [zetabar, zeta, v, u]

(3)

PDEtools:-declare(`#msup(mi("ds",mathcolor = "#af00af"),mn("2",mathcolor = "#af00af"))`)

Ybar(zetabar, zeta, v, u)*`will now be displayed as`*Ybar

(4)

Setup(automaticsimplification = true, coordinatesystems = (X = [zetabar, zeta, v, u]), metric = 2*dzeta*dzetabar+2*du*dv+2*H(zetabar, zeta, v, u)*(du+Ybar(zetabar, zeta, v, u)*dzeta+Y(zetabar, zeta, v, u)*dzetabar-Y(zetabar, zeta, v, u)*Ybar(zetabar, zeta, v, u)*dv)^2)

[automaticsimplification = true, coordinatesystems = {X}, metric = {(1, 1) = 2*H(X)*Y(X)^2, (1, 2) = 1+2*H(X)*Y(X)*Ybar(X), (1, 3) = -2*H(X)*Y(X)^2*Ybar(X), (1, 4) = 2*H(X)*Y(X), (2, 2) = 2*H(X)*Ybar(X)^2, (2, 3) = -2*H(X)*Ybar(X)^2*Y(X), (2, 4) = 2*H(X)*Ybar(X), (3, 3) = 2*H(X)*Y(X)^2*Ybar(X)^2, (3, 4) = 1-2*H(X)*Y(X)*Ybar(X), (4, 4) = 2*H(X)}]

(5)

g_[]

g_[mu, nu] = (Matrix(4, 4, {(1, 1) = 2*H(X)*Y(X)^2, (1, 2) = 1+2*H(X)*Y(X)*Ybar(X), (1, 3) = -2*H(X)*Y(X)^2*Ybar(X), (1, 4) = 2*H(X)*Y(X), (2, 1) = 1+2*H(X)*Y(X)*Ybar(X), (2, 2) = 2*H(X)*Ybar(X)^2, (2, 3) = -2*H(X)*Ybar(X)^2*Y(X), (2, 4) = 2*H(X)*Ybar(X), (3, 1) = -2*H(X)*Y(X)^2*Ybar(X), (3, 2) = -2*H(X)*Ybar(X)^2*Y(X), (3, 3) = 2*H(X)*Y(X)^2*Ybar(X)^2, (3, 4) = 1-2*H(X)*Y(X)*Ybar(X), (4, 1) = 2*H(X)*Y(X), (4, 2) = 2*H(X)*Ybar(X), (4, 3) = 1-2*H(X)*Y(X)*Ybar(X), (4, 4) = 2*H(X)}))

(6)

``

NULL

NULL

eqn3 := SumOverRepeatedIndices(Physics:-`*`(d_[mu](Y(X)), e_[1, `~mu`]))

((Y(X)*Ybar(X)-1)*(diff(Y(X), zetabar))+(Y(X)*Ybar(X)-1)*(diff(Y(X), zeta))+(diff(Y(X), u)+diff(Y(X), v))*(Y(X)+Ybar(X)))*2^(1/2)/((-(Ybar(X)^2+1)*(Y(X)^2+1)/(Y(X)+Ybar(X))^2)^(1/2)*(2*Y(X)+2*Ybar(X)))

(7)

NULL

``

NULL


Download Directional_Derivative.mw

Dear All,

I'm trying to solve the following in Maple.

minimize(int(0.1e-3+.5*t+0.2e-2*t^2-b*t-a, t = 0 .. 300), location = true)

But Maple told me that the answer is

Float(-infinity), {[{a = Float(infinity), b = Float(infinity)}, Float(-infinity)]}.

I really need to get a kind of numerical answer. Would it be possible? Please Help me!!

https://social.msdn.microsoft.com/Forums/vstudio/en-US/cc2a85ad-30ec-44ed-8c75-636ff71eade2/how-to-convert-integer-or-decimal-number-into-any-base-number?forum=csharpgeneral

1. for example how to convert decimal or integer number into base 3 number, base 5 number etc.

2.how to do logical operation with custom logic table for example,

 

120 special operator 235 

01111000

11101011

 

special operator according to logical table is

1st op 2nd op output
0 0 1
0 1 0
1 0 1
1 1 0

 

  01111000

  11101011

=00010100 = 20

Hi,

I'm trying to solve the following non-linear ODE numerically:

by ececuting

but maple gives me this error-message:

"Error, (in dsolve/numeric/make_proc) Could not convert to an explicit first order system due to 'RootOf'"

I couldnt find any useful information in the manual. What does this error mean? Is there something wrong with my maple code or is there just no solution for this particulare differential equation?

 

Thanks in advance

Hi

I have this PDE and was wondering how I can get Maple to solve it

utt+2ut-uxx=18sin(3πx/l)

with conditions u(0,t)=u(l,t)=0 and u(x,0)=ut(x,0)=0

Thanks

James

 

 

 

Could anyone assist in rectifying this error ''Error, (in fsolve) {f[1], f[2], f[3], f[4], f[5], f[6], f[7], f[8], f[9], f[10], f[11], theta[11]} are in the equation, and are not solved for''. Here is the worksheet FDM_Revisit_1.mw

I am using the SumOverRepeatedIndices and get a Length of Output Exceeded error. Sometimes if I close the file and restart the program then I get a result and no error.  However, if I recalculate then I get the error.

 


restart; with(Physics); with(Tetrads)

[e_, eta_, gamma_, l_, lambda_, m_, mb_, n_]

(1)

`#msup(mi("ds",mathcolor = "#af00af"),mn("2",mathcolor = "#af00af"))` := Physics:-`*`(Physics:-`*`(2, dzeta), dzetabar)+Physics:-`*`(Physics:-`*`(2, du), dv)+Physics:-`*`(Physics:-`*`(2, H(zeta, zetabar, u, v)), (du+Physics:-`*`(Ybar(zeta, zetabar, u, v), dzeta)+Physics:-`*`(Y(zeta, zetabar, u, v), dzetabar)-Physics:-`*`(Physics:-`*`(Y(zeta, zetabar, u, v), Ybar(zeta, zetabar, u, v)), dv))^2)

2*dzeta*dzetabar+2*du*dv+2*H(zeta, zetabar, u, v)*(du+Ybar(zeta, zetabar, u, v)*dzeta+Y(zeta, zetabar, u, v)*dzetabar-Y(zeta, zetabar, u, v)*Ybar(zeta, zetabar, u, v)*dv)^2

(2)

X = [zeta, zetabar, u, v]

X = [zeta, zetabar, u, v]

(3)

PDEtools:-declare(`#msup(mi("ds",mathcolor = "#af00af"),mn("2",mathcolor = "#af00af"))`)

Ybar(zeta, zetabar, u, v)*`will now be displayed as`*Ybar

(4)

Setup(coordinates = (X = [zeta, zetabar, u, v]), metric = 2*dzeta*dzetabar+2*du*dv+2*H(zeta, zetabar, u, v)*(du+Ybar(zeta, zetabar, u, v)*dzeta+Y(zeta, zetabar, u, v)*dzetabar-Y(zeta, zetabar, u, v)*Ybar(zeta, zetabar, u, v)*dv)^2)

[coordinatesystems = {X}, metric = {(1, 1) = 2*H(X)*Ybar(X)^2, (1, 2) = 1+2*H(X)*Ybar(X)*Y(X), (1, 3) = 2*H(X)*Ybar(X), (1, 4) = -2*H(X)*Ybar(X)^2*Y(X), (2, 2) = 2*H(X)*Y(X)^2, (2, 3) = 2*H(X)*Y(X), (2, 4) = -2*H(X)*Y(X)^2*Ybar(X), (3, 3) = 2*H(X), (3, 4) = 1-2*H(X)*Ybar(X)*Y(X), (4, 4) = 2*H(X)*Y(X)^2*Ybar(X)^2}]

(5)

g_[]

g_[mu, nu] = (Matrix(4, 4, {(1, 1) = 2*H(X)*Ybar(X)^2, (1, 2) = 1+2*H(X)*Ybar(X)*Y(X), (1, 3) = 2*H(X)*Ybar(X), (1, 4) = -2*H(X)*Ybar(X)^2*Y(X), (2, 1) = 1+2*H(X)*Ybar(X)*Y(X), (2, 2) = 2*H(X)*Y(X)^2, (2, 3) = 2*H(X)*Y(X), (2, 4) = -2*H(X)*Y(X)^2*Ybar(X), (3, 1) = 2*H(X)*Ybar(X), (3, 2) = 2*H(X)*Y(X), (3, 3) = 2*H(X), (3, 4) = 1-2*H(X)*Ybar(X)*Y(X), (4, 1) = -2*H(X)*Ybar(X)^2*Y(X), (4, 2) = -2*H(X)*Y(X)^2*Ybar(X), (4, 3) = 1-2*H(X)*Ybar(X)*Y(X), (4, 4) = 2*H(X)*Y(X)^2*Ybar(X)^2}))

(6)

NULL

``

eqn1 := SumOverRepeatedIndices(Physics:-`*`(d_[mu](Y(Zeta, zetabar, u, v)), e_[1, `~mu`])) = 0

`[Length of output exceeds limit of 1000000]`

(7)

eqn2 := SumOverRepeatedIndices(Physics:-`*`(d_[mu](Y(Zeta, zetabar, u, v)), e_[2, `~mu`])) = 0

`[Length of output exceeds limit of 1000000]`

(8)

eqn3 := SumOverRepeatedIndices(Physics:-`*`(d_[mu](Y(Zeta, zetabar, u, v)), e_[4, `~mu`])) = x

(1/2)*(-(diff(Y(Zeta, zetabar, u, v), zetabar))*(Y(X)*Ybar(X)+1)*2^(1/2)+(diff(Y(Zeta, zetabar, u, v), u))*2^(1/2)*(Y(X)-Ybar(X))-(diff(Y(Zeta, zetabar, u, v), v))*2^(1/2)*(Y(X)-Ybar(X)))/((-(Ybar(X)^2+1)*(Y(X)^2+1)/(Y(X)-Ybar(X))^2)^(1/2)*(Y(X)-Ybar(X))) = x

(9)

eqn1 := `[Length of output exceeds limit of 1000000]` = 0

`[Length of output exceeds limit of 1000000]` = 0

(10)

algsubs(`[Length of output exceeds limit of 1000000]` = 0, `[Length of output exceeds limit of 1000000]`)

0

(11)

``

simplify(`[Length of output exceeds limit of 1000000]`)

``


Download Derive_Eq_4.4.mw

In brief the problem can be stated as follows:

 

Given dependent variables Qi i=1,...,N and independent variables xi, yi, and zi i=1,...,N

which are related via the following system of N linear equations with parameters P1, P2 and P3 :

Qi = P1xi+P2yi+P3zi   i=1,...,N

How to find the optimal values of  P1, P2 and P3 which satisfy the above system of linear equations subject to the following constraints:

Pi>=0   i=1,2,3

and  P1>=P2P3

 Without the requirement of P1>=P2P3, the problem can be solved with the Non-negative Least Squares Method of Lawson and Hanson.  But with this additional constraint, I am stuck.  

 

Your suggestions are welcome.

 

 

 

 

I am trying to find the root of an equation that involves a procedure and a definite integral (solved numerically). Of course, I don't need the root to be found symbolically, but numerically would be fine. The problem is, I keep getting the error

"Error, (in fsolve) Can't handle expressions with typed procedures"

whenever I try to solve it. Anyone have any ideas? My worksheet is here: Table-1-duplication-mapleprimes.mw 

 

Hi,

I'm trying to fit some parameters to the data of multiple slightly different experiments.

I've written a function which returns the sums of the error for all experiments.

 

The function which has to be minimized does the following:

1) Set parameters for the model ODE based on the input parameters

2) calculate the difference (numeric solution of ODE <-> experimentel data)

3) repeat step 1+2 for all experiments

4) return the sum of all differences

 

The function works as aspected. But when I try to minimize it by calling:

I got an error:

"Error, (in dsolve/numeric/process_parameters) parameter values must evaluate to numeric, got A = A"

For some reason Maple isnt able to set the new parameters of the ODE.

Anyone got an idea how to fix this?

 

Thanks in advance.

 

 

 

I got a problem with solving a second order ODE. 

The ODE is :

-V(xi)+(1/2)*xi*(diff(V(xi), xi))+(1/4)*(diff(V(xi), xi, xi))=-(1/2)*k2*(diff(H(xi), xi))-k1*n*X/E+1+k2

where k1,k2,n,X,E  all are constant.

the condition is :

V(xi) tends to 2*xi^2 as xi tends to infinity.

I used 'dsolve' to solve the equation firstly, and got a solution with two constant C1 and C2, I want to use the condition to elimilate C2, so I used limit(sol,xi=infinity)=2*xi^2. But when I used the command 'limit', I can't get the answer.

Could any one help me? 

Many thanks!!!

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