MathAction SandBoxDifferentialEquations

	1
Editor: Bill Page Time: 2009/10/16 19:44:19 GMT-7
Note: merge

changed:
-
Differential Equations

\begin{axiom}
y := operator y
deq := D(y(x),x)=(1+y(x)^2)/(1+2*x*y(x))
solve(deq,y,x)
deq := D(y(x),x,3)+D(y(x),x,2)+D(y(x),x,1)+y(x)=x*sin(x)
solve(deq,y,x)
\end{axiom}

test solution

\begin{axiom}
deq := D(y(x),x)=(a+y(x)^2)/(b+2*x^2*y(2*x))
solve(deq,y,x)
deq := D(y(x),x,3)+D(y(x),x,2)+D(y(x),x,1)+y(x)=x*sin(x)
solve(deq,y,x)
\end{axiom}

Linear non homog equation
\begin{axiom}
deq := D(y(x),x) = y(x)^2 + 2 * y(x) + 1
solve(deq,y,x)
z := operator z
deq2 := D(z(x),x) = exp(y(x))/(1 - d * y(x))
solve(deq2,z,x)
\end{axiom}

Differential Equations

fricas

(1) -> y := operator y

$\label{eq1}y$

(1)

Type: BasicOperator?

fricas

deq := D(y(x),x)=(1+y(x)^2)/(1+2*x*y(x))

$\label{eq2}{{y^{\prime}}\left({x}\right)}={\frac{{{y \left({x}\right)}^{2}}+ 1}{{2 \ x \ {y \left({x}\right)}}+ 1}}$

(2)

Type: Equation(Expression(Integer))

fricas

solve(deq,y,x)

$\label{eq3}\frac{{{\left({{y \left({x}\right)}^{2}}+ 1 \right)}\ {\arctan \left({y \left({x}\right)}\right)}}+{y \left({x}\right)}-{2 \ x}}{{2 \ {{y \left({x}\right)}^{2}}}+ 2}$

(3)

Type: Union(Expression(Integer),...)

fricas

deq := D(y(x),x,3)+D(y(x),x,2)+D(y(x),x,1)+y(x)=x*sin(x)

$\label{eq4}\begin{array}{@{}l} \displaystyle {{{y^{\prime \prime \prime}}\left({x}\right)}+{{y^{\prime \prime}}\left({x}\right)}+{{y^{\prime}}\left({x}\right)}+{y \left({x}\right)}}= \ \ \displaystyle {x \ {\sin \left({x}\right)}}$

(4)

Type: Equation(Expression(Integer))

fricas

solve(deq,y,x)

$\label{eq5}\begin{array}{@{}l} \displaystyle \left[{ \begin{array}{@{}l} \displaystyle particular = \ \ \displaystyle {\frac{{{\left({2 \ x}+ 1 \right)}\ {\cos \left({x}\right)}\ {{\sin \left({x}\right)}^{2}}}+{{\left(-{{x}^{2}}+{3 \ x}\right)}\ {\sin \left({x}\right)}}+{{\left({2 \ x}+ 1 \right)}\ {{\cos \left({x}\right)}^{3}}}+{{\left(-{{x}^{2}}-{3 \ x}+ 2 \right)}\ {\cos \left({x}\right)}}}{8}}$

(5)

Type: Union(Record(particular: Expression(Integer),basis: List(Expression(Integer))),...)

test solution

fricas

deq := D(y(x),x)=(a+y(x)^2)/(b+2*x^2*y(2*x))

$\label{eq6}{{y^{\prime}}\left({x}\right)}={\frac{{{y \left({x}\right)}^{2}}+ a}{{2 \ {{x}^{2}}\ {y \left({2 \ x}\right)}}+ b}}$

(6)

Type: Equation(Expression(Integer))

fricas

solve(deq,y,x)

   >> Error detected within library code:
   getfreelincoeff: not a linear ordinary differential equation

Linear non homog equation

fricas

deq := D(y(x),x) = y(x)^2 + 2 * y(x) + 1

$\label{eq7}{{y^{\prime}}\left({x}\right)}={{{y \left({x}\right)}^{2}}+{2 \ {y \left({x}\right)}}+ 1}$

(7)

Type: Equation(Expression(Integer))

fricas

solve(deq,y,x)

$\label{eq8}\frac{-{x \ {y \left({x}\right)}}- x - 1}{{y \left({x}\right)}+ 1}$

(8)

Type: Union(Expression(Integer),...)

fricas

z := operator z

$\label{eq9}z$

(9)

Type: BasicOperator?

fricas

deq2 := D(z(x),x) = exp(y(x))/(1 - d * y(x))

$\label{eq10}{{z^{\prime}}\left({x}\right)}= -{\frac{{e}^{y \left({x}\right)}}{{d \ {y \left({x}\right)}}- 1}}$

(10)

Type: Equation(Expression(Integer))

fricas

solve(deq2,z,x)

$\label{eq11}\begin{array}{@{}l} \displaystyle \left[{particular ={\int^{ \displaystyle x}-{{\frac{{e}^{y \left({\%D}\right)}}{{d \ {y \left({\%D}\right)}}- 1}}\ {d \%D}}}}, \: \right. \ \ \displaystyle \left.{basis ={\left[ 1 \right]}}\right]$

(11)

Type: Union(Record(particular: Expression(Integer),basis: List(Expression(Integer))),...)

Edit detail for SandBoxDifferentialEquations revision 1 of 1