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
