Metodo Isoclinas
Metodo Isoclinas
Metodo Isoclinas
En el punto (3, 1), @y @x = 3 En el punto (1, 1), @y @x = 1 Sea m = cte y x y = m entonces x = my =) 8< : y = 1 m x si m 6= 0 x = 0 si m = 0 Si m = 1 =) y = x entonces son las rectas de pendiente m = 1, ahora si m = 2 =) y = 1 2x entonces son las rectas de pendiente m = 2.
de pendiente k). Dibujando las distintas isoclinas se obtiene una representaci_on similar a la del campo de direcciones. Puede tener inter_es identi_car la isoclina para la pendiente 0 pues las soluciones tendr_an generalmente un m_aximo o un m__nimo al pasar por esta isoclina. Ejemplo 1.10 1 Representar las isoclinas de la ED y0 = x + y2. 2 Representar el gr_a_co de contorno (curvas de nivel o isoclinas) de la super_cie z = x + y2.
Clase Pr_actica 1 9
Desarrollo de la soluci_on del ejemplo 1.10 >>[x,y]=meshgrid(0:0.05:3,-2:0.05:2); >> z=x+y.^2; >> isoclinas=contour(x,y,z,20) El resultado en pantalla es el que ofrece la figura 1.6.
Ejemplo 1.11 Representar las isoclinas de la ED y0 = x2 + y2. Soluci_on, paso a paso, del ejemplo 1.11 (ver _g. 1.7) >>[x,y]=meshgrid(-4:0.05:4); >>z=x.^2+y.^2; >>isoclinas=contour(x,y,z,20); Ejemplo 1.12 (a) Representar las isoclinas de la ED y0 = 2x y. (b) >Qu_e tipo de curvas son dichas isoclinas? (c) Representar las isoclinas correspondientes a k = 0, y k = 2. (d) >Qu_e particularidad tiene la correspondiente a k = 2?
Ejemplo 1.13 Construir el campo de direcciones y las curvas de nivel de la ED y0 = sen(x) + y: Desarrollo, paso a paso, de la soluci_on del ejemplo 1.13 (_gura 1.10) >>f=inline('sin(x)+y','x','y'); >>paso=0.5;iz=-3;der=3; >>[x,y]=meshgrid(iz:paso:der,iz:paso:der); >>[n,m]=size(x);dx=ones(n,m); >>z=f(x,y);dy=z; >>hold on,contour(x,y,z,20); >>quiver(x,y,dx,dy);