Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- S4 :=
- Function[{x, w, y1, y2},
- Piecewise[{{(11*(y1 + y2))/24,
- w == 1/2 + x && w + x == 1/2 && 4*w > 1}, {(y1 + y2)/2,
- Element[x, Reals] && w == 0}, {(2*(y1 + y2))/3,
- w == 1/2 + x && w + x == 1/2 && 0 < w < 1/4},
- {(1/120)*(4*y1 + y2),
- 4*w == 1 && x == -(3/4)}, {((-(11 + 20*x))*y1 + y2)/(480*w),
- w + x == -(1/2) && w >= 1 + x && 4*w > 1},
- {(11/120)*(3*y1 + 2*y2), (w == 1/2 + x && w + x == 0 &&
- 4*w > 1) || (4*w == 1 && x == -(1/4))}, {(11/120)*(2*y1 +
- 3*y2), (w == x && w + x == 1/2 && 4*w > 1) || (4*w == 1 &&
- 4*x == 1)},
- {(1/120)*(y1 + 4*y2),
- 4*w == 1 && 4*x == 3}, {(2/15)*(3*y1 + 7*y2),
- w == x && w + x == 1 &&
- 0 < w < 1/4}, {-((7*y1 + 60*x*y1 - 7*y2)/(120*w)),
- w + x == 0 && w >= 1 + x},
- {(1/30)*(-57*y1 + 92*y2),
- w == 1/2 + x && w + x == 1 &&
- 0 < w < 1/4}, {(109*y1 - 460*x*y1 + 121*y2)/(480*w),
- w + x == 1/2 && w >= 1 + x},
- {(1/90)*(-181*y1 + 286*y2),
- w == 1/2 + x && w + x == 1 &&
- 0 < w < 1/4}, {-(((5 - 4*x)^4*((-5 + 4*x)*y1 - (5 + 4*x)*y2))/
- 3840), 4*w == 1 && 3/4 < x < 5/4},
- {-(((5 + 4*x)^4*((-5 + 4*x)*y1 - (5 + 4*x)*y2))/3840),
- 4*w == 1 && -(5/4) <
- x < -(3/4)}, {(1/960)*(-55 +
- 8*x*(-5 + 4*x*(15 + 4*x*(-5 + 2*x))))*((-5 + 4*x)*
- y1 - (5 + 4*x)*y2),
- 4*w == 1 &&
- 1/4 < x <
- 3/4}, {(1/960)*(-55 +
- 8*x*(5 + 4*x*(15 + 4*x*(5 + 2*x))))*((-5 + 4*x)*
- y1 - (5 + 4*x)*y2), 4*w == 1 && -(3/4) < x < -(1/4)},
- {-(((115 +
- 96*x^2*(-5 + 8*x^2))*((-5 + 4*x)*y1 - (5 + 4*x)*y2))/
- 1920), 4*w == 1 && -(1/4) < x <
- 1/4}, {(x*(-y1 + y2) + w*(y1 + y2))/(2*w),
- w >= 1 + x && w + x >= 1},
- {((1 + w + x)^4*((-1 + 9*w - x)*y1 + (1 + w + x)*y2))/(15*w),
- 1 + w + x >
- 0 && ((4*w > 1 && 1/2 + w + x < 0) || (w >= 1 + x &&
- 4*w < 1))},
- {((1 + w - x)^4*((1 + w - x)*y1 + (-1 + 9*w + x)*y2))/(15*
- w), (4*w < 1 && 1 + w > x && w + x > 1) || (4*w > 1 &&
- x < 1 + w && 1/2 + w < x)},
- {(y1 - 11*y2 + 20*x*y2)/(480*w),
- 1/2 + w == x && w + x >= 1 &&
- 4*w > 1}, {(7*y1 - 7*y2 + 60*x*y2)/(120*x),
- w == x && w + x >= 1},
- {(121*y1 + 109*y2 + 460*x*y2)/(480*w),
- w == 1/2 + x && w + x >= 1}, {(8/15)*
- w*(-5*w*(-1 + x)^2*(y1 - y2) + w^3*(-y1 + y2) -
- 5*w^2*(-1 + x)*(y1 + y2) - 5*(-1 + x)^3*(y1 + y2)),
- 0 < w <
- 1/4 && ((w + x == 1 &&
- 1/2 + w <= x) || (Inequality[1/2, Less, w + x, LessEqual,
- 1] && 1/2 + w < x))},
- {(8/15)*
- w*(w^3*(y1 - y2) + 5*w*(1 + x)^2*(y1 - y2) +
- 5*w^2*(1 + x)*(y1 + y2) + 5*(1 + x)^3*(y1 + y2)),
- w > 0 && w < 1 + x && w + x <= -(1/2)},
- {(4/15)*
- x*(5*y1 + 5*y2 - 20*x^2*(3*y1 + y2) + 24*x^3*(4*y1 + y2)),
- w == x && 4*w < 1 && Inequality[0, Less, w + x, LessEqual, 1/2]},
- {(4/15)*
- w*(6*w^3*(y1 - y2) + 10*w*x*(-2 + 3*x)*(y1 - y2) +
- 10*w^2*(-1 + 3*x)*(y1 + y2) +
- 5*(1 + 6*(-1 + x)*x^2)*(y1 + y2)),
- 0 < w < 1/4 &&
- x <= 1/2 +
- w && ((w + x == 1/2 &&
- w <= x) || (Inequality[0, Less, w + x, LessEqual, 1/2] &&
- w < x))},
- {(4/15)*
- w*(-10*w*x*(2 + 3*x)*(y1 - y2) + 6*w^3*(-y1 + y2) -
- 10*w^2*(1 + 3*x)*(y1 + y2) -
- 5*(-1 + 6*x^2*(1 + x))*(y1 + y2)),
- 0 < w < 1/4 &&
- x <= w && ((w + x == 0 &&
- w <= 1/2 + x) || (Inequality[-(1/2), Less, w + x,
- LessEqual, 0] && w < 1/2 + x))},
- {-(((1 - 2*x)^2*((-71 + 8*x*(32 + x*(-39 + 16*x)))*
- y1 + (-29 + 4*x*(21 - 22*x + 8*x^2))*y2))/(60*w)),
- w + x > 1/2 && 1/2 + w == x && 4*w < 1},
- {-(((1 + 2*x)^2*((-7 + 8*x*(16 + 57*x + 48*x^2))*
- y1 + (-13 + 4*x*(3 + 26*x + 24*x^2))*y2))/(60*w)),
- w == 1/2 + x && 0 < w < 1/4 && w + x <= 0},
- {((55 + 4*x*(-55 + 4*x*(5 + 4*x*(5 + 2*(5 - 12*x)*x))))*
- y1 + (55 - 16*x*(5 + x*(25 + 8*x*(5 - 45*x + 48*x^2))))*
- y2)/(240*w), w + x == 1/2 && w < 1/2 + x && 4*w > 1},
- {((55 + 16*x*(5 + x*(-25 + 8*x*(5 + (5 - 16*x)*x))))*
- y1 + (55 + 4*x*(55 + 4*x*(5 - 4*x*(5 + 2*x*(-5 + 4*x)))))*
- y2)/(240*w), w == 1/2 + x && 1/2 < w + x < 1},
- {((-1 - 64*x^2*(-5 + 2*x*(10 + x*(-15 + 8*x))))*y1 + y2 -
- 4*x*(5 + 16*x*(-5 + 2*x*(5 + x*(-5 + 2*x))))*y2)/(120*x),
- w == x && 1/2 < w + x < 1},
- {-((1/(120*
- w))*((224*w^5 + 10*w*(1 + 2*x)^4 - (1 + 2*x)^5 +
- 80*w^4*(3 + 10*x) + 80*w^3*(1 + 4*x*(3 + 4*x)) +
- 40*w^2*(-5 + 2*x*(-3 + 6*x + 8*x^2)))*y1 +
- (96*w^5 + 10*w*(1 + 2*x)^4 + (1 + 2*x)^5 -
- 80*w^4*(3 + 2*x) + 80*w^3*(5 + 4*x*(3 + 2*x)) -
- 40*w^2*(7 + 6*x*(3 + 2*x)))*y2)),
- (w < 1 + x && w + x < 0 && 4*w > 1) || (w > 1/2 + x &&
- w + x > -(1/2) && 4*w < 1)},
- {-((1/(120*
- w))*((96*w^5 + 10*w*(1 - 2*x)^4 +
- 80*w^4*(-3 + 2*x) - (-1 + 2*x)^5 +
- 80*w^3*(5 + 4*x*(-3 + 2*x)) -
- 40*w^2*(7 + 6*x*(-3 + 2*x)))*y1 +
- (224*w^5 + 80*w^4*(3 - 10*x) +
- 10*w*(1 - 2*x)^4 + (-1 + 2*x)^5 +
- 80*w^3*(1 + 4*x*(-3 + 4*x)) -
- 40*w^2*(5 + 2*x*(-3 - 6*x + 8*x^2)))*y2)),
- (4*w < 1 && 1/2 + w > x && w + x > 1/2) || (4*w > 1 &&
- w < x && w + x < 1)},
- {-((1/(120*
- w))*((7 + 216*w^5 + 280*w^4*(1 + 3*x) +
- 400*w^3*x*(2 + 3*x) +
- 10*w*(-1 + 2*x^2)*(3 + 8*x + 6*x^2) +
- 120*w^2*(-1 + 6*x^2*(1 + x)) +
- 2*x*(15 + 20*x - 4*x^3*(5 + 3*x)))*y1 - 7*y2 +
- 2*(w + x)*(-5*(3 + 4*x) +
- 4*(3*w^4 + x^3*(5 + 3*x) + 3*w^2*x*(5 + 6*x) +
- w^3*(5 + 12*x) + w*(-5 + 3*x^2*(5 + 4*x))))*y2)), -(1/
- 2) < w + x < 0 && w >= 1 + x},
- {(1/(120*
- w))*((7 + 2*w*(15 + 20*w - 4*w^3*(5 + 3*w)) - 30*x +
- 40*w*(-2 + w^2*(4 + 3*w))*x -
- 40*(-1 + 6*w^2*(1 + w))*x^2 + 80*w*(2 + 3*w)*x^3 -
- 40*(1 + 3*w)*x^4 + 24*x^5)*y1 +
- (-7 + 2*w*(15 + 60*w - 4*w^3*(35 + 27*w)) + 30*x +
- 40*w*(-2 + w^2*(20 + 21*w))*x -
- 40*(1 + 6*w^2*(3 + 5*w))*x^2 + 80*w*(2 + 9*w)*x^3 +
- 40*(1 - 3*w)*x^4 - 24*x^5)*y2),
- x < 1/2 + w && w < x &&
- w + x >=
- 1}, {(1/(30*w))*(x*(5 + 40*x^2 + 4*x^4)*(y1 - y2) -
- 20*w^5*(y1 + y2) - 5*w*(1 + 4*x^2*(6 + x^2))*(y1 + y2) +
- 20*w^4*((4 - 3*x)*y1 + (4 + 3*x)*y2) +
- 40*w^2*(2*y1 - x*(3 + (-6 + x)*x)*y1 + 2*y2 +
- x*(3 + x*(6 + x))*y2) -
- 40*w^3*((3 + x*(-4 + 3*x))*y1 + (3 + x*(4 + 3*x))*y2)),
- 1/2 < w - x < 1 && 1/2 < w + x < 1},
- {(1/(30*
- w))*((2 - 2*w^5 + 10*w^4*(1 + x) - 20*w^3*(1 + x)^2 +
- 20*w^2*(1 + x)^3 -
- 5*w*(-1 + 2*x*(2 + x)*(2 + x*(2 + x))) +
- x*(-5 + 2*x*(10 + x*(10 + x*(5 + x)))))*y1 -
- (2 + 18*w^5 - 70*w^4*(1 + x) + 100*w^3*(1 + x)^2 -
- 60*w^2*(1 + x)^3 +
- 5*w*(-1 + 2*x*(2 + x)*(2 + x*(2 + x))) +
- x*(-5 + 2*x*(10 + x*(10 + x*(5 + x)))))*y2),
- 1/2 + x < w && w < 1 + x && w + x >= 1},
- {-((1/(30*
- w))*((2 + 18*w^5 + 70*w^4*(-1 + x) + 100*w^3*(-1 + x)^2 +
- 60*w^2*(-1 + x)^3 +
- 5*w*(-1 + 2*(-2 + x)*x*(2 + (-2 + x)*x)) +
- x*(5 - 2*x*(-10 + x*(10 + (-5 + x)*x))))*y1 - 2*y2 +
- (w + x)*(-5*(1 + 4*x) +
- 2*(w*(-10 + w*(10 + (-5 + w)*w)) +
- w*(20 + w*(-15 + 4*w))*x + (10 + 3*w*(-5 + 2*w))*
- x^2 + (-5 + 4*w)*x^3 + x^4))*y2)),
- 1/2 < w + x < 1 && w >= 1 + x},
- {(2*x^5*(-y1 + y2))/(5*w) + 2*w^4*(y1 + y2) +
- 2*x^4*(y1 + y2) + (4/3)*
- w*(y1 + 3*(-2 + x)*x^2*y1 + y2 - 3*x^2*(2 + x)*y2) + (2/3)*
- w^3*((-4 + 9*x)*y1 - (4 + 9*x)*y2) +
- (4/3)*w^2*
- x*((-4 + 9*x)*y1 + (4 + 9*x)*y2), (x < w && w + x > 0 &&
- 4*w < 1) || (w < 1/2 + x && w + x < 1/2 && 4*w > 1)},
- {(1/(120*
- w))*((1 + 208*w^5 + 80*w^4*(-3 + 11*x) +
- 80*w^3*(-1 + 2*x*(-6 + 7*x)) +
- 40*w^2*(5 + 2*x*(3 + 2*x*(-3 + 5*x))) +
- 2*x*(5 + 4*x*(5 - 2*x*(-5 + (-5 + x)*x))) +
- 10*w*(-1 + 8*x*(-1 + x*(-3 + (-4 + x)*x))))*
- y1 - (1 + 48*w^5 - 80*w^4*(3 + 5*x) +
- 80*w^3*(5 + 2*x*(6 + x)) -
- 40*w^2*(7 + 6*x*(3 + 2*x*(1 + x))) +
- 2*x*(5 + 4*x*(5 - 2*x*(-5 + (-5 + x)*x))) +
- 10*w*(1 + 8*x*(1 + x*(3 - (-4 + x)*x))))*y2),
- 0 < w + x < 1/2 && 1/2 < w - x < 1},
- {(1/(120*
- w))*((-(1 + 48*w^5 + 80*w^4*(-3 + 5*x) +
- 80*w^3*(5 + 2*(-6 + x)*x) +
- 40*w^2*(-7 + 6*x*(3 + 2*(-1 + x)*x)) -
- 10*w*(-1 + 8*x*(1 + x*(-3 + x*(4 + x)))) +
- 2*x*(-5 + 4*x*(5 + 2*x*(-5 + x*(5 + x))))))*
- y1 + y2 +
- 2*(104*w^5 - 40*w^4*(3 + 11*x) +
- 40*w^3*(-1 + 2*x*(6 + 7*x)) -
- 20*w^2*(-5 + 2*x*(3 + 2*x*(3 + 5*x))) +
- 5*w*(-1 + 8*x*(1 + x*(-3 + x*(4 + x)))) +
- x*(-5 + 4*x*(5 + 2*x*(-5 + x*(5 + x)))))*y2),
- x < w && w < 1/2 + x && 1/2 < w + x < 1},
- {(1/(120*
- w))*((-7 + 216*w^5 + 400*w^3*x*(-2 + 3*x) +
- 280*w^4*(-1 + 3*x) + 120*w^2*(1 + 6*(-1 + x)*x^2) +
- 10*w*(3 + 8*x + 4*x^3*(-4 + 3*x)) -
- 2*x*(15 + 4*x*(5 + x^2*(-5 + 3*x))))*y1 + 7*y2 +
- 2*(w + x)*(5*(3 + 4*x) +
- 4*(3*w^4 + x^3*(-5 + 3*x) + 3*w^2*x*(-5 + 6*x) +
- w^3*(-5 + 12*x) + w*(5 + 3*x^2*(-5 + 4*x))))*y2),
- 0 < w + x < 1/2 && w >= 1 + x},
- {(1/(120*
- w))*((7 + 24*w^5 - 40*w^4*(1 + 3*x) + 80*w^3*x*(2 + 3*x) -
- 40*w^2*(-1 + 6*x^2*(1 + x)) +
- 10*w*(3 - 8*x + 4*x^3*(4 + 3*x)) -
- 2*x*(15 + 4*x*(-5 + x^2*(5 + 3*x))))*y1 +
- (-7 + 216*w^5 - 280*w^4*(1 + 3*x) +
- 400*w^3*x*(2 + 3*x) - 120*w^2*(-1 + 6*x^2*(1 + x)) +
- 10*w*(3 - 8*x + 4*x^3*(4 + 3*x)) +
- 2*x*(15 + 4*x*(-5 + x^2*(5 + 3*x))))*y2),
- x < w && w < 1/2 + x &&
- w + x >= 1}, {(y1 +
- 4*x*(5 + 16*x*(5 + 2*x*(5 + x*(5 + 2*x))))*
- y1 + (-1 + 64*x^2*(5 + 2*x*(10 + x*(15 + 8*x))))*y2)/(120*
- w), w + x == 0 && w < 1 + x && 4*w > 1},
- {((55 + 16*x*(5 + x*(-25 + 8*x*(5 + 45*x + 48*x^2))))*
- y1 + (55 + 4*x*(55 + 4*x*(5 + 4*x*(-5 + 2*x*(5 + 12*x)))))*
- y2)/(240*w), w == 1/2 + x && 0 < w + x < 1/2},
- {((55 + 4*x*(-55 + 4*x*(5 + 4*x*(5 + 2*x*(5 + 4*x)))))*
- y1 + (55 + 16*x*(-5 + x*(-25 + 8*x*(-5 + x*(5 + 16*x)))))*
- y2)/(240*w), w + x == 1/2 && 1/2 + x < w && w < 1 + x}}, 0]]
- In[241]:= piecewiseFunction = S4[x, w, y1, y2]
- In[239]:= oneOfTheConditions = w == 1/2 + x && w + x == 1/2 && 4 w > 1
- Out[239]= w == 1/2 + x && w + x == 1/2 && 4 w > 1
- In[240]:= Reduce[oneOfTheConditions, {x, w}]
- Out[240]= x == 0 && w == 1/2
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement