Items tagged with fsolve fsolve Tagged Items Feed

Hello all,

I try to solve this system using maple 18 by "fsolve", but I don't get the solution, I don't Know what is the problem or What this mean.

Do you have any idea?

 

Best Regards L.Sn=10_R=23.5.mw

Hey there,

I am using the fsolve command in order to solve numerically a system of equations with N equations and N unknowns. According to my discretization the number of equations changes. If I have a small number of equations it all works out perfectly. But if I increase the number of equations I just get something like that:

Sorry, for the long post, but for a small number of unknowns Ai it works. It seems that maple doesnt try to compute? Has anyone encountered the same problems?

 

Any help is appreciated.

Jens

Table_1_for_example_1.mw

I am try to find root by using fsolve. But I am not get solution.

Please help me to solve this problem?

I have been attached the program above.

Thank You.

Best Regards.

Velmurugan G

 

 

In functions such as fsolve, there is an optional parameter that allows one to specify the interval to perform the function on. Additionally, sometimes, if a solution is left out, one can specify an interval to search on to obtain the missing solution.

How does Maple determine the interval to search on if this is not specified?

 

(Additional Question - you don't have to answer this)

Ultimately, I am asking this question because I have a function for Newton's method; however, it requires an interval to run. I have read that fsolve uses Newton's method, so I am curious how to automatically select such an interval. Does anyone know how to implement such a thing?

I am curious about the numerical method(s) used by Maple to calculate fsolve. I've looked at the documentation (https://www.maplesoft.com/support/help/maple/view.aspx?path=fsolve%2fdetails), but the method(s) used are not stated on there. Does anyone know which method Maple uses in fsolve? Additionally, does Maple use multiple methods and if so, how does it determine which one to use? Thanks!

Hi Maple community

I'm running an algorithm where a non-linear equation system must be solved, in this case is a 26x26 system.

After 16116 succesful previous computations, fsolve stops giving me results.
I checked why and I was first expecting that, for some reason, the 26x26 system had an error and I ended with something like 25x26 or vice versa. But that was not the case.

So I tried the command solve and it not only worked fine but also gave me two results, but I only need one. I guess I could check for the wrong solution and discard it, but I still wondering why fsolve is failing and if there is anything to help fsolve not to fail.

These are the set of equations if somebody wants to check them:

EQ[16117][1] := W[1, 16117]*(-0.3860115660e-1*HRa[1, 16117]-0.1876793978e-1*ga[1, 16117]+0.7836678184e-1) = 2.040147478*10^6*SR[1, 16118], W[1, 16117]*(-0.3915554290e-1*HRa[1, 16117]-0.1903748329e-1*ga[1, 16117]+0.8260795999e-1) = 3.876387504, W[1, 16117]*(-0.1876794098e-1*HRa[1, 16117]-0.9892449327e-2*ga[1, 16117]+0.3810204607e-1) = 2.040147478*10^6*v[1, 16118], HLa[1, 16117] = .9724029753*ga[1, 16117]+HRa[1, 16117], NRa[1, 16117] = 0.7006679273e-1*HRa[1, 16117]-.1803623678*ga[1, 16117]+1.002451672, NLa[1, 16117] = 0.7006679273e-1*HRa[1, 16117]+.2484955248*ga[1, 16117]+1.002451672, SL[2, 16118] = SR[1, 16118], fra[1, 16117] = HRa[1, 16117]-HLa[2, 16117], fra[1, 16117] = .25*NRa[1, 16117]+.25*NLa[2, 16117], ga[1, 16117] = 0.;

EQ[16117][2] := W[2, 16117]*(-0.3860115660e-1*HRa[2, 16117]-0.1876793978e-1*ga[2, 16117]+0.7836678184e-1) = -2.040147478*10^6*SL[2, 16118]+7.152482840, W[2, 16117]*(-0.3915554290e-1*HRa[2, 16117]-0.1903748329e-1*ga[2, 16117]+0.8260795999e-1) = 3.876387504, W[2, 16117]*(-0.1876794098e-1*HRa[2, 16117]-0.9892449327e-2*ga[2, 16117]+0.3810204607e-1) = -1.983845478*10^6*SL[2, 16118]+5.221405977, HLa[2, 16117] = .9724029753*ga[2, 16117]+HRa[2, 16117], NRa[2, 16117] = 0.7006679273e-1*HRa[2, 16117]-.1803623678*ga[2, 16117]+1.002451672, NLa[2, 16117] = 0.7006679273e-1*HRa[2, 16117]+.2484955248*ga[2, 16117]+1.002451672, SL[3, 16118] = 0.3505865589e-5, fra[2, 16117] = HRa[2, 16117]-HLa[3, 16117];

EQ[16117][3] := W[3, 16117]*(-0.3860115660e-1*HRa[3, 16117]-0.1876793978e-1*ga[3, 16117]+0.7836678184e-1) = -2.040147478*10^6*SL[3, 16118]+10.82168541, W[3, 16117]*(-0.3915554290e-1*HRa[3, 16117]-0.1903748329e-1*ga[3, 16117]+0.8260795999e-1) = 3.876387504, W[3, 16117]*(-0.1876794098e-1*HRa[3, 16117]-0.9892449327e-2*ga[3, 16117]+0.3810204607e-1) = -1.983845478*10^6*SL[3, 16118]+8.751240594, HLa[3, 16117] = .9724029753*ga[3, 16117]+HRa[3, 16117], NRa[3, 16117] = 0.7006679273e-1*HRa[3, 16117]-.1803623678*ga[3, 16117]+1.002451672, NLa[3, 16117] = 0.7006679273e-1*HRa[3, 16117]+.2484955248*ga[3, 16117]+1.002451672, SL[4, 16118] = 0.5304364281e-5, fra[3, 16117] = HRa[3, 16117];

And after these the solving command that I used was:

SOL[j]:=fsolve({seq(EQ[j][n],n=1..N)},indets({entries(EQ[j],nolist)},assignable(name)));

Which returns

SOL[j]:=

As I said, then I tried the solve command:

SOL[j]:=solve({seq(EQ[j][n],n=1..N)},indets({entries(EQ[j],nolist)},assignable(name)));

which returns:

SOL[16117] :=

{HLa[1, 16117] = 1.011251860, HLa[2, 16117] = .5007913055, HLa[3, 16117] = -0.4240068535e-1, HRa[1, 16117] = 1.011251860, HRa[2, 16117] = .8728245835, HRa[3, 16117] = .2686716410, NLa[1, 16117] = 1.073306847, NLa[2, 16117] = .9685353734, NLa[3, 16117] = .9417827567, NRa[1, 16117] = 1.073306847, NRa[2, 16117] = 1.132612831, NRa[3, 16117] = 1.078974668, SL[2, 16118] = 0.1737463747e-5, SL[3, 16118] = 0.3505865589e-5, SL[4, 16118] = 0.5304364281e-5, SR[1, 16118] = 0.1737463747e-5, W[1, 16117] = 90.12372195, W[2, 16117] = 69.57451714, W[3, 16117] = 49.58407210, fra[1, 16117] = .5104605550, fra[2, 16117] = .9152252689, fra[3, 16117] = .2686716410, ga[1, 16117] = 0., ga[2, 16117] = -.3825916698, ga[3, 16117] = -.3199006320, v[1, 16118] = 8.447574110*10^(-7)},

{HLa[1, 16117] = 3.043461992, HLa[2, 16117] = 2.386862361, HLa[3, 16117] = -0.4240068535e-1, HRa[1, 16117] = 3.043461992, HRa[2, 16117] = 1.087485894, HRa[3, 16117] = .2686716410, NLa[1, 16117] = 1.215697293, NLa[2, 16117] = 1.410701230, NLa[3, 16117] = .9417827567, NRa[1, 16117] = 1.215697293, NRa[2, 16117] = .8376385519, NRa[3, 16117] = 1.078974668, SL[2, 16118] = 0.2032780481e-5, SL[3, 16118] = 0.3505865589e-5, SL[4, 16118] = 0.5304364281e-5, SR[1, 16118] = 0.2032780481e-5, W[1, 16117] = -106.0268094, W[2, 16117] = 265.7250566, W[3, 16117] = 49.58407210, fra[1, 16117] = .6565996307, fra[2, 16117] = 1.129886580, fra[3, 16117] = .2686716410, ga[1, 16117] = 0., ga[2, 16117] = 1.336253076, ga[3, 16117] = -.3199006320, v[1, 16118] = 9.883410782*10^(-7)}

Thanks in advance for any recommendations and suggestions.
 


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

Hi, I have the equations (below) with different parameters. I'd like to find out if is it possitble to:

(1) Create a loop that solves for values of on of the parameters, T, in the interval -0.3<T<0.3 and obtain the solutions allowing for a grid of 0.01.

That is, if one can first solve for T= -0.3 and then use the loop with increments that raise T by 0.01 each time until T=0.3. Other parameters are held constant at their assigned values.

(2) Generates the corresponding vectors for T and all the endogenous variables.

eq1 := PI_T = alpha*M*(kappa/(1+kappa-sigma)-1):
eq2 := l = alpha*((sigma-1)*kappa/(1+kappa-sigma)+1):
eq3 := M = phi_c^(-kappa)*F:
eq4 := e*F = PI_T*v+T:
eq5 := M*l+L_A = L_s:
eq6 := F*e+A = A_s:
eq7 := A = (1-beta)*(L_s+(1-v)*PI_T-T):
eq8 := A_s = L_A:
eq9 := p = sigma/((sigma-1)*phi):
eq10 := P_Y = M^(lambda/(1-sigma))*p:
eq11 := L_s = N*(theta*P_Y^beta)^(-1/delta):
eq12 := phi = (kappa/(1+kappa-sigma))^(1/(sigma-1))*phi_c:
eq13 := U = delta*theta*L_s/((1+delta)*P_Y^beta)+((1-v)*PI_T-T)/P_Y^beta:

Params:= [ T=0, N = 1, v = 1, beta = .75, lambda = 1, sigma = 5, kappa = 4.8, delta = 2, theta = 2.591350635, alpha = 0.3998699153e-1, e = 0.8333333332e-1]:

Init_Values:= {A_s = .2500000000, l = .9996747882, PI_T = 0.8333333332e-1, L_A = .2500000000, phi = 1.878101496, M = .4168022157, A = .1666666667, p = .6655657336, P_Y = .8283396488, L_s = 2/3, phi_c = 1.2, F = 1, U = 1.326439132 }:


SOL:= fsolve(eval({eq||(1..13)}, Params), Init_Values);

 

Many thanks!!

Hi all,

I seem to be quite stuck on figuring out how to leave certain letters (e.g. planck's constant h) inside the equation without having to assign it as some particular number. 

What I am trying to do is find the value of a when the following equation is at a minimum:

E = (a*(h^2)/2m) + 0.3989422804/sqrt(a)

Here h and m are what I want to set as constants without actually setting them to h := 1 because I want a in terms of h and m. I have already found the derivative dE/da:

((h^2)/2m) - 0.1994711402/a^(3/2)

But I cannot use fsolve to find the value of a at the minimum because it keeps saying that h and m are variables and unsolved for.

Any help would be greatly appreciated.

I am trying to find the root of an equation. The problem is, I keep getting the error

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

and 

Warning, solutions may have been lost


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

hi.i am a problem for solving this non linear algebric equation.

please help me...thanks

FSOLVE.mw

FSOLVE.mw

 

Dear all,

I have a question: how to compute the roots of exp(z) = -1 with z in C? 

I tried: 

fsolve( exp(z) = -1, z, complex );

But it only gives one root (0.1671148658e-3+4.934802220*10^9*I) which does not even seem to be correct. I would prefere smth like z_n = I*(2*n-1)*pi or at least multiple roots...

By using

solve(exp(x) = -1, x);

it returns I*Pi.

 

MATLAB MuPAD gives the desired result:


solve(exp(x) = -1, x)

(PI*I + 2*PI*k*I, k in Z)

 

 

Thanks!

How I can solve it ? If I want a solution dependent of a. With fsolve? But how?

-x3+ax2-lnx=0
 

I want to solve one equation with one variable and the variable is also in definite integral delimiter. When trying fsolve, I get the error:

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

code

Here is worksheet.mw

How can I obtain solution with method other from simple manual testing Te values?

hi.please help me for solve algebric equations below?

solve or fsolve dos not any answer?

i think use Newton  Raphson, because solve or fsolve not work.

thankscomparision.mw

restart; mu1 := .1; mu2 := .1; sigma1 := -40; F := 25; upsilon1 := (1/2)*sigma2; upsilon2 := (3/2)*sigma2-sigma1; sigma2 := 100

100

(1)

gamma11 := 2.686901; -1; gamma12 := 7.175339; -1; gamma21 := 2.436735; -1; gamma22 := 12.94855; -1; gamma22 := 12.94855; -1; delta1 := .928207; -1; delta2 := .105073; -1; s11 := .629894; -1; s12 := .217601; -1; s21 := 0.73897e-1; -1; s22 := .805815

.805815

(2)

 

 

 

 
Q1 := -mu1*p1-upsilon1*q1+gamma11*q1*(p1^2+q1^2)+gamma12*q1*(p2^2+q2^2)-delta1*(2*p1*q1*p2-q2*(p1^2-q1^2))-s11*F*q1+s12*F*q2 = 0

-.1*p1-65.747350*q1+2.686901*q1*(p1^2+q1^2)+7.175339*q1*(p2^2+q2^2)-1.856414*p1*q1*p2+.928207*q2*(p1^2-q1^2)+5.440025*q2 = 0

(3)

Q2 := -mu1*q1+upsilon1*p1-gamma11*p1*(p1^2+q1^2)-gamma12*p1*(p2^2+q2^2)-delta1*(2*p1*q1*q2+p2*(p1^2-q1^2))-s11*F*p1-s12*F*p2 = 0

-.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0

(4)

Q3 := -mu2*p2-upsilon2*q2+gamma21*q2*(p1^2+q1^2)+gamma22*q2*(p2^2+q2^2)+delta2*(3*p1^2*q1-q1^3)+s21*F*q1 = 0

-.1*p2-190*q2+2.436735*q2*(p1^2+q1^2)+12.94855*q2*(p2^2+q2^2)+.315219*p1^2*q1-.105073*q1^3+1.847425*q1 = 0

(5)

Q4 := -mu2*q2+upsilon2*p2-gamma21*p2*(p1^2+q1^2)-gamma22*p2*(p2^2+q2^2)+delta2*(-p1^3+3*p1*q1^2)-s21*F*p1 = 0

-.1*q2+190*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1 = 0

(6)

 

-.1*q2+((3/2)*sigma2-40)*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1

(7)

NULL

fsolve({Q1, Q2, Q3, Q4}, {p1, p2, q1, q2})

fsolve({-.1*p1-65.747350*q1+2.686901*q1*(p1^2+q1^2)+7.175339*q1*(p2^2+q2^2)-1.856414*p1*q1*p2+.928207*q2*(p1^2-q1^2)+5.440025*q2 = 0, -.1*p2-190*q2+2.436735*q2*(p1^2+q1^2)+12.94855*q2*(p2^2+q2^2)+.315219*p1^2*q1-.105073*q1^3+1.847425*q1 = 0, -.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0, -.1*q2+190*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1 = 0}, {p1, p2, q1, q2})

(8)

solve(Q1, Q2, Q3, Q4)

Error, invalid input: too many and/or wrong type of arguments passed to solve; first unused argument is -.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0

 

Sol := [fsolve(ZZ, omega)]



Download comparision.mw

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