i have system of 15 equations and 15 unknown when i use (vpasolve) after 45 minit i have solutions with Warning: Solutions might be lost. [solvelib:​​:allvalue​s​_warning​]

2 views (last 30 days)
sadeem alqarni
sadeem alqarni on 6 Sep 2018
Commented: mbvoyager on 7 Sep 2018
syms y1 y2 y3 y4 yn1 yn11 yr yr11 yr111 ys ys1 ys11 yng yng1 yng11 format longe h=0.01 r=1.8594945668873855539277962650337 %1694/911=1.859495060373216e+00 s=1.5291603115332358052897544638033 %1075/703=1.529160739687056e+00 g=0.41887850354106247407749038324581 y0=0 y01=0; y011=1; fnr= sym(-1/2*(yr111)*(yr )) fns=sym(-1/2*ys11*ys) fng=sym(-1/2*yng11*yng) fnx=sym(-1/2*y011*y0) fnxh=sym(-1/2*yn11*y1) fn2=sym(-1/2*y4*y2)
eq1=-y2-(r - 2)/r*y0+ 2/(r*(r - 1))*yr+ (2*r - 4)/(r - 1)*y1+ (h^3*(r - 2)*(10*g + 19*r - 48*g*r + 24*g*r^3 - 6*g*r^4 + 3*r^2 - 5*r^3 - 9*r^4 + 3*r^5 + 9))/(840*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(10*g + 9*r + 10*s - 6*g*r - 42*g*s - 6*r*s - 14*g*r^2 - 18*g*r^3 + 8*g*r^4 - 14*r^2*s - 18*r^3*s + 8*r^4*s + 9*r^2 + 9*r^3 + 9*r^4 - 5*r^5 + 42*g*r^2*s - 14*g*r^3*s + 14*g*r*s + 9))/(840*r*(g - r)*(r - s)*(r - 1))*fnr -(h^3*(r - 2)*(172*g + 87*r + 172*s - 106*g*r - 182*g*s - 106*r*s - 42*g*r^2 - 10*g*r^3 + 6*g*r^4 - 42*r^2*s - 10*r^3*s + 6*r^4*s + 39*r^2 + 15*r^3 + 3*r^4 - 3*r^5 + 42*g*r^2*s - 14*g*r^3*s + 154*g*r*s - 181))/(840*(g - 1)*(r - 1)*(s - 1))*fnxh -(h^3*(r - 2)*(10*g + 19*r + 10*s - 48*g*r - 42*g*s - 48*r*s + 24*g*r^3 - 6*g*r^4 + 24*r^3*s - 6*r^4*s + 3*r^2 - 5*r^3 - 9*r^4 + 3*r^5 - 84*g*r^2*s + 14*g*r^3*s + 140*g*r*s + 9))/(1680*g*r*s)*fnx+ (h^3*(r - 2)*(19*r + 10*s - 48*r*s + 24*r^3*s - 6*r^4*s + 3*r^2 - 5*r^3 - 9*r^4 + 3*r^5 + 9))/(840*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng -(h^3*(38*g + 3*r + 38*s + 8*g*r - 14*g*s + 8*r*s - 4*g*r^3 - 6*g*r^4 - 4*r^3*s - 6*r^4*s + 3*r^2 + 3*r^3 + 3*r^4 + 3*r^5 + 14*g*r^3*s - 28*g*r*s - 67))/(1680*(g - 2)*(s - 2))*fn2 eq2=-h*y3-(r - 3)/r*y0+ 3/(r*(r - 1))*yr+ (r - 4)/(r - 1)*y1+ (h^3*(88*g + 61*r + 234*g*r - 144*g*r^2 - 144*g*r^3 + 108*g*r^4 - 18*g*r^5 + 39*r^2 + 39*r^3 + 39*r^4 - 45*r^5 + 9*r^6 - 292))/(1680*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(88*g - 27*r + 88*s + 66*g*r + 168*g*s + 66*r*s + 66*g*r^2 + 66*g*r^3 - 102*g*r^4 + 24*g*r^5 + 66*r^2*s + 66*r^3*s - 102*r^4*s + 24*r^5*s - 27*r^2 - 27*r^3 - 27*r^4 + 57*r^5 - 15*r^6 - 210*g*r^2*s + 210*g*r^3*s - 42*g*r^4*s - 210*g*r*s - 292))/(1680*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr+ (h^3*(1712*g + 1745*r + 1712*s - 1684*g*r - 1624*g*s - 1684*r*s + 66*g*r^2 + 66*g*r^3 + 66*g*r^4 - 18*g*r^5 + 66*r^2*s + 66*r^3*s + 66*r^4*s - 18*r^5*s - 27*r^2 - 27*r^3 - 27*r^4 - 27*r^5 + 9*r^6 - 210*g*r^2*s - 210*g*r^3*s + 42*g*r^4*s + 1918*g*r*s - 2004))/(1680*(g - 1)*(r - 1)*(s - 1))*fnxh -(h^3*(88*g + 61*r + 88*s + 234*g*r + 168*g*s + 234*r*s - 144*g*r^2 - 144*g*r^3 + 108*g*r^4 - 18*g*r^5 - 144*r^2*s - 144*r^3*s + 108*r^4*s - 18*r^5*s + 39*r^2 + 39*r^3 + 39*r^4 - 45*r^5 + 9*r^6 + 924*g*r^2*s - 336*g*r^3*s + 42*g*r^4*s - 938*g*r*s - 292))/(3360*g*r*s)*fnx+ (h^3*(61*r + 88*s + 234*r*s - 144*r^2*s - 144*r^3*s + 108*r^4*s - 18*r^5*s + 39*r^2 + 39*r^3 + 39*r^4 - 45*r^5 + 9*r^6 - 292))/(1680*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+(h^3*(920*g + 911*r + 920*s - 486*g*r - 504*g*s - 486*r*s - 24*g*r^2 - 24*g*r^3 - 24*g*r^4 + 18*g*r^5 - 24*r^2*s - 24*r^3*s - 24*r^4*s + 18*r^5*s + 9*r^2 + 9*r^3 + 9*r^4 + 9*r^5 - 9*r^6 + 84*g*r^2*s + 84*g*r^3*s - 42*g*r^4*s + 126*g*r*s - 1548))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq3=-h^2*y4+ 2/r*y0+ 2/(r*(r - 1))*yr -2/(r - 1)*y1+ (h^3*(246*g + 237*r - 48*g*r - 48*g*r^2 - 48*g*r^3 + 36*g*r^4 - 6*g*r^5 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 - 457))/(840*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(246*g - 9*r + 246*s + 22*g*r - 70*g*s + 22*r*s + 22*g*r^2 + 22*g*r^3 - 34*g*r^4 + 8*g*r^5 + 22*r^2*s + 22*r^3*s - 34*r^4*s + 8*r^5*s - 9*r^2 - 9*r^3 - 9*r^4 + 19*r^5 - 5*r^6 - 70*g*r^2*s + 70*g*r^3*s - 14*g*r^4*s - 70*g*r*s - 457))/(840*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr+ (h^3*(1324*g + 1335*r + 1324*s - 1098*g*r - 1078*g*s - 1098*r*s + 22*g*r^2 + 22*g*r^3 + 22*g*r^4 - 6*g*r^5 + 22*r^2*s + 22*r^3*s + 22*r^4*s - 6*r^5*s - 9*r^2 - 9*r^3 - 9*r^4 - 9*r^5 + 3*r^6 - 70*g*r^2*s - 70*g*r^3*s + 14*g*r^4*s + 1050*g*r*s - 1781))/(840*(g - 1)*(r - 1)*(s - 1))*fnxh+(h^3*(48*g*r - 237*r - 246*s - 246*g + 70*g*s + 48*r*s + 48*g*r^2 + 48*g*r^3 - 36*g*r^4 + 6*g*r^5 + 48*r^2*s + 48*r^3*s - 36*r^4*s + 6*r^5*s - 13*r^2 - 13*r^3 - 13*r^4 + 15*r^5 - 3*r^6 - 308*g*r^2*s + 112*g*r^3*s - 14*g*r^4*s + 252*g*r*s + 457))/(1680*g*r*s)*fnx+ (h^3*(237*r + 246*s - 48*r*s - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 - 457))/(840*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+ (h^3*(2022*g + 2019*r + 2022*s - 1128*g*r - 1134*g*s - 1128*r*s - 8*g*r^2 - 8*g*r^3 - 8*g*r^4 + 6*g*r^5 - 8*r^2*s - 8*r^3*s - 8*r^4*s + 6*r^5*s + 3*r^2 + 3*r^3 + 3*r^4 + 3*r^5 - 3*r^6 + 28*g*r^2*s + 28*g*r^3*s - 14*g*r^4*s + 588*g*r*s - 3587))/(1680*(g - 2)*(r - 2)*(s - 2))*fn2 eq4=-h*yn1-(r - 1)/r*y0+ 1/(r*(r - 1))*yr+ (r - 2)/(r - 1)*y1+ (h^3*(r - 1)*(40*g + 27*r - 66*g*r - 18*g*r^2 + 30*g*r^3 - 6*g*r^4 + 14*r^2 + r^3 - 12*r^4 + 3*r^5 - 22))/(1680*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns+ (h^3*(13*r - 40*g - 40*s - 18*g*r + 84*g*s - 18*r*s + 4*g*r^2 + 26*g*r^3 - 8*g*r^4 + 4*r^2*s + 26*r^3*s - 8*r^4*s + 4*r^2 - 5*r^3 - 14*r^4 + 5*r^5 - 56*g*r^2*s + 14*g*r^3*s + 14*g*r*s + 22))/(1680*r*(g - r)*(r - s)*(r - 2))*fnr -(h^3*(72*g + 33*r + 72*s - 60*g*r - 112*g*s - 60*r*s - 38*g*r^2 - 16*g*r^3 + 6*g*r^4 - 38*r^2*s - 16*r^3*s + 6*r^4*s + 24*r^2 + 15*r^3 + 6*r^4 - 3*r^5 + 56*g*r^2*s - 14*g*r^3*s + 126*g*r*s - 50))/(1680*(g - 1)*(s - 1))*fnxh -(h^3*(r - 1)*(40*g + 27*r + 40*s - 66*g*r - 84*g*s - 66*r*s - 18*g*r^2 + 30*g*r^3 - 6*g*r^4 - 18*r^2*s + 30*r^3*s - 6*r^4*s + 14*r^2 + r^3 - 12*r^4 + 3*r^5 - 98*g*r^2*s + 14*g*r^3*s + 210*g*r*s - 22))/(3360*g*r*s)*fnx+ (h^3*(r - 1)*(27*r + 40*s - 66*r*s - 18*r^2*s + 30*r^3*s - 6*r^4*s + 14*r^2 + r^3 - 12*r^4 + 3*r^5 - 22))/(1680*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng -(h^3*(r - 1)*(18*g*r - 9*r - 16*s - 16*g + 28*g*s + 18*r*s + 10*g*r^2 + 2*g*r^3 - 6*g*r^4 + 10*r^2*s + 2*r^3*s - 6*r^4*s - 6*r^2 - 3*r^3 + 3*r^5 - 14*g*r^2*s + 14*g*r^3*s - 42*g*r*s + 10))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq5=-h^2*yn11+ 2/r*y0+ 2/(r*(r - 1))*yr -2/(r - 1)*y1+ (h^3*(162*g*r - 85*r - 76*g - 48*g*r^2 - 48*g*r^3 + 36*g*r^4 - 6*g*r^5 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 47))/(840*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns+ (h^3*(76*g + 9*r + 76*s - 22*g*r - 140*g*s - 22*r*s - 22*g*r^2 - 22*g*r^3 + 34*g*r^4 - 8*g*r^5 - 22*r^2*s - 22*r^3*s + 34*r^4*s - 8*r^5*s + 9*r^2 + 9*r^3 + 9*r^4 - 19*r^5 + 5*r^6 + 70*g*r^2*s - 70*g*r^3*s + 14*g*r^4*s + 70*g*r*s - 47))/(840*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr+ (h^3*(232*g + 243*r + 232*s - 328*g*r - 308*g*s - 328*r*s + 22*g*r^2 + 22*g*r^3 + 22*g*r^4 - 6*g*r^5 + 22*r^2*s + 22*r^3*s + 22*r^4*s - 6*r^5*s - 9*r^2 - 9*r^3 - 9*r^4 - 9*r^5 + 3*r^6 - 70*g*r^2*s - 70*g*r^3*s + 14*g*r^4*s + 490*g*r*s - 185))/(840*(g - 1)*(r - 1)*(s - 1))*fnxh+ (h^3*(76*g + 85*r + 76*s - 162*g*r - 140*g*s - 162*r*s + 48*g*r^2 + 48*g*r^3 - 36*g*r^4 + 6*g*r^5 + 48*r^2*s + 48*r^3*s - 36*r^4*s + 6*r^5*s - 13*r^2 - 13*r^3 - 13*r^4 + 15*r^5 - 3*r^6 - 308*g*r^2*s + 112*g*r^3*s - 14*g*r^4*s + 392*g*r*s - 47))/(1680*g*r*s)*fnx+ (h^3*(162*r*s - 76*s - 85*r - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 47))/(840*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng -(h^3*(36*g + 39*r + 36*s - 62*g*r - 56*g*s - 62*r*s + 8*g*r^2 + 8*g*r^3 + 8*g*r^4 - 6*g*r^5 + 8*r^2*s + 8*r^3*s + 8*r^4*s - 6*r^5*s - 3*r^2 - 3*r^3 - 3*r^4 - 3*r^5 + 3*r^6 - 28*g*r^2*s - 28*g*r^3*s + 14*g*r^4*s + 112*g*r*s - 25))/(1680*(g - 2)*(r - 2)*(s - 2))*fn2 eq6=-h*yr11+ (r - 1)/r*y0+ (2*r - 1)/(r*(r - 1))*yr -r/(r - 1)*y1 -(h^3*r*(r - 1)*(22*g + 4*r - 26*g*r - 74*g*r^2 + 74*g*r^3 - 16*g*r^4 + 17*r^2 + 30*r^3 - 41*r^4 + 10*r^5 - 9))/(1680*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns+(h^3*(22*g - 18*r + 22*s + 44*g*r - 70*g*s + 44*r*s + 66*g*r^2 - 136*g*r^3 + 40*g*r^4 + 66*r^2*s - 136*r^3*s + 40*r^4*s - 27*r^2 - 36*r^3 + 95*r^4 - 30*r^5 + 210*g*r^2*s - 56*g*r^3*s - 140*g*r*s - 9))/(1680*(g - r)*(r - s)*(r - 2))*fnr -(h^3*r*(2*r - 20*g - 20*s + 2*g*r + 42*g*s + 2*r*s + 24*g*r^2 + 46*g*r^3 - 16*g*r^4 + 24*r^2*s + 46*r^3*s - 16*r^4*s - 7*r^2 - 16*r^3 - 25*r^4 + 10*r^5 - 98*g*r^2*s + 28*g*r^3*s - 28*g*r*s + 11))/(1680*(g - 1)*(s - 1))*fnxh+ (h^3*(r - 1)*(22*g + 4*r + 22*s - 26*g*r - 70*g*s - 26*r*s - 74*g*r^2 + 74*g*r^3 - 16*g*r^4 - 74*r^2*s + 74*r^3*s - 16*r^4*s + 17*r^2 + 30*r^3 - 41*r^4 + 10*r^5 - 154*g*r^2*s + 28*g*r^3*s + 238*g*r*s - 9))/(3360*g*s)*fnx -(h^3*r*(r - 1)*(4*r + 22*s - 26*r*s - 74*r^2*s + 74*r^3*s - 16*r^4*s + 17*r^2 + 30*r^3 - 41*r^4 + 10*r^5 - 9))/(1680*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+(h^3*r*(r - 1)*(2*g*r - 6*s - 6*g + 14*g*s + 2*r*s + 10*g*r^2 + 18*g*r^3 - 16*g*r^4 + 10*r^2*s + 18*r^3*s - 16*r^4*s - 3*r^2 - 6*r^3 - 9*r^4 + 10*r^5 - 42*g*r^2*s + 28*g*r^3*s - 14*g*r*s + 3))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq7=-h^2*yr111 + 2/r*y0+ 2/(r*(r - 1))*yr -2/(r - 1)*y1+ (h^3*(22*g + 13*r - 48*g*r - 48*g*r^2 + 232*g*r^3 - 174*g*r^4 + 36*g*r^5 + 13*r^2 + 13*r^3 - 127*r^4 + 111*r^5 - 25*r^6 - 9))/(840*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(22*g - 9*r + 22*s + 22*g*r - 70*g*s + 22*r*s + 22*g*r^2 - 538*g*r^3 + 596*g*r^4 - 160*g*r^5 + 22*r^2*s - 538*r^3*s + 596*r^4*s - 160*r^5*s - 9*r^2 - 9*r^3 + 411*r^4 - 485*r^5 + 135*r^6 + 770*g*r^2*s - 770*g*r^3*s + 196*g*r^4*s - 70*g*r*s - 9))/(840*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr -(h^3*(20*g + 9*r + 20*s - 22*g*r - 42*g*s - 22*r*s - 22*g*r^2 - 22*g*r^3 + 118*g*r^4 - 36*g*r^5 - 22*r^2*s - 22*r^3*s + 118*r^4*s - 36*r^5*s + 9*r^2 + 9*r^3 + 9*r^4 - 75*r^5 + 25*r^6 + 70*g*r^2*s - 210*g*r^3*s + 56*g*r^4*s + 70*g*r*s - 11))/(840*(g - 1)*(r - 1)*(s - 1))*fnxh -(h^3*(22*g + 13*r + 22*s - 48*g*r - 70*g*s - 48*r*s - 48*g*r^2 + 232*g*r^3 - 174*g*r^4 + 36*g*r^5 - 48*r^2*s + 232*r^3*s - 174*r^4*s + 36*r^5*s + 13*r^2 + 13*r^3 - 127*r^4 + 111*r^5 - 25*r^6 - 532*g*r^2*s + 308*g*r^3*s - 56*g*r^4*s + 308*g*r*s - 9))/(1680*g*r*s)*fnx+ (h^3*(13*r + 22*s - 48*r*s - 48*r^2*s + 232*r^3*s - 174*r^4*s + 36*r^5*s + 13*r^2 + 13*r^3 - 127*r^4 + 111*r^5 - 25*r^6 - 9))/(840*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+ (h^3*(6*g + 3*r + 6*s - 8*g*r - 14*g*s - 8*r*s - 8*g*r^2 - 8*g*r^3 + 62*g*r^4 - 36*g*r^5 - 8*r^2*s - 8*r^3*s + 62*r^4*s - 36*r^5*s + 3*r^2 + 3*r^3 + 3*r^4 - 39*r^5 + 25*r^6 + 28*g*r^2*s - 112*g*r^3*s + 56*g*r^4*s + 28*g*r*s - 3))/(1680*(g - 2)*(r - 2)*(s - 2))*fn2 eq8=-ys -((r - s)*(s - 1))/r*y0+ (s*(s - 1))/(r*(r - 1))*yr+ (s*(r - s))/(r - 1)*y1+ (h^3*(22*g + 13*r - 9*s - 15*r^2*s^2 + 3*r^2*s^3 + 3*r^3*s^2 - 48*g*r + 22*g*s + 13*r*s - 48*g*r^2 + 36*g*r^3 - 6*g*r^4 + 22*g*s^2 - 34*g*s^3 + 8*g*s^4 + 13*r*s^2 + 13*r^2*s - 15*r*s^3 - 15*r^3*s + 3*r*s^4 + 3*r^4*s + 13*r^2 + 13*r^3 - 15*r^4 + 3*r^5 - 9*s^2 - 9*s^3 + 19*s^4 - 5*s^5 + 36*g*r*s^2 + 36*g*r^2*s - 6*g*r*s^3 - 6*g*r^3*s - 6*g*r^2*s^2 - 48*g*r*s - 9))/(1680*(g - s)*(s - 2))*fns -(h^3*s*(s - 1)*(22*g - 9*r + 13*s - 15*r^2*s^2 + 3*r^2*s^3 + 3*r^3*s^2 + 22*g*r - 48*g*s + 13*r*s + 22*g*r^2 - 34*g*r^3 + 8*g*r^4 - 48*g*s^2 + 36*g*s^3 - 6*g*s^4 + 13*r*s^2 + 13*r^2*s - 15*r*s^3 - 15*r^3*s + 3*r*s^4 + 3*r^4*s - 9*r^2 - 9*r^3 + 19*r^4 - 5*r^5 + 13*s^2 + 13*s^3 - 15*s^4 + 3*s^5 + 36*g*r*s^2 + 36*g*r^2*s - 6*g*r*s^3 - 6*g*r^3*s - 6*g*r^2*s^2 - 48*g*r*s - 9))/(1680*r*(g - r)*(r^2 - 3*r + 2))*fnr -(h^3*s*(r - s)*(20*g + 9*r + 9*s - 13*r^2*s^2 + 3*r^2*s^3 + 3*r^3*s^2 - 22*g*r - 22*g*s - 13*r*s - 22*g*r^2 - 22*g*r^3 + 6*g*r^4 - 22*g*s^2 - 22*g*s^3 + 6*g*s^4 - 13*r*s^2 - 13*r^2*s - 13*r*s^3 - 13*r^3*s + 3*r*s^4 + 3*r^4*s + 9*r^2 + 9*r^3 + 9*r^4 - 3*r^5 + 9*s^2 + 9*s^3 + 9*s^4 - 3*s^5 + 48*g*r*s^2 + 48*g*r^2*s - 8*g*r*s^3 - 8*g*r^3*s - 8*g*r^2*s^2 + 48*g*r*s - 11))/(1680*(g - 1)*(r - 1))*fnxh -(h^3*(r - s)*(s - 1)^2*(s - 2)*(22*g + 13*r + 13*s + 21*r^2*s^2 - 3*r^2*s^3 - 3*r^3*s^2 - 48*g*r - 48*g*s - 35*r*s - 48*g*r^2 + 36*g*r^3 - 6*g*r^4 - 48*g*s^2 + 36*g*s^3 - 6*g*s^4 - 35*r*s^2 - 35*r^2*s + 21*r*s^3 + 21*r^3*s - 3*r*s^4 - 3*r^4*s + 13*r^2 + 13*r^3 - 15*r^4 + 3*r^5 + 13*s^2 + 13*s^3 - 15*s^4 + 3*s^5 - 76*g*r*s^2 - 76*g*r^2*s + 8*g*r*s^3 + 8*g*r^3*s + 8*g*r^2*s^2 + 260*g*r*s - 9))/(3360*g*r*(s^2 - 3*s + 2))*fnx -(h^3*s*(s - 1)*(- 3*r^6 + 6*r^5*s + 15*r^5 - 36*r^4*s - 13*r^4 + 48*r^3*s - 13*r^3 + 48*r^2*s - 13*r^2 - 6*r*s^5 + 36*r*s^4 - 48*r*s^3 - 48*r*s^2 + 9*r + 3*s^6 - 15*s^5 + 13*s^4 + 13*s^3 + 13*s^2 - 9*s))/(1680*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+ (h^3*s*(r - s)*(s - 1)*(6*g + 3*r + 3*s - 5*r^2*s^2 + 3*r^2*s^3 + 3*r^3*s^2 - 8*g*r - 8*g*s - 5*r*s - 8*g*r^2 - 8*g*r^3 + 6*g*r^4 - 8*g*s^2 - 8*g*s^3 + 6*g*s^4 - 5*r*s^2 - 5*r^2*s - 5*r*s^3 - 5*r^3*s + 3*r*s^4 + 3*r^4*s + 3*r^2 + 3*r^3 + 3*r^4 - 3*r^5 + 3*s^2 + 3*s^3 + 3*s^4 - 3*s^5 + 20*g*r*s^2 + 20*g*r^2*s - 8*g*r*s^3 - 8*g*r^3*s - 8*g*r^2*s^2 + 20*g*r*s - 3))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq9=-h*ys1 -(r - 2*s + 1)/r*y0+ (2*s - 1)/(r*(r - 1))*yr+ (r - 2*s)/(r - 1)*y1+ (h^3*(9*r - 18*s - 22*g*r + 44*g*s + 26*r*s + 48*g*r^2 + 48*g*r^3 - 36*g*r^4 + 6*g*r^5 - 280*g*s^4 + 252*g*s^5 - 56*g*s^6 + 26*r^2*s + 26*r^3*s - 280*r*s^4 + 26*r^4*s + 252*r*s^5 - 30*r^5*s - 56*r*s^6 + 6*r^6*s - 13*r^2 - 13*r^3 - 13*r^4 + 15*r^5 - 3*r^6 + 168*s^5 - 168*s^6 + 40*s^7 - 96*g*r^2*s + 560*g*r*s^3 - 96*g*r^3*s - 420*g*r*s^4 + 72*g*r^4*s + 84*g*r*s^5 - 12*g*r^5*s - 96*g*r*s))/(1680*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns-(h^3*(9*r - 18*s + 44*r^2*s^2 + 44*r^3*s^2 - 68*r^4*s^2 + 16*r^5*s^2 - 22*g*r + 44*g*s - 40*r*s - 22*g*r^2 - 22*g*r^3 + 34*g*r^4 - 8*g*r^5 - 140*g*s^2 + 280*g*s^4 - 168*g*s^5 + 28*g*s^6 + 44*r*s^2 - 40*r^2*s - 40*r^3*s + 16*r^4*s + 30*r^5*s - 10*r^6*s + 9*r^2 + 9*r^3 + 9*r^4 - 19*r^5 + 5*r^6 + 44*s^2 - 112*s^5 + 84*s^6 - 16*s^7 - 140*g*r*s^2 + 114*g*r^2*s - 26*g*r^3*s - 54*g*r^4*s + 16*g*r^5*s - 140*g*r^2*s^2 + 140*g*r^3*s^2 - 28*g*r^4*s^2 + 114*g*r*s))/(1680*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr+ (h^3*(22*s - 11*r + 44*r^2*s^2 + 44*r^3*s^2 + 44*r^4*s^2 - 12*r^5*s^2 + 20*g*r - 40*g*s + 2*r*s - 22*g*r^2 - 22*g*r^3 - 22*g*r^4 + 6*g*r^5 + 84*g*s^2 - 112*g*s^5 + 28*g*s^6 + 44*r*s^2 - 40*r^2*s - 40*r^3*s - 40*r^4*s - 112*r*s^5 - 12*r^5*s + 28*r*s^6 + 6*r^6*s + 9*r^2 + 9*r^3 + 9*r^4 + 9*r^5 - 3*r^6 - 40*s^2 + 56*s^6 - 16*s^7 - 140*g*r*s^2 + 114*g*r^2*s + 114*g*r^3*s + 280*g*r*s^4 + 30*g*r^4*s - 56*g*r*s^5 - 12*g*r^5*s - 140*g*r^2*s^2 - 140*g*r^3*s^2 + 28*g*r^4*s^2 + 2*g*r*s))/(1680*(g - 1)*(r - 1)*(s - 1))*fnxh+ (h^3*(18*s - 9*r + 96*r^2*s^2 + 96*r^3*s^2 - 72*r^4*s^2 + 12*r^5*s^2 + 22*g*r - 44*g*s - 4*r*s - 48*g*r^2 - 48*g*r^3 + 36*g*r^4 - 6*g*r^5 + 140*g*s^2 - 280*g*s^4 + 168*g*s^5 - 28*g*s^6 + 96*r*s^2 - 74*r^2*s - 74*r^3*s - 280*r*s^4 + 10*r^4*s + 168*r*s^5 + 24*r^5*s - 28*r*s^6 - 6*r^6*s + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 - 44*s^2 + 112*s^5 - 84*s^6 + 16*s^7 - 616*g*r*s^2 + 404*g*r^2*s + 1120*g*r*s^3 - 16*g*r^3*s - 420*g*r*s^4 - 58*g*r^4*s + 56*g*r*s^5 + 12*g*r^5*s - 616*g*r^2*s^2 + 224*g*r^3*s^2 - 28*g*r^4*s^2 + 26*g*r*s))/(3360*g*r*s)*fnx-(h^3*(- 6*r^6*s + 3*r^6 + 12*r^5*s^2 + 24*r^5*s - 15*r^5 - 72*r^4*s^2 + 10*r^4*s + 13*r^4 + 96*r^3*s^2 - 74*r^3*s + 13*r^3 + 96*r^2*s^2 - 74*r^2*s + 13*r^2 - 28*r*s^6 + 168*r*s^5 - 280*r*s^4 + 96*r*s^2 - 4*r*s - 9*r + 16*s^7 - 84*s^6 + 112*s^5 - 44*s^2 + 18*s))/(1680*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng -(h^3*(6*s - 3*r + 16*r^2*s^2 + 16*r^3*s^2 + 16*r^4*s^2 - 12*r^5*s^2 + 6*g*r - 12*g*s - 8*g*r^2 - 8*g*r^3 - 8*g*r^4 + 6*g*r^5 + 28*g*s^2 - 56*g*s^5 + 28*g*s^6 + 16*r*s^2 - 14*r^2*s - 14*r^3*s - 14*r^4*s - 56*r*s^5 + 28*r*s^6 + 6*r^6*s + 3*r^2 + 3*r^3 + 3*r^4 + 3*r^5 - 3*r^6 - 12*s^2 + 28*s^6 - 16*s^7 - 56*g*r*s^2 + 44*g*r^2*s + 44*g*r^3*s + 140*g*r*s^4 + 2*g*r^4*s - 56*g*r*s^5 - 12*g*r^5*s - 56*g*r^2*s^2 - 56*g*r^3*s^2 + 28*g*r^4*s^2 + 2*g*r*s))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq10=-h^2*ys11 + 2/r*y0+ 2/(r*(r - 1))*yr -2/(r - 1)*y1+ (h^3*(22*g + 13*r - 48*g*r - 48*g*r^2 - 48*g*r^3 + 36*g*r^4 - 6*g*r^5 - 560*g*s^3 + 630*g*s^4 - 168*g*s^5 - 560*r*s^3 + 630*r*s^4 - 168*r*s^5 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 420*s^4 - 504*s^5 + 140*s^6 + 840*g*r*s^2 - 840*g*r*s^3 + 210*g*r*s^4 - 9))/(840*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(22*g - 9*r + 22*s + 22*g*r - 70*g*s + 22*r*s + 22*g*r^2 + 22*g*r^3 - 34*g*r^4 + 8*g*r^5 + 280*g*s^3 - 210*g*s^4 + 42*g*s^5 + 22*r^2*s + 22*r^3*s - 34*r^4*s + 8*r^5*s - 9*r^2 - 9*r^3 - 9*r^4 + 19*r^5 - 5*r^6 - 140*s^4 + 126*s^5 - 28*s^6 - 70*g*r^2*s + 70*g*r^3*s - 14*g*r^4*s - 70*g*r*s - 9))/(840*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr-(h^3*(20*g + 9*r + 20*s - 22*g*r - 42*g*s - 22*r*s - 22*g*r^2 - 22*g*r^3 - 22*g*r^4 + 6*g*r^5 + 140*g*s^4 - 42*g*s^5 - 22*r^2*s - 22*r^3*s + 140*r*s^4 - 22*r^4*s - 42*r*s^5 + 6*r^5*s + 9*r^2 + 9*r^3 + 9*r^4 + 9*r^5 - 3*r^6 - 84*s^5 + 28*s^6 + 70*g*r^2*s - 280*g*r*s^3 + 70*g*r^3*s + 70*g*r*s^4 - 14*g*r^4*s + 70*g*r*s - 11))/(840*(g - 1)*(r - 1)*(s - 1))*fnxh -(h^3*(22*g + 13*r + 22*s - 48*g*r - 70*g*s - 48*r*s - 48*g*r^2 - 48*g*r^3 + 36*g*r^4 - 6*g*r^5 + 280*g*s^3 - 210*g*s^4 + 42*g*s^5 - 48*r^2*s + 280*r*s^3 - 48*r^3*s - 210*r*s^4 + 36*r^4*s + 42*r*s^5 - 6*r^5*s + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 - 140*s^4 + 126*s^5 - 28*s^6 - 840*g*r*s^2 + 308*g*r^2*s + 420*g*r*s^3 - 112*g*r^3*s - 70*g*r*s^4 + 14*g*r^4*s + 308*g*r*s - 9))/(1680*g*r*s)*fnx+ (h^3*(3*r^6 - 6*r^5*s - 15*r^5 + 36*r^4*s + 13*r^4 - 48*r^3*s + 13*r^3 - 48*r^2*s + 13*r^2 + 42*r*s^5 - 210*r*s^4 + 280*r*s^3 - 48*r*s + 13*r - 28*s^6 + 126*s^5 - 140*s^4 + 22*s - 9))/(840*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+ (h^3*(6*g + 3*r + 6*s - 8*g*r - 14*g*s - 8*r*s - 8*g*r^2 - 8*g*r^3 - 8*g*r^4 + 6*g*r^5 + 70*g*s^4 - 42*g*s^5 - 8*r^2*s - 8*r^3*s + 70*r*s^4 - 8*r^4*s - 42*r*s^5 + 6*r^5*s + 3*r^2 + 3*r^3 + 3*r^4 + 3*r^5 - 3*r^6 - 42*s^5 + 28*s^6 + 28*g*r^2*s - 140*g*r*s^3 + 28*g*r^3*s + 70*g*r*s^4 - 14*g*r^4*s + 28*g*r*s - 3))/(1680*(g - 2)*(r - 2)*(s - 2))*fn2 eq11=-h*y01-(r + 1)/r*y0 -1/(r*(r - 1))*yr+ r/(r - 1)*y1 -(h^3*r*(22*g + 13*r - 48*g*r - 48*g*r^2 + 36*g*r^3 - 6*g*r^4 + 13*r^2 + 13*r^3 - 15*r^4 + 3*r^5 - 9))/(1680*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns+ (h^3*(22*g - 9*r + 22*s + 22*g*r - 70*g*s + 22*r*s + 22*g*r^2 - 34*g*r^3 + 8*g*r^4 + 22*r^2*s - 34*r^3*s + 8*r^4*s - 9*r^2 - 9*r^3 + 19*r^4 - 5*r^5 + 70*g*r^2*s - 14*g*r^3*s - 70*g*r*s - 9))/(1680*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr+ (h^3*r*(20*g + 9*r + 20*s - 22*g*r - 42*g*s - 22*r*s - 22*g*r^2 - 22*g*r^3 + 6*g*r^4 - 22*r^2*s - 22*r^3*s + 6*r^4*s + 9*r^2 + 9*r^3 + 9*r^4 - 3*r^5 + 70*g*r^2*s - 14*g*r^3*s + 70*g*r*s - 11))/(1680*(g - 1)*(r - 1)*(s - 1))*fnxh+ (h^3*(22*g + 13*r + 22*s - 48*g*r - 70*g*s - 48*r*s - 48*g*r^2 + 36*g*r^3 - 6*g*r^4 - 48*r^2*s + 36*r^3*s - 6*r^4*s + 13*r^2 + 13*r^3 - 15*r^4 + 3*r^5 - 112*g*r^2*s + 14*g*r^3*s + 308*g*r*s - 9))/(3360*g*s)*fnx -(h^3*r*(13*r + 22*s - 48*r*s - 48*r^2*s + 36*r^3*s - 6*r^4*s + 13*r^2 + 13*r^3 - 15*r^4 + 3*r^5 - 9))/(1680*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng -(h^3*r*(6*g + 3*r + 6*s - 8*g*r - 14*g*s - 8*r*s - 8*g*r^2 - 8*g*r^3 + 6*g*r^4 - 8*r^2*s - 8*r^3*s + 6*r^4*s + 3*r^2 + 3*r^3 + 3*r^4 - 3*r^5 + 28*g*r^2*s - 14*g*r^3*s + 28*g*r*s - 3))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq12=-h^2*y011 + 2/r*y0+ 2/(r*(r - 1))*yr -2/(r - 1)*y1+ (h^3*(22*g + 13*r - 48*g*r - 48*g*r^2 - 48*g*r^3 + 36*g*r^4 - 6*g*r^5 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 - 9))/(840*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(22*g - 9*r + 22*s + 22*g*r - 70*g*s + 22*r*s + 22*g*r^2 + 22*g*r^3 - 34*g*r^4 + 8*g*r^5 + 22*r^2*s + 22*r^3*s - 34*r^4*s + 8*r^5*s - 9*r^2 - 9*r^3 - 9*r^4 + 19*r^5 - 5*r^6 - 70*g*r^2*s + 70*g*r^3*s - 14*g*r^4*s - 70*g*r*s - 9))/(840*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr -(h^3*(20*g + 9*r + 20*s - 22*g*r - 42*g*s - 22*r*s - 22*g*r^2 - 22*g*r^3 - 22*g*r^4 + 6*g*r^5 - 22*r^2*s - 22*r^3*s - 22*r^4*s + 6*r^5*s + 9*r^2 + 9*r^3 + 9*r^4 + 9*r^5 - 3*r^6 + 70*g*r^2*s + 70*g*r^3*s - 14*g*r^4*s + 70*g*r*s - 11))/(840*(g - 1)*(r - 1)*(s - 1))*fnxh -(h^3*(22*g + 13*r + 22*s - 48*g*r - 70*g*s - 48*r*s - 48*g*r^2 - 48*g*r^3 + 36*g*r^4 - 6*g*r^5 - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 308*g*r^2*s - 112*g*r^3*s + 14*g*r^4*s + 308*g*r*s - 9))/(1680*g*r*s)*fnx+ (h^3*(13*r + 22*s - 48*r*s - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 - 9))/(840*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+ (h^3*(6*g + 3*r + 6*s - 8*g*r - 14*g*s - 8*r*s - 8*g*r^2 - 8*g*r^3 - 8*g*r^4 + 6*g*r^5 - 8*r^2*s - 8*r^3*s - 8*r^4*s + 6*r^5*s + 3*r^2 + 3*r^3 + 3*r^4 + 3*r^5 - 3*r^6 + 28*g*r^2*s + 28*g*r^3*s - 14*g*r^4*s + 28*g*r*s - 3))/(1680*(g - 2)*(r - 2)*(s - 2))*fn2 eq13=-yng + ((g - r)*(g - 1))/r*y0+ (g*(g - 1))/(r*(r - 1))*yr -(g*(g - r))/(r - 1)*y1+ (g*h^3*(g - 1)*(- 3*g^6 + 6*g^5*r + 15*g^5 - 36*g^4*r - 13*g^4 + 48*g^3*r - 13*g^3 + 48*g^2*r - 13*g^2 - 6*g*r^5 + 36*g*r^4 - 48*g*r^3 - 48*g*r^2 + 9*g + 3*r^6 - 15*r^5 + 13*r^4 + 13*r^3 + 13*r^2 - 9*r))/(1680*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns+(g*h^3*(g - 1)*(13*g - 9*r + 22*s + 13*g*r - 48*g*s + 22*r*s + 13*g*r^2 + 13*g^2*r - 15*g*r^3 - 15*g^3*r + 3*g*r^4 + 3*g^4*r - 48*g^2*s + 36*g^3*s - 6*g^4*s + 22*r^2*s - 34*r^3*s + 8*r^4*s + 13*g^2 + 13*g^3 - 15*g^4 + 3*g^5 - 9*r^2 - 9*r^3 + 19*r^4 - 5*r^5 - 15*g^2*r^2 + 3*g^2*r^3 + 3*g^3*r^2 + 36*g*r^2*s + 36*g^2*r*s - 6*g*r^3*s - 6*g^3*r*s - 6*g^2*r^2*s - 48*g*r*s - 9))/(1680*r*(r - s)*(r^2 - 3*r + 2))*fnr+ (g*h^3*(g - r)*(9*g + 9*r + 20*s - 13*g*r - 22*g*s - 22*r*s - 13*g*r^2 - 13*g^2*r - 13*g*r^3 - 13*g^3*r + 3*g*r^4 + 3*g^4*r - 22*g^2*s - 22*g^3*s + 6*g^4*s - 22*r^2*s - 22*r^3*s + 6*r^4*s + 9*g^2 + 9*g^3 + 9*g^4 - 3*g^5 + 9*r^2 + 9*r^3 + 9*r^4 - 3*r^5 - 13*g^2*r^2 + 3*g^2*r^3 + 3*g^3*r^2 + 48*g*r^2*s + 48*g^2*r*s - 8*g*r^3*s - 8*g^3*r*s - 8*g^2*r^2*s + 48*g*r*s - 11))/(1680*(r - 1)*(s - 1))*fnxh+ (h^3*(g - r)*(g - 1)^2*(g - 2)*(13*g + 13*r + 22*s - 35*g*r - 48*g*s - 48*r*s - 35*g*r^2 - 35*g^2*r + 21*g*r^3 + 21*g^3*r - 3*g*r^4 - 3*g^4*r - 48*g^2*s + 36*g^3*s - 6*g^4*s - 48*r^2*s + 36*r^3*s - 6*r^4*s + 13*g^2 + 13*g^3 - 15*g^4 + 3*g^5 + 13*r^2 + 13*r^3 - 15*r^4 + 3*r^5 + 21*g^2*r^2 - 3*g^2*r^3 - 3*g^3*r^2 - 76*g*r^2*s - 76*g^2*r*s + 8*g*r^3*s + 8*g^3*r*s + 8*g^2*r^2*s + 260*g*r*s - 9))/(3360*r*s*(g^2 - 3*g + 2))*fnx -(h^3*(13*r - 9*g + 22*s + 13*g*r + 22*g*s - 48*r*s + 13*g*r^2 + 13*g^2*r - 15*g*r^3 - 15*g^3*r + 3*g*r^4 + 3*g^4*r + 22*g^2*s - 34*g^3*s + 8*g^4*s - 48*r^2*s + 36*r^3*s - 6*r^4*s - 9*g^2 - 9*g^3 + 19*g^4 - 5*g^5 + 13*r^2 + 13*r^3 - 15*r^4 + 3*r^5 - 15*g^2*r^2 + 3*g^2*r^3 + 3*g^3*r^2 + 36*g*r^2*s + 36*g^2*r*s - 6*g*r^3*s - 6*g^3*r*s - 6*g^2*r^2*s - 48*g*r*s - 9))/(1680*(g - s)*(g - 2))*fng -(g*h^3*(g - r)*(g - 1)*(3*g + 3*r + 6*s - 5*g*r - 8*g*s - 8*r*s - 5*g*r^2 - 5*g^2*r - 5*g*r^3 - 5*g^3*r + 3*g*r^4 + 3*g^4*r - 8*g^2*s - 8*g^3*s + 6*g^4*s - 8*r^2*s - 8*r^3*s + 6*r^4*s + 3*g^2 + 3*g^3 + 3*g^4 - 3*g^5 + 3*r^2 + 3*r^3 + 3*r^4 - 3*r^5 - 5*g^2*r^2 + 3*g^2*r^3 + 3*g^3*r^2 + 20*g*r^2*s + 20*g^2*r*s - 8*g*r^3*s - 8*g^3*r*s - 8*g^2*r^2*s + 20*g*r*s - 3))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq14=-h*yng1 -(r - 2*g + 1)/r*y0+ (2*g - 1)/(r*(r - 1))*yr -(2*g - r)/(r - 1)*y1 -(h^3*(16*g^7 - 28*g^6*r - 84*g^6 + 168*g^5*r + 112*g^5 - 280*g^4*r + 12*g^2*r^5 - 72*g^2*r^4 + 96*g^2*r^3 + 96*g^2*r^2 + 96*g^2*r - 44*g^2 - 6*g*r^6 + 24*g*r^5 + 10*g*r^4 - 74*g*r^3 - 74*g*r^2 - 4*g*r + 18*g + 3*r^6 - 15*r^5 + 13*r^4 + 13*r^3 + 13*r^2 - 9*r))/(1680*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(9*r - 18*g - 40*g*r + 44*g*s - 22*r*s - 40*g*r^2 + 44*g^2*r - 40*g*r^3 + 16*g*r^4 + 30*g*r^5 - 10*g*r^6 - 140*g^2*s + 280*g^4*s - 168*g^5*s + 28*g^6*s - 22*r^2*s - 22*r^3*s + 34*r^4*s - 8*r^5*s + 44*g^2 - 112*g^5 + 84*g^6 - 16*g^7 + 9*r^2 + 9*r^3 + 9*r^4 - 19*r^5 + 5*r^6 + 44*g^2*r^2 + 44*g^2*r^3 - 68*g^2*r^4 + 16*g^2*r^5 + 114*g*r^2*s - 140*g^2*r*s - 26*g*r^3*s - 54*g*r^4*s + 16*g*r^5*s - 140*g^2*r^2*s + 140*g^2*r^3*s - 28*g^2*r^4*s + 114*g*r*s))/(1680*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr+ (h^3*(22*g - 11*r + 2*g*r - 40*g*s + 20*r*s - 40*g*r^2 + 44*g^2*r - 40*g*r^3 - 40*g*r^4 - 12*g*r^5 - 112*g^5*r + 6*g*r^6 + 28*g^6*r + 84*g^2*s - 112*g^5*s + 28*g^6*s - 22*r^2*s - 22*r^3*s - 22*r^4*s + 6*r^5*s - 40*g^2 + 56*g^6 - 16*g^7 + 9*r^2 + 9*r^3 + 9*r^4 + 9*r^5 - 3*r^6 + 44*g^2*r^2 + 44*g^2*r^3 + 44*g^2*r^4 - 12*g^2*r^5 + 114*g*r^2*s - 140*g^2*r*s + 114*g*r^3*s + 30*g*r^4*s + 280*g^4*r*s - 12*g*r^5*s - 56*g^5*r*s - 140*g^2*r^2*s - 140*g^2*r^3*s + 28*g^2*r^4*s + 2*g*r*s))/(1680*(g - 1)*(r - 1)*(s - 1))*fnxh+ (h^3*(18*g - 9*r - 4*g*r - 44*g*s + 22*r*s - 74*g*r^2 + 96*g^2*r - 74*g*r^3 + 10*g*r^4 - 280*g^4*r + 24*g*r^5 + 168*g^5*r - 6*g*r^6 - 28*g^6*r + 140*g^2*s - 280*g^4*s + 168*g^5*s - 28*g^6*s - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s - 44*g^2 + 112*g^5 - 84*g^6 + 16*g^7 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 96*g^2*r^2 + 96*g^2*r^3 - 72*g^2*r^4 + 12*g^2*r^5 + 404*g*r^2*s - 616*g^2*r*s - 16*g*r^3*s + 1120*g^3*r*s - 58*g*r^4*s - 420*g^4*r*s + 12*g*r^5*s + 56*g^5*r*s - 616*g^2*r^2*s + 224*g^2*r^3*s - 28*g^2*r^4*s + 26*g*r*s))/(3360*g*r*s)*fnx -(h^3*(18*g - 9*r - 26*g*r - 44*g*s + 22*r*s - 26*g*r^2 - 26*g*r^3 - 26*g*r^4 + 280*g^4*r + 30*g*r^5 - 252*g^5*r - 6*g*r^6 + 56*g^6*r + 280*g^4*s - 252*g^5*s + 56*g^6*s - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s - 168*g^5 + 168*g^6 - 40*g^7 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 96*g*r^2*s + 96*g*r^3*s - 560*g^3*r*s - 72*g*r^4*s + 420*g^4*r*s + 12*g*r^5*s - 84*g^5*r*s + 96*g*r*s))/(1680*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng -(h^3*(6*g - 3*r - 12*g*s + 6*r*s - 14*g*r^2 + 16*g^2*r - 14*g*r^3 - 14*g*r^4 - 56*g^5*r + 6*g*r^6 + 28*g^6*r + 28*g^2*s - 56*g^5*s + 28*g^6*s - 8*r^2*s - 8*r^3*s - 8*r^4*s + 6*r^5*s - 12*g^2 + 28*g^6 - 16*g^7 + 3*r^2 + 3*r^3 + 3*r^4 + 3*r^5 - 3*r^6 + 16*g^2*r^2 + 16*g^2*r^3 + 16*g^2*r^4 - 12*g^2*r^5 + 44*g*r^2*s - 56*g^2*r*s + 44*g*r^3*s + 2*g*r^4*s + 140*g^4*r*s - 12*g*r^5*s - 56*g^5*r*s - 56*g^2*r^2*s - 56*g^2*r^3*s + 28*g^2*r^4*s + 2*g*r*s))/(3360*(g - 2)*(r - 2)*(s - 2))*fn2 eq15=-h^2*yng11+ 2/r*y0+ 2/(r*(r - 1))*yr -2/(r - 1)*y1+ (h^3*(- 28*g^6 + 42*g^5*r + 126*g^5 - 210*g^4*r - 140*g^4 + 280*g^3*r - 6*g*r^5 + 36*g*r^4 - 48*g*r^3 - 48*g*r^2 - 48*g*r + 22*g + 3*r^6 - 15*r^5 + 13*r^4 + 13*r^3 + 13*r^2 + 13*r - 9))/(840*s*(g - s)*(r - s)*(s^2 - 3*s + 2))*fns -(h^3*(22*g - 9*r + 22*s + 22*g*r - 70*g*s + 22*r*s + 22*g*r^2 + 22*g*r^3 - 34*g*r^4 + 8*g*r^5 + 280*g^3*s - 210*g^4*s + 42*g^5*s + 22*r^2*s + 22*r^3*s - 34*r^4*s + 8*r^5*s - 140*g^4 + 126*g^5 - 28*g^6 - 9*r^2 - 9*r^3 - 9*r^4 + 19*r^5 - 5*r^6 - 70*g*r^2*s + 70*g*r^3*s - 14*g*r^4*s - 70*g*r*s - 9))/(840*r*(g - r)*(r - s)*(r^2 - 3*r + 2))*fnr -(h^3*(20*g + 9*r + 20*s - 22*g*r - 42*g*s - 22*r*s - 22*g*r^2 - 22*g*r^3 - 22*g*r^4 + 140*g^4*r + 6*g*r^5 - 42*g^5*r + 140*g^4*s - 42*g^5*s - 22*r^2*s - 22*r^3*s - 22*r^4*s + 6*r^5*s - 84*g^5 + 28*g^6 + 9*r^2 + 9*r^3 + 9*r^4 + 9*r^5 - 3*r^6 + 70*g*r^2*s + 70*g*r^3*s - 280*g^3*r*s - 14*g*r^4*s + 70*g^4*r*s + 70*g*r*s - 11))/(840*(g - 1)*(r - 1)*(s - 1))*fnxh -(h^3*(22*g + 13*r + 22*s - 48*g*r - 70*g*s - 48*r*s - 48*g*r^2 - 48*g*r^3 + 280*g^3*r + 36*g*r^4 - 210*g^4*r - 6*g*r^5 + 42*g^5*r + 280*g^3*s - 210*g^4*s + 42*g^5*s - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s - 140*g^4 + 126*g^5 - 28*g^6 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 308*g*r^2*s - 840*g^2*r*s - 112*g*r^3*s + 420*g^3*r*s + 14*g*r^4*s - 70*g^4*r*s + 308*g*r*s - 9))/(1680*g*r*s)*fnx+ (h^3*(13*r + 22*s - 48*r*s - 560*g^3*r + 630*g^4*r - 168*g^5*r - 560*g^3*s + 630*g^4*s - 168*g^5*s - 48*r^2*s - 48*r^3*s + 36*r^4*s - 6*r^5*s + 420*g^4 - 504*g^5 + 140*g^6 + 13*r^2 + 13*r^3 + 13*r^4 - 15*r^5 + 3*r^6 + 840*g^2*r*s - 840*g^3*r*s + 210*g^4*r*s - 9))/(840*g*(g - r)*(g - s)*(g^2 - 3*g + 2))*fng+ (h^3*(6*g + 3*r + 6*s - 8*g*r - 14*g*s - 8*r*s - 8*g*r^2 - 8*g*r^3 - 8*g*r^4 + 70*g^4*r + 6*g*r^5 - 42*g^5*r + 70*g^4*s - 42*g^5*s - 8*r^2*s - 8*r^3*s - 8*r^4*s + 6*r^5*s - 42*g^5 + 28*g^6 + 3*r^2 + 3*r^3 + 3*r^4 + 3*r^5 - 3*r^6 + 28*g*r^2*s + 28*g*r^3*s - 140*g^3*r*s - 14*g*r^4*s + 70*g^4*r*s + 28*g*r*s - 3))/(1680*(g - 2)*(r - 2)*(s - 2))*fn2 [y1,y2,y3,y4,yn1,yn11,yr,yr11,yr111,ys,ys1,ys11,yng,yng1,yng11]=vpasolve([eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10,eq11,eq12,eq13,eq14,eq15])
  5 Comments

Sign in to comment.

Answers (1)

Alex Sha
Alex Sha on 10 Dec 2019
Hi, it is easy to get the results below:
y1: 4.99999995833333E-5
y2: 0.000199999986666668
y3: 0.0199999966666674
y4: 0.999999333333578
yn1: 0.00999999979166667
yn11: 0.99999991666667
yr: 0.000172885992950938
yr11: 0.0185949431780768
yr111: 0.999999464199188
ys: 0.000116916559434617
ys1: 0.015291601976209
ys11: 0.999999702026436
yng: 8.77296003106683E-6
yng1: 0.00418878502899688
yng11: 0.999999993875326

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!