3d6c0455-6df8-4b91-bc37-b9179b8e11aa 7 overlord 3 5 false false 0 0 decimal
&[DATE] &[TIME] - &[FILENAME]
&[PAGENUM] / &[COUNT]
 Restart complete. n is assumed to be positive. ε_%_co is assumed to be positive. 0.85 n * ε.co * n 1 + / 3 1 sys

correct result

 n is assumed to be positive. ε_%_co is assumed to be positive. 17 ε.co n 1 + ^ n - 1 - ( ε.co * + ( * 20 n 1 + ( * / - 3 1 sys

incorrect result

 int() handled by Maxima n is assumed to be positive. ε_%_co is assumed to be positive. 17 ε.co n 1 + ^ n - 1 - ( ε.co * + ( * 20 n 1 + ( * / - 3 1 sys ε.co 0 > n 0 > 2 1 sys 17 n * ε.co * 20 n 1 + ( * / pid=20208 Maxima 5.44.0 http://maxima.sourceforge.net using Lisp SBCL 2.0.0 Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) display2d:false$ (%i2) linel:700$ (%i3) simpsum:true$ (%i4) logsimp:true$ (%i5) keepfloat:true$ (%i6) ratprint:false$ (%i7) load(vect)$ (%i8) load(abs_integrate)$ WARNING: redefining MAXIMA::OPAPPLY in DEFMACRO WARNING: redefining MAXIMA::OPCONS in DEFMACRO (%i9) setup_autoload(mnewton, mnewton)$ (%i10) setup_autoload(draw, draw2d, draw3d)$ (%i11) load(lsquares)$ (%i12) logb(x,y):=log(x)/log(y)$ (%i13) ntewurzel(x,n):=x^(1/n)$ (%i14) el(x,y):=x[y]$ (%i15) el(x,[y]):=if equal(length(y),1) then x[y[1]] else x[y[1],y[2]]$ (%i16) Fit(Data,vars,fun,par,init):=lsquares_estimates_approximate(lsquares_mse(Data,vars,fun),par,initial=init,iprint=[-1,0])$ (%i17) Fit6(Data,vars,fun,par,init,t):=lsquares_estimates_approximate(lsquares_mse(Data,vars,fun),par,initial=init,tol=float(t),iprint=[-1,0])[1]$ (%i18) Residuals(Data,vars,fun,[par]):=float(transpose(matrix(lsquares_residuals(Data,vars,fun,par))))$ (%i19) MSE(Data,vars,fun,[par]):=float(lsquares_residual_mse(Data,vars,fun,par))$ (%i20) round2(x,n):=round(x*10^n)/10^n$ (%i21) check; (%o21) check (%i22) integrate(0.85*(1-(1-x/ε_%_co)^n),x,0,ε_%_co);Is n positive, negative or zero? positive;Is ε_%_co positive, negative or zero? positive;(%o22) (0.85*n*ε_%_co)/(n+1) (%i23) integrate((17*(1-((ε_%_co-x)/ε_%_co)^n))/20,x,0,ε_%_co);Is n positive, negative or zero? positive;Is ε_%_co positive, negative or zero? positive;(%o23) -(17*(ε_%_co^(n+1)+((-n)-1)*ε_%_co))/(20*(n+1)) (%i24) integrate((17*(1-((ε_%_co-x)/ε_%_co)^n))/20,x,0,ε_%_co);Is n positive, negative or zero? positive;Is ε_%_co positive, negative or zero? positive;(%o24) -(17*(ε_%_co^(n+1)+((-n)-1)*ε_%_co))/(20*(n+1)) (%i25) assume(ε_%_co>0,n>0); (%o25) [ε_%_co > 0,n > 0] (%i26) integrate((17*(1-((ε_%_co-x)/ε_%_co)^n))/20,x,0,ε_%_co); (%o26) (17*n*ε_%_co)/(20*(n+1)) (%i27)