var('t')
y=function('y',t)
deq=diff(y,t)-10-0.025*y==0
plot(desolve(deq,[y,t],ics=(0,40)),(t,0,10))

Kreyszig-1.4-2 (last edited 2010-12-18 06:49:21 by manishadeepthi)