var('x')
y=function('y',x)
deq=diff(y)-(-0.046)*(y-32)==0
plot(desolve(deq,[y,x],ics=(0,66)),(x,0,10),ymin =0, ymax=67)+plot(32,(x,0,10))

Kreyszig-1.4-3 (last edited 2010-12-18 06:28:48 by manishadeepthi)