Octave - Ayuda Metodo Euler

 
Vista:

Ayuda Metodo Euler

Publicado por Walter Amarilla (1 intervención) el 15/02/2020 21:34:40
Buenas:
Consulto por un ejercicio de metodo de Euler para programar que no me sale.

Me piden graficar simultaneamente en la region [0.10]x[0,10] las soluciones del problema
y'(t)=(y(t)-5)*(cos^2(t)-0.5)
con y(0)=k
para k de 0,1, hasta 10 en pasos h=0.01

Tengo mi programa de euler asi:

1
2
3
4
5
6
7
8
9
function [y]= Euler(to, tf, yo, h, f)
  t = to:h:tf;
  m = length (t);
  y(1)= yo;
  for i=1:(m-1)
    y(i+1)= y(i)+h*f(t(i),y(i));
  endfor
  plot(t,y)
endfunction

en comand windows cargo los valores para mis variables pero no se como modificar el programa para todas las k y como que se grafiquen simultaneamentes.

intente haciendo esto en mi programa pero lo que me grafica no es lo que me piden.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
function [y]= Eulerk(to, tf, yo, h, f)
  t = to:h:tf;
  m = length (t);
  for k = 0:10;
      y = zeros(m,1);
      y(1)= k;
      for i=2:(m-1)
          y(i+1)= y(i)+h*f(t(i),y(i));
      end
      plot(t,y)
      hold on
  end
 
endfunction



muchas graciasss.
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder