/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 0.8.3a ] */

/* [wxMaxima: input   start ] */
H:R*(1-exp(-a*N))/N;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(H,R=1,a=1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([ev(H,R=20,a=.1),ev(H,R=10,a=.1)],[N,0,50]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
limit(H,N,0);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
solve(ev(H,R=20,a=.1,N=n1)-ev(H,R=10,a=.1,N=n2),n1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
load("mnewton")$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(H,R=20,a=.1,N=n1)-a*10;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
mnewton(ev(H,R=20,a=.1,N=n1)-.1*10,n1,10);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ev(H,R=20,a=.1,N=15.94);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([ev(H,R=20,a=.1),ev(H,R=10,a=.1), [discrete,[[0,1],[50,1]]]],[N,0,50]);
/* [wxMaxima: input   end   ] */

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