• Book	Advanced Engineering Mathematics
  Author	Erwin Kreyszig
  Edition	8th Edition

code_begins

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))

code_ends

• Solution by:

  • Aashita Kesarwani , Student, IIT Roorkee
  • Swati Ahuja, Student, IIT Roorkee
  • Praveen Kumar, Student, IIT Roorkee.