Thanks for helping. I used a mathematical approach for 4th order derivative i.e substituing double derivative as a variable, solving for 2nd order differential equation and substituting back and again solved for 2nd order differentiation. here is code: var('w,x,E,L,k1,k2') y = function('y', x) w= function('w' , x) q = function('q', x) assume(L>0) assume(E>0) q=x de=E*L*diff(y,x,2)==q y_res=desolve(de,y,ivar=x,ics=[L,0,0]) des=diff(w,x,x)-y_res==0 dess=desolve(des,w,ivar=x,ics=[0,0,0]) print "Solution of bernoulli's equation:",dess #####Remeber plot can't be formed without giving values of constant############### E=6 L=10 p=plot( 1/120*(20*L^3*x^2 - 10*L^2*x^3 + x^5)/(E*L),(x,0,1),thickness=3) p.show()
-- Priyanka Kapoor priyankacool10.wordpress.com Linux User Group, Ludhiana -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org