您可以给变量数值
kill(nmax, dp, offset, c1);
numerval(nmax, 1);
numerval(dp, 1);
numerval(offset, 1);
/* fill a list with dummy values */
c1: makelist(1/m+m*%i, m, 1, ev(nmax, numer)+22);
Ucr(y):=offset+c1[21]+2*sum(realpart(c1[m+20])*cos(2*%pi*(m-1)*y/dp) - imagpart(c1[m+20])*sin(2*%pi*(m-1)*y/dp),m,2,nmax+1);
ev(Ucr(1), numer);
或者只是在最后结合数值
kill(nmax, dp, offset, c1);
Ucr(y):=offset+c1[21]+2*sum(realpart(c1[m+20])*cos(2*%pi*(m-1)*y/dp) - imagpart(c1[m+20])*sin(2*%pi*(m-1)*y/dp),m,2,nmax+1);
env: [dp=1, offset=1, nmax=1, c1=makelist(1/m+m*%i, m, 1, 23)];
ev(Ucr(1), env);