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

我想定义一个分段函数4 E8 _) w6 X: S$ O5 L! ]
P=0 (u<=4);4 \$ d- W7 q- ?) n) Z
P=166.67*u-666.67 (u>4 &u<=13)" S9 \" t6 h. ^
P=1500 (u>13 & u<=20)
6 Q, P. p3 e, T! j7 [P=0 (u>20)

我在matlab中定义的函数如下

function [ P ] = Power_curve( u )% c# f" A! }: V
P=eps.*(u<=4)+(166.67*u-666.67).*(u>4 & u<=13)+1500.*(u>13 & u<=20+eps.*(u>20));
/ L$ B" c2 g$ b- ]# F2 q6 s% m$ |end

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

syms u;& x6 z5 F& v( X- t" ], r
int(Power_curve(u),1,2)

然后被报错了

Error using symengine
* s' X; J5 b& S7 rThe integrand is invalid.

Error in sym/int (line 151)
" N* i5 a4 S) q/ x5 z' JrSym = mupadmex('symobj::intdef',f.s,x.s,a.s,b.s,options);

Error in eg1 (line 3)
" m+ Y. j6 i% q5 H9 N$ z0 F0 R% jint(Power_curve(u),0,25)

请问大神这怎么解决？