Advanced Search

Question: how i can fsolve 8 equations....

Question:how i can fsolve 8 equations....

9009134 110 Maple 2016
fsolve
December 28 2016
0

hi .how i can fsolve 8 equations, in which obtained from

for k from 2 to m do eq1[k_] end do

AND

for k from 2 to m do eq2[k_] end do  and so on........

And gain unknown functions as f3[2], f4[4],....

thanks?

fdm-maple.mw
 

 

 ############################Define some parameters

 

 
restart; Digits := 15; n := 1; m := 3; len := 1; h := len/m; nn := m+1
 ############################Define some equation

eq1[k_] := -3.0*h*(-f2[k]*f1[k-1]+f2[k]*f1[k+1]+f1[k]*(-f2[k+1]+f2[k-1]))*f4[k]^2+((-8.0*f1[k]+4.0*f1[k-1]+4.0*f1[k+1])*f3[k]+(-f1[k+1]+f1[k-1])*(-f3[k+1]+f3[k-1]))*f4[k]-f3[k]*(-f1[k+1]+f1[k-1])*(-f4[k+1]+f4[k-1]):

 

 

 

 

                                     ######################################  APPLY BOUNDARY CONDITIONS


f2[0] := f2[2];

1.0

 (1)


``for k from 2 to m do eq1[k_] end do

-1.00000000000000*(-f1[2]*f2[3]+f1[3]*(-f2[4]+f2[2]))*f4[3]^2+((-8.0*f1[3]+4.0*f1[2])*f3[3]+f1[2]*(-f3[4]+f3[2]))*f4[3]-f3[3]*f1[2]*(-f4[4]+f4[2])

 (2)

for k from 2 to m do eq2[k_] end do

(-4.16666666666667*f2[2]^2+(50.0000000000000*f2[3]+8.33333333333333*f2[4])*f2[2]-100.000000000000*f2[3]^2+50.0000000000000*f2[3]*f2[4]-4.16666666666667*f2[4]^2+1/27)*f4[3]^2+((f2[2]-1.*f2[4]-.5*f2[3])*f3[3]+(-50.*f3[2]+50.*f3[4])*f2[2]+(-50.*f3[2]+50.*f3[4])*f2[4]+(100.*f3[2]-100.*f3[4])*f2[3])*f4[3]+((50.*f4[2]-50.*f4[4])*f2[2]+(50.*f4[2]-50.0*f4[4])*f2[4]+(-100.0*f4[2]+100.*f4[4])*f2[3])*f3[3]

 (3)

for k from 2 to m do eq3[k_] end do

-0.123456790123457e-1*f4[3]^3*f3[3]+((1/27)*(.5*f2[2]-.5*f2[4])*f3[3]-0.277777777777778e-1*f2[3]*f3[2]+0.277777777777778e-1*f2[3]*f3[4])*f4[3]^2+(-.222222222222222*f3[3]^2+((1/9)*f3[2]+(1/9)*f3[4]+4.*f2[3]^2+(-4.*f2[2]-4.*f2[4])*f2[3]+f2[2]^2+2.*f2[4]*f2[2]+f2[4]^2)*f3[3]+0.277777777777778e-1*f3[2]^2-0.555555555555556e-1*f3[2]*f3[4]+0.277777777777778e-1*f3[4]^2)*f4[3]+(1/9)*((-.25*f4[2]+.25*f4[4])*f3[2]+(.25*f4[2]-.25*f4[4])*f3[4])*f3[3]

 (4)

for k from 2 to m do eq4[k_] end do

-0.493827160493827e-1*f4[3]^4-(1/27)*(-f2[4]+f2[2])*f4[3]^3+(.111111111111111*f2[3]*(f4[4]-f4[2])-.888888888888889*f3[3]+16.0*(f2[3]-.500000000000000*f2[2]-.500000000000000*f2[4])^2)*f4[3]^2+.444444444444444*((f3[3]+.250000000000000*f3[2]-.250000000000000*f3[4])*f4[2]+f4[4]*(f3[3]-.250000000000000*f3[2]+.250000000000000*f3[4]))*f4[3]-(1/9)*f3[3]*(f4[4]-f4[2])^2

 (5)

``


 

Download fdm-maple.mw

 
Please Wait...