thanks! my eyes are so trained to mathematica that i didn't even notice i was swapping it there
On Apr 11, 12:52 pm, William Stein <[email protected]> wrote: > 2009/4/11 Chris Chiasson <[email protected]>: > > > > > > > > > does anyone know what could be causing this? (pasted notebook below - > > i would publish this on sagenb.org, but its publishing feature isn't > > working for me right now (keeps giving an invalid url)) > > > ---------------begin notebook after this line----------------- > > Lift Coefficient > > system:sage > > > {{{id=21| > > var('Cv Vr k g SWL Design') > > assume(Cv>1,Vr>=0,k>0,g>0,SWL>0,Design>0) > > eqns={1:Cv==1+Vr*sqrt(k/(g*SWL)),2:SWL*Cv==Design} > > eqns > > /// > > > {1: Cv == sqrt(k)*Vr/(sqrt(g)*sqrt(SWL)) + 1, 2: Cv*SWL == Design} > > }}} > > > {{{id=34| > > power_eq = lambda eq, n: eq.operator()(eq.lhs()^n, eq.rhs()^n) #from > > mhansen on sage-devel > > subs_eq = lambda eq, *args: eq.operator()(eq.lhs().subs_expr > > (*args),eq.rhs().subs_expr(*args)) > > /// > > }}} > > > {{{id=38| > > var('alpha') > > assume(alpha>0) > > eqns[3]=alpha==k*Vr^2/g/Design > > /// > > }}} > > > {{{id=25| > > eqns[4]=subs_eq(power_eq(eqns[1]-1,2),*solve(map(eqns.get,range > > (2,3+1)),k,SWL)[0]);eqns[4] > > /// > > > (Cv - 1)^2 == alpha*Cv > > }}} > > > {{{id=37| > > eqns[5]=(eqns[4]-alpha*Cv).expand();eqns[5] > > /// > > > Cv^2 - alpha*Cv - 2*Cv + 1 == 0 > > }}} > > > {{{id=40| > > solve[eqns[5],Cv] > > You are using solve[...] instead of solve(...). Note the square > versus round brackets. > > William- Hide quoted text - > > - Show quoted text - --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
