Matlab - Metodo de Euler (como gráfico la respuesta?)

Problema: Dada la ecuacion diferencial y':sqrt(x^2+y^2) use el metodo de euler para aproximar y(2.3) tomando como numero de pasos n=3 para el proceso iterativo si la condicion inicial es y(2)=0.5
Respuesta y(2.3)=1.166470

Ya esta hecho el problema, lo que me falta es hacer que se grafique.

Código

function f
fprintf('\n \tRESOLUCION DE ECUACIONES DIFERENCIALES POR MEDIO METODO DE EULER\n')
f=input('\nIngrese la ecuacion diferencial de la forma: dy/dx=f(x,y)\n','s');
x0=input('\nIngrese el primer punto x0:\n');
x1=input('\nIngrese el segundo punto x1:\n');
y0=input('\nIngrese la condicion inicial y(x0):\n');
n=input('\nIngrese el numero de pasos n:\n');
h=(x1-x0)/n;
xs=x0:h:x1;
y1=y0;
fprintf('\n''it x0 x1 y1');
for i=1:n
it=i-1;
x0=xs(i);
x=x0;
x1=xs(i+1);
y=y0;
y1=y0+h*eval(f);
fprintf('\n%2.0f%10.6f%10.6f%10.6f\n',it,x0,x1,y1);
y0=y1;
end
fprintf('\n El punto aproximado y(x1) es = %10.6f\n',y1);

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

%function euler

 fprintf('\n \tRESOLUCION DE ECUACIONES DIFERENCIALES POR MEDIO METODO DE EULER\n')

 f=input('Ingrese la ecuacion diferencial de la forma: dy/dx=f(x,y): ','s');

 x0=input('Ingrese el primer punto x0:');

 x1=input('Ingrese el segundo punto x1:');

 y0=input('Ingrese la condicion inicial y(x0):');

 n=input('Ingrese el numero de pasos n:');

 h=(x1-x0)/n;

 fprintf('\n i     xi      yi');

  x=x0;

  X(1)=x;

  y=y0;

  Y(1)=y;

     fprintf('\n%2d  %10.4f  %10.6f%10.6f\n',0,x,y);

 for i=1:n

     y=y+h*eval(f);

     x=x0+i*h;

     X(i+1)=x;

     Y(i+1)=y;

     fprintf('\n%2d  %10.4f  %10.6f%10.6f',i,x,y);

 end

 fprintf('\nEl punto aproximado y(x1) es = %10.6f\n',y1);

 plot(X,Y)

Saludos.
JOSE JEREMÍAS CABALLERO
Asesoría online y Presencial en Matlab
programador en matlab
Servicios de programación matlab
[email protected]
Estimado Usuario de Matlab, el correo es para servicios de programación, toda ayuda gratuita es vía foro.

http://matlabcaballero.blogspot.com

http://www.lawebdelprogramador.com/foros/Matlab/1371532-FORMA_DE_APRENDER_MATLAB.html

fprintf('\nEl punto aproximado y(x1) es = %10.6f\n',(y1));
para la linea 22