t = var('t') # time in hrs  
a = var('a') # Avg. secretion rate
b = var('b')
k = var('k')
y = function('y',t) # Amount of a certain hormone in blood at time = 't'
sol=desolve (diff(y,t) - a + b*cos(2*pi*t/24) + k*y, [y,t])
print "General solution= ", sol
a=1
b=1
k=1
sol1=desolve (diff(y,t) - a + b*cos(2*pi*t/24) + k*y, [y,t], ics=(0,2))
print "Specific solution= ", sol1
plot(sol1,(t,0,20))

SageDays/Kreyszig-1.6-40 (last edited 2010-08-12 10:31:37 by Aashita)