{VERSION 1 0 "X11/Motif" "1.0"}{GLOBALS 3 1}{FONT 0 "-adobe-helve
tica-bold-r-normal--*-140-*" "helvetica" "Helvetica-Bold" 8 14 0 
"Helvetica-Bold" 12}{FONT 1 "-adobe-times-medium-r-normal--*-180-
*-*-*-*-*-*" "times" "Times-Roman" 4 18 64 "Times-Roman" 18}{FONT
 2 "-adobe-courier-medium-r-normal--*-140-*" "courier" "Courier" 
4 14 192 "Courier" 12}{FONT 3 "-adobe-helvetica-medium-r-normal--
*-180-*-*-*-*-*-*" "helvetica" "Helvetica" 4 18 0 "Helvetica" 18}
{FONT 4 "-adobe-helvetica-medium-r-normal--*-180-*-*-*-*-*-*" "he
lvetica" "Helvetica" 4 18 0 "Helvetica" 18}{SCP_R 1 0 30{INP_R 2 
0 "> "{TEXT 0 29 "with(DEtools): with(plots):  "}}{SEP_R 3 0}
{COM_R 4 0{TEXT 1 0 ""}}{INP_R 5 0 "> "{TEXT 0 94 "p1:=DEplot(dif
f(y(x),x)=2*y(x)-1,[x,y],0..0.5,\{[0,1]\},stepsize=.1,method=`Eul
er`,thickness=1);"}}{SEP_R 6 0}{COM_R 7 0{TEXT 1 0 ""}}{INP_R 8 0
 "> "{TEXT 0 95 "p2:=DEplot(diff(y(x),x)=2*y(x)-1,[x,y],0..0.5,\{
[0,1]\},stepsize=.05,method=`Euler`,thickness=2):"}}{SEP_R 9 0}
{INP_R 10 0 "> "{TEXT 0 95 "p3:=DEplot(diff(y(x),x)=2*y(x)-1,[x,y
],0..0.5,\{[0,1]\},stepsize=.01,method=`Euler`,thickness=3):"}}
{COM_R 11 0{TEXT 1 0 ""}}{INP_R 12 0 "> "{TEXT 0 46 "diffeq1:=dif
f(y(x),x)=2*y(x)-1; inits:=y(0)=1;"}}{OUT_R 13 0 12{DAG :3n4\`dif
feq1`=3(3n4\`diff`,3(3n3\`y`,2n3\`x`pE+5pAj2x0002i2x0001j2x0001}}
{OUT_R 14 0 12{DAG :3n4\`inits`=3(3n3\`y`,2j2x0000j2x0001}}{INP_R
 15 0 "> "{TEXT 0 35 "sol1:=dsolve(\{diffeq1,inits\},y(x));"}}
{OUT_R 16 0 15{DAG :3n4\`sol1`=3(3n3\`y`,2n3\`x`+5/3j2x0001j2x000
2pD(3n3\`exp`,2+3p9pFpC}}{INP_R 17 0 "> "{TEXT 0 29 "p4:=plot(rhs
(sol1),x=0..0.4):"}}{INP_R 18 0 "> "{TEXT 0 0 ""}}{INP_R 19 0 "> \+
"{TEXT 0 23 "display([p1,p2,p3,p4]);"}}{INP_R 20 0 "> "{TEXT 0 96
 "p5:=DEplot(diff(y(x),x)=y(x)^2+x^2,[x,y],0..0.5,\{[0,1]\},steps
ize=.1,method=`Euler`,thickness=1):"}}{INP_R 21 0 "> "{TEXT 0 97 
"p6:=DEplot(diff(y(x),x)=y(x)^2+x^2,[x,y],0..0.5,\{[0,1]\},stepsi
ze=.05,method=`Euler`,thickness=2):"}}{INP_R 22 0 "> "{TEXT 0 97 
"p7:=DEplot(diff(y(x),x)=y(x)^2+x^2,[x,y],0..0.5,\{[0,1]\},stepsi
ze=.01,method=`Euler`,thickness=3):"}}{INP_R 23 0 "> "{TEXT 0 48 
"diffeq2:=diff(y(x),x)=y(x)^2+x^2; inits:=y(0)=1;"}}{OUT_R 24 0 
23{DAG :3n4\`diffeq2`=3(3n4\`diff`,3(3n3\`y`,2n3\`x`pE+5*3pAj2x00
02j2x0001*3pEp14p16}}{OUT_R 25 0 23{DAG :3n4\`inits`=3(3n3\`y`,2j
2x0000j2x0001}}{INP_R 26 0 "> "{TEXT 0 35 "sol2:=dsolve(\{diffeq2
,inits\},y(x));"}}{OUT_R 27 0 26{DAG :3n4\`sol2`=3(3n3\`y`,2n3\`x
`+5*Bp9j2x0001+5*7(3n4\`GAMMA`,2/3j2x0003j2x0004j2x0002+5*3p12p1C
pEn3\`Pi`i2x0001p25(3n4\`BesselY`,3/3pEp1A+3*3p9p1C/3pEp1CpEp25(3
n4\`BesselJ`p2CpEp25p12p1Cp1Ep25(3p29,3/3i2x0003p1Ap30pEpE*7p9pEp
10p25(3p3Ap46pEp25}}{INP_R 28 0 "> "{TEXT 0 29 "p8:=plot(rhs(sol2
),x=0..0.5):"}}{INP_R 29 0 "> "{TEXT 0 23 "display([p5,p6,p7,p8])
;"}}{SEP_R 30 0}{INP_R 31 0 "> "{TEXT 0 0 ""}}}{END}