Maybe (I don‘t have my pc on me right now) the “de = L*E*diff(y,x,2)==q“ is
making some problems because that is also a comparing method (and so, de is
True or False)
El 29/05/2012 07:29, "Priyanka Kapoor" <anjalicool.kapoor...@gmail.com>
escribió:
> 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
>
--
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
