求助：matlab中如何定义分段函数、定积分和不定积分

我想定义一个分段函数2 X' K8 [. q* B4 N
P=0 (u<=4);- l" ~  p1 h2 P% |8 M( C
P=166.67*u-666.67 (u>4 &u<=13)& v3 p# G7 a& j: k3 S* \
P=1500 (u>13 & u<=20)
7 o  z& r9 R  }$ PP=0 (u>20)

我在matlab中定义的函数如下

function [ P ] = Power_curve( u )' Z* v( I# M7 w4 q" |
P=eps.*(u<=4)+(166.67*u-666.67).*(u>4 & u<=13)+1500.*(u>13 & u<=20+eps.*(u>20));
" z" m0 J+ U1 Bend

之后我想求函数0到25区间内的定积分

syms u;" B# x8 O% {, _! W! N* X, j
int(Power_curve(u),1,2)

然后被报错了

Error using symengine
. L) r* }6 t# k7 H+ q1 j' k: cThe integrand is invalid.

Error in sym/int (line 151)  X) {5 ?1 L% }) _3 ]
rSym = mupadmex('symobj::intdef',f.s,x.s,a.s,b.s,options);

Error in eg1 (line 3)
6 ~; a. L2 x3 {8 c7 w: |+ \" r2 `- nint(Power_curve(u),0,25)

请问大神这怎么解决？