In [7]:
from sympy import *
init_printing() #muestra símbolos más agradab
R=lambda n,d: Rational(n,d)
In [8]:
x,y,t=symbols('x,y,t')
In [2]:
Si uno plantea el sistema queda de esta forma
In [5]:
t=symbols('t')
A=Matrix([[-R(1,50),0,0],[R(1,50),-R(2,50),0],[0, +R(2,50),-R(3,50)]])*t
A
x0=Matrix([[10],[10],[10]])
A, x0
Out[5]:
"Resolvamoslo"
In [6]:
n=10
g=Matrix([[0],[0],[0]])
for i in range(n+1):
g=g+A**i*x0
g
Out[6]:
"las soluciones son iguales"
Reslover el problema
\begin{equation*} \left\{ \begin{array}{cc} -u_{x}+u_{y}=u^2 & (x,y)\in\mathbb{R}^{2}\\ u(x,0)=\frac{1}{2}e^{-x} & x\in \mathbb{R} \end{array} \right. \end{equation*}La solucion del problema es $u(x,t)=\frac{1}{2e^{x+y}-y}$. Verifiquemoslo
In [9]:
u=(2*exp(x+y)-y)**(-1)
u
Out[9]:
In [11]:
((-u.diff(x)+u.diff(y))-u**2).factor()
Out[11]:
In [12]:
u.subs(y,0)
Out[12]:
Sea $y=ax$ una recta que pasa por el origen. Sea $(x_0,y_0):=(x_0,y (x_0))$ el punto de corte entre la recta $y=ax$ y una solución de la ecuación $y'(x)=f(x,y)$, y sea $\alpha$ el ángulo que forman las rectas $y=ax$ con la recta tangente $y=y'(x_0)(x-x_0)+y_0$ a la curva $(x,y(x))$ en el punto $(x_0,y_0)$. Luego, el $\cos (\alpha)$ se puede calcular de esta manera
$$ \cos (\alpha)=\frac{u}{|| u||} \cdot \frac{v}{|| v||}$$donde $u$ y $v$ son vectores paralelos a las rectas $y=ax$ e $y=y'(x_0)(x-x_0)+y_0$ respectivamente.
En la igualdad anterior podemos tomar $u=(x-x0,ax-ax_0)$ y $v=(x-x_0,y'(x_0)(x-x_0)+y_0-y_0)$, para algún $x\neq x_0$. Luego,
\begin{equation*} \begin{split} \cos (\alpha)&=\frac{(x-x0,ax-ax_0)}{\sqrt{1+a^2}|x-x_0 |}\cdot \frac{(x-x_0,y'(x_0)(x-x_0))}{\sqrt{1+(y'(x_0))^2}|x-x_0 |}\\ &=\frac{(x-x_0)^2+(a+y'(x_0))(x-x_0)^2}{\sqrt{(1+a^2)(1+(y'(x_0))^2)}(x-x_0)^2}\\ &=\frac{1+(a+y'(x_0))}{\sqrt{(1+a^2)(1+(y'(x_0))^2)}}\\ &=\frac{1+(a+f(x_0,y_0))}{\sqrt{(1+a^2)(1+(f(x_0,y_0))^2)}}\\ &=\frac{1+(a+f(1,\frac{y_0}{x_0}))}{\sqrt{(1+a^2)(1+(f(1,\frac{y_0}{x_0}))^2)}}\\ &=\frac{1+(a+f(1,a))}{\sqrt{(1+a^2)(1+(f(1,a))^2)}}. \end{split} \end{equation*}Y esta última expresión es independiente del punto $(x_0,y_0)$ donde la recta interseca a la curva solución de $y'=f(x,y)$.
In [ ]:
In [ ]:
In [ ]:
In [ ]: