Items tagged with rootfinding

Hello, 

I have just started using Maple, and it seems very powerful. I am trying to solve trigonometric equations and get all the solutions in a range, but when I use fsolve I only get one solution. 

Is this by design of the function or is there another way to do this? 

 

Tom


 

restart

G := 6.6743*10^(-8); 1; R := 1336599.126346; 1; rho := 2.2450*10^14; 1; c := 2.9799*10^10; 1; a := 1/(8/3*(6.67*10^(-8)*Pi*rho/c^2))^.5

0.6674300000e-7

 

1336599.126346

 

0.2245000000e15

 

0.2979900000e11

 

4715700.713/Pi^.5

(1)

y(x) = rho*c((1-(x/a)^2)^(1/2)-(1-(R/a)^2)^(1/2))^2/(3*(1-(R/a)^2)^(1/2)-(1-(x/a)^2)^(1/2))

y(x) = 0.1993516000e36/(3*(1-0.8033593953e-1*Pi^1.0)^(1/2)-(1-0.4496840993e-13*x^2*Pi^1.0)^(1/2))

(2)

``


 

Download y(x).mw

Hi!

I am an error with the use of the function "Analytic" of the packpage RootFinding. These are the procedures:

 

CreaCos := proc (C, n, m, t) local k, F; F := C[1][1]+(C[1][2]-C[1][1])*t; for k to n-1 do F := F, C[k+1][1]+((1/2)*C[k+1][2]-(1/2)*C[k+1][1])*(1-cos(Pi*m^k*t)) end do; return F end proc;

 

Then, for k=50, 100, 150... the instruction

works correctly. However, for higher values of k (for instance, k=250) returns the below error. Some idea or suggets about occurs this error?

Many thanks for your time! 

Error, (in RootFinding:-Analytic) unable to evaluate `@`(evalf, proc (x) option remember; table( [( 0.524900000000000000000000000000e-1+Float(undefined)*I ) = Float(undefined)+Float(undefined)*I ] ) 31250*Pi*sin(62500*Pi*x)/(7/18-(1/2)*cos(62500*Pi*x)) end proc) at the value 0.524900000000000000000000000000e-1+Float(undefined)*I. The expression to be solved was probably not analytic.

 

 

 

hy 
need help 
i made this code but i can not get the answer ,help me to find out where i did wrong.

thanx in advance




restart;
f:=x->(x^3+3*x^2-1);
n:=30;
tol:=1e-9;
a[0]:=0;
b[0]:=10;
Digits :=15;

 

printf("No root F(x) abs(x[i+1]-x[i])\n");

for i from 1 to n do
t[i-1] :=evalf( (b[i-1]-a[i-1])/(f(b[i-1])-f(a[i-1])));
c[i-1] := evalf((a[i-1]*f(b[i-1])-b[i-1]*f(a[i-1]))/(f(b[i-1])-f(a[i-1])));
x[i] :=evalf( x[i-1]-t[i-1]*f(x[i-1])^2/(f(x[i-1])-f(c[i-1])));

printf("%d %10.15f %10.15f %10.15e \n",i,x[i],f(x[i]),abs(x[i]-x[i-1]));
if f(a[i-1])*f(c[i-1])<0 then
a[i]:=a[i-1];
b[i]:=c[i-1];
else
a[i]:=c[i-1];
b[i]:=b[i-1];
if abs(f(x[i]))<tol then
print("approximate solution"= x[i]);
print("No of iterations"= i);
break;
end if;
end if;
end do:

Running the code costs a lot of time, I need some suggestions to make faster and more accurate. Thanks!

sonkok.mw

how i can find order of convergence of newton method by expanding taylor series?? plz send me code???

 

Hi,

I have to find the root of an equation corresponding to the maximum absolute value. I am using root finding package to get all the roots. But after getting all the roots i am not able to apply abs function. Maple sheet is attached.

restart

with(plots):

with(LinearAlgebra):

with(DEtools):

with(ColorTools):

Digits := 30

30

(1)

x := proc (t) options operator, arrow; x0*exp(lambda*t) end proc:

phi := proc (t) options operator, arrow; phi0*exp(lambda*t) end proc:

eqm1 := collect(simplify(coeff(expand(diff(x(t), `$`(t, 2))+(2*0)*beta*(diff(x(t), t))+0*x(t)+n*psi*(-v*(phi(t)-phi(t-2*Pi/(n*omega0)))+x(t)-x(t-2*Pi/(n*omega0)))), exp(lambda*t))), {phi0, x0})

(-n*psi*v+n*psi*v*exp(-2*lambda*Pi/(n*omega0)))*phi0+(lambda^2+n*psi-n*psi*exp(-2*lambda*Pi/(n*omega0)))*x0

(2)

eqm2 := collect(simplify(coeff(expand(diff(phi(t), `$`(t, 2))+(2*0)*(diff(phi(t), t))+phi(t)+n*(-v*(phi(t)-phi(t-2*Pi/(n*omega0)))+x(t)-x(t-2*Pi/(n*omega0)))), exp(lambda*t))), {phi0, x0})

(-n*v+n*v*exp(-2*lambda*Pi/(n*omega0))+lambda^2+1)*phi0+(n-n*exp(-2*lambda*Pi/(n*omega0)))*x0

(3)

mode := simplify(evalc(Re(evalc(subs(lambda = I*Omega, solve(subs(x0 = m*phi0, eqm1), m)))))^2+evalc(Im(evalc(subs(lambda = I*Omega, solve(subs(x0 = m*phi0, eqm1), m)))))^2)

-2*n^2*psi^2*v^2*(-1+cos(2*Omega*Pi/(n*omega0)))/(Omega^4-2*Omega^2*n*psi+2*Omega^2*n*psi*cos(2*Omega*Pi/(n*omega0))+2*n^2*psi^2-2*n^2*psi^2*cos(2*Omega*Pi/(n*omega0)))

(4)

A, b := GenerateMatrix([eqm1, eqm2], [x0, phi0])

A, b := Matrix(2, 2, {(1, 1) = lambda^2+n*psi-n*psi*exp(-2*lambda*Pi/(n*omega0)), (1, 2) = -n*psi*v+n*psi*v*exp(-2*lambda*Pi/(n*omega0)), (2, 1) = n-n*exp(-2*lambda*Pi/(n*omega0)), (2, 2) = -n*v+n*v*exp(-2*lambda*Pi/(n*omega0))+lambda^2+1}), Vector(2, {(1) = 0, (2) = 0})

(5)

with(RootFinding):

eq := subs(n = 6, psi = 1000, omega0 = 1.15, v = 0.1e-1, Determinant(A))

6000.94*lambda^2-5999.94*exp(-.289855072463768115942028985507*lambda*Pi)*lambda^2+lambda^4+6000-6000*exp(-.289855072463768115942028985507*lambda*Pi)

(6)

zeros := RootFinding:-Analytic(eq, lambda, re = 0 .. 400, im = -200 .. 200)

0.899769545162895563524511282265e-56, 0.813609592584011756247655681635e-1-20.6993361029378520006643410260*I, .242743338419727199544214811606-34.4961764258358768825593120288*I, .440964962950043888796944083291-100.074138054178692973033664525*I, .107710271188082726666762251538-106.954651646879437684160623413*I, 1.12290283496379505456476079030-62.0290638297730162295171014475*I, .879463466045683309032252293625-93.2168861049771086211729407830*I, 2.54860869821265794971735119535-80.1919866273564551209847942490*I, 1.52678990439144770439544731898-86.4450560720567958301493690195*I, 2.62945288424037545703549470125-75.0161229879790946191171617450*I, 1.68779005203728587549371003511-68.8012471850312399391042105550*I, .776570081405504740452992339900-55.1681878011205261920670466495*I, 0.851171007270465178285429398270e-9+1.00000500045406723708450960132*I, 0.851171007270465178285445699470e-9-1.00000500045406723708450960133*I, 0.874874719902730972066854301075e-2-6.89997772561385443312823760560*I, 0.354201863215292148351069041542e-1-13.7998152076043523748759861636*I, .369195444156713173497807954493-41.3921704506707022569621870947*I, .540047057129385026999638567235-48.2843908783769449582520027744*I, .149078330738225743331408017894-27.5982749361891156626731068484*I, .369195444156713173497807954500+41.3921704506707022569621870948*I, .440964962950043888796944083291+100.074138054178692973033664525*I, .107710271188082726666762251538+106.954651646879437684160623413*I, 1.12290283496379505456476079030+62.0290638297730162295171014475*I, .879463466045683309032252293625+93.2168861049771086211729407830*I, 2.54860869821265794971735119535+80.1919866273564551209847942490*I, 1.52678990439144770439544731898+86.4450560720567958301493690195*I, 2.62945288424037545703549470125+75.0161229879790946191171617450*I, 1.68779005203728587549371003511+68.8012471850312399391042105550*I, .776570081405504740452992339900+55.1681878011205261920670466495*I, 0.813609592584011756247655681660e-1+20.6993361029378520006643410260*I, 0.354201863215292148351069041261e-1+13.7998152076043523748759861634*I, 0.874874719902730972066854301075e-2+6.89997772561385443312823760560*I, .540047057129385026999638567235+48.2843908783769449582520027744*I, .242743338419727199544214811602+34.4961764258358768825593120288*I, .149078330738225743331408017894+27.5982749361891156626731068484*I

(7)

"zeros.select(int 1)"

Error, missing operation

"zeros.select(int 1)"

 

``


Download question.mw

I will be really thankful for the help.

Regards

Sunit


Here, I attached my maple code. I need to find root. I am using fsolve. But I am not geting the root. Please any one help me... to find the root.

reatart:NULL``

m1 := 0.3e-1;

0.3e-1

(1)

m2 := .4;

.4

(2)

m3 := 2.5;

2.5

(3)

m4 := .3;

.3

(4)

be := .1;

.1

(5)

rho := .1;

.1

(6)

ga := 25;

25

(7)

a := 3.142;

3.142

(8)

q := .5;

.5

(9)

z[0] := 3;

3

(10)

x[0] := 1.5152;

1.5152

(11)

w[0] := 1.1152;

1.1152

(12)

a1 := be*z[0];

.3

(13)

a2 := be*x[0];

.15152

(14)

a3 := rho*w[0];

.11152

(15)

a4 := rho*z[0];

.3

(16)

a5 := rho*w[0];

.11152

(17)

a6 := rho*z[0];

.3

(18)

b1 := a1*a4*ga+a4*ga*m1;

2.475

(19)

D1 := a1+m1+m2+m3+m4;

3.53

(20)

D2 := a1*m2+a1*m3+a1*m4-a2*ga+a3*ga+m1*m2+m1*m3+m1*m4+m2*m3+m2*m4+m3*m4;

1.92600

(21)

D3 := a1*a3*ga+a1*m2*m3+a1*m2*m4+a1*m3*m4-a2*ga*m1-a2*ga*m4+a3*ga*m1+a3*ga*m4+m1*m2*m3+m1*m3*m4+m2*m3*m4+m1*m2*m3;

1.4499000

(22)

D4 := a1*a3*a4*ga+a1*m2*m3*m4-a2*ga*m1*m4+a3*ga*m1*m4+m1*m2*m3*m4;

.3409200

(23)

G1 := -a1*a6-a6*m1-a6*m2-a6*m3;

-.969

(24)

G2 := -a1*a6*m2-a1*a6*m3+a2*a6*ga-a3*a6*ga+a4*a5*ga-a6*m1*m2-a6*m1*m3-a6*m2*m3;

.549300

(25)

G3 := -a1*a3*a6*ga-a1*a6*m2*m3+a2*a6*ga*m1-a3*a6*ga*m1-a6*m1*m2*m3;

-.3409200

(26)

A1 := w^(4*q)*cos(4*q*a*(1/2))+D1*w^(3*q)*cos(3*q*a*(1/2))+D2*w^(2*q)*cos(2*q*a*(1/2))+D3*w^q*cos((1/2)*q*a)+D4;

-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200

(27)

B1 := w^(4*q)*sin(4*q*a*(1/2))+D1*w^(3*q)*sin(3*q*a*(1/2))+D2*w^(2*q)*sin(2*q*a*(1/2))+D3*w^q*sin((1/2)*q*a);

-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5

(28)

A2 := -w^(3*q)*a6*cos(3*q*a*(1/2))+G1*w^(2*q)*cos(2*q*a*(1/2))+G2*w^q*cos((1/2)*q*a)+G3;

.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200

(29)

B2 := -w^(3*q)*a6*sin(3*q*a*(1/2))+G1*w^(2*q)*sin(2*q*a*(1/2))+G2*w^q*sin((1/2)*q*a);

-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5

(30)

C := .27601200;

.27601200

(31)

Q1 := 4*C^2*(A2^2+B2^2);

.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2

(32)

Q2 := -4*C*A2*(A1^2-A2^2+B1^2-B2^2-C^2);

-1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)

(33)

Q3 := (A1^2-A2^2+B1^2-B2^2-C^2)^2-4*C^2*B2^2;

((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)^2-.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2

(34)

V := simplify(-4*Q1*Q3+Q2^2);

-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2)

(35)

x := (-Q2+sqrt(V))/(2*Q1);

(1/2)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)

(36)

E := -2*A1*C*x-A1^2+A2^2-B1^2+B2^2-C^2;

-.2760120000*(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)-(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2+(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2-(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2+(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1

(37)

y := -E/(2*C*B1);

-1.811515442*(-.2760120000*(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)-(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2+(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2-(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2+(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)/(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)

(38)

``

fsolve(x^2+y^2 = 1, w)

fsolve((1/4)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))^2/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)^2+3.281588197*(-.2760120000*(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)-(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2+(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2-(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2+(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)^2/(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2 = 1, w)

(39)

``

 

Download root.mw

 For solving polynomial systems I used RootFinding[Isolate]. But after discussing the question http://www.mapleprimes.com/questions/211774-Roots-Of--Expz--1
I decided to compare Isolate and evalf(solve ([...], [...])). It seemed to me that solve some convenient. The only if in the equation there are integers as a real, they should be recorded with a decimal point. (For real solutions of this procedure should be used with (RealDomain).)  Examples:

SOLVE_ISOLATE.mw

I wonder why then the need Root Finding [Isolate]?

I want to set the determinant of the coefficient matrix equal to zero and then solving for the roots. But I could not achieve it via Maple. Can you help me please? 

You can reach two examples in the following file.  Yeni_Microsoft_Word_Belgesi_(2).docx

 

Besides. how can i compute the following transcental equations via maple 

 

sinh(beta*L)*sin(beta*L)=0

 

cosh(beta*L)*cos(beta*L)-1=0

 

cosh(beta*L)*cos(beta*L)+1=0

 

regards

mehmet

How do I write a procedure to find a root of f(x)=0 in the vicinity of a given value x0. The procedure should initially use the rearrangement method to produce a linearly convergent sequence of values, and should, when appropriate, switch to Aitkin's Method. The input for the procedure should be the re-arranged function and the velue for x0. The output should be the root and the number of iteration taken. The procedure should check that re-arrangement will converge. This program should do in Maple V Release5.

Thank you for your help.

 

I have the following function

where A,B,Ψ, K1,K2,K3,α,β are all constants.

How to find the value of m for which the above expression is 0 or approximate to 0 for different values fo the constants.

e.g., Fixing all the parameters except A, I want to find the values of m for different values of A. How to do that in maple?

 

hi every one..

how i solve numerically  couple equations which attached below .in solve this equation we must  starting from a very small value of V(voltage) with initial guesses for x1 and x3

near zero and using find root method.it is noted that  the solution at this voltage step are used as initial guesses

for the next voltage step, and the process is repeated..

thanks

2DOF.mw

Hello,

Today I've playied a bit with CellDecomposition from the RootFinding package and for one of the systems with which I've playied I got an error which seems to me to be a bug related.

In particular, 

with(RootFinding[Parametric]):

m := CellDecomposition([x^3-y^2 = 0, x^2+y^2-1 < 0], [x, y])

Error, (in RootFinding:-Parametric:-CellDecomposition) Segmentation Violation occurred in external routine

 

Did I make a mistake somewhere or Maple 2015.1 faild?

1 2 3 4 Page 1 of 4