/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 17.10.1 ] */
/* [wxMaxima: title   start ]
Procedury - horni a dolni soucet
   [wxMaxima: title   end   ] */


/* [wxMaxima: input   start ] */
RiemannSum(fn,a,b,n,opt):=
block([xx,s],
xx(i,n):= a+(i+opt/2)*(b-a)/n,
s: sum(ev(fn,x=xx(i,n)),i,0,n-1)*(b-a)/n,
float(s)
)$
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannRectangles(fn,a,b,n,opt):=
block([xi,xx,rects,wxd],
xi(i,n):= a+i*(b-a)/n,
xx(i,n):= a+(i+opt/2)*(b-a)/n,
rects(n):=makelist(rectangle([xi(k,n),0], [xi(k+1,n),ev(fn,x=xx(k,n))]),k,0,n-1),
wxd(n):=apply(wxdraw2d,
append([xrange=[a-(b-a)/4, b+(b-a)/4],color=blue],
rects(n),
[transparent=true, explicit(fn,x,a, b)]
)
),
wxd(n)
)$
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
Priklad procedur na konkretni funkci - postupna aproximace
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
float(integrate(6*x-x^2,x,0,2));
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannRectangles(6*x-x^2,0,2,20,0);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannSum(6*x-x^2,0,2,20,0);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannRectangles(6*x-x^2,0,2,50,0);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannSum(6*x-x^2,0,2,50,0);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannRectangles(6*x-x^2,0,2,20,2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannSum(6*x-x^2,0,2,20,2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannRectangles(6*x-x^2,0,2,50,2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
RiemannSum(6*x-x^2,0,2,50,2);
/* [wxMaxima: input   end   ] */



/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$