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

code_begins

x = var('x')
y = function('y',x)
deq1 = diff(y,x,3)-3*diff(y,x,2)+3*diff(y,x)-y==0
soln = desolve_laplace(deq1 ,  [y,x] , ics=[2,2,10])
soln
plot(y)

code_ends

• Solution by:

  • V Krishna Harsha , Student , ASR College Of Enginnering