Differences between revisions 1 and 4 (spanning 3 versions)
⇤ ← Revision 1 as of 2010-08-09 14:09:43 →
Size: 452
Editor: SWATI AHUJA
Comment:
← Revision 4 as of 2010-08-12 10:31:37 → ⇥
Size: 739
Editor: Aashita
Comment:
Deletions are marked like this. Additions are marked like this.
Line 8: Line 8:
var('t')
var('a')
var('b')
var('k')
y=function('y',t)
print desolve(diff(y,t)-a + b*cos(2*pi*t/24)+k*y,[y,t])		 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))
Line 19: Line 26:
   * Aashita Kesarwani , Student, IIT Roorkee
Line 20: Line 28:
Kesarwani Aashita Chandrashekhar, Student, IIT Roorkee
Praveen Kumar, Student, IIT Roorkee		    * Praveen Kumar, Student, IIT Roorkee.

  • Book

    Advanced Engineering Mathematics

    Author

    Erwin Kreyszig

    Edition

    8th Edition 

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

  • Solution by:

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

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