On Oct 23, 6:12 pm, Martin Rubey <[email protected]> wrote: > Hm, strange enough it seems that this is not implemented. Does the > following workaround help? > > r := solve([eq1,eq2], [x,y], t) > solve(eval(matrix(r.basis), t=0), [1,0]) > > (18) [particular= [1,0],basis= []]
I'm not sure. r.basis does not seem to contain any arbitrary constants to be solved for. In this simple case it should just pick out the first soln in r.basis (x(t) = cos(t) and y(t) = sin(t)) which satisfy the initial conditions. > > subst(inverseLaplace((R*g2849^3+(G*R+G)*g2849^2+(R^2+R > > +G^2)*g2849+2*G*R)/(10*(R*g2849^4+(G*R > > +G)*g2849^3+subst(inverseLaplace((R*g2849^3+(G*R+G)*g2849^2+(R^2+R > > +G^2)*g2849+2*G*R)/(10*(R*g2849^4+(G*R+G)*g2849^3+(2*R^2+2*R > > +G^2)*g2849^2+(3*G*R+G)*g2849+R^2+R)),g2849,t), [G=2,R=100]) > > there are some parens missing, could you please clean up? My apologises; I accidentally copied the expression from a FriCAS session; try: subst(inverseLaplace((R*g2849^3+(G*R+G)*g2849^2+(R^2+R +G^2)*g2849+2*G*R)/(10*(R*g2849^4+(G*R+G)*g2849^3+(2*R^2+2*R +G^2)*g2849^2+(3*G*R+G)*g2849+R^2+R)),g2849,t),[G=2,R=200]) That should work. Regards, Freddie. -- You received this message because you are subscribed to the Google Groups "FriCAS - computer algebra system" group. 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/fricas-devel?hl=en.
