%Euler
clear all
clc
syms x y
f=cos(x)^2-tan(x)*y;
f0=4;
y0=-1;
h=0.1;

%Iteraciones
%Para x=4.1
f1=subs(f,0);
y1=y0+h*f1;
y1=vpa(y1,4)

%Para x=4.1
f2=subs(f,0.1);
y2=y1+h*f2;
y2=vpa(y2,4)

%Para x=4.2
f3=subs(f,0.2);
y3=y2+h*f3;
y3=vpa(y3,4)

%Para x=4.3
f4=subs(f,0.3);
y4=y3+h*f4;
y4=vpa(y4,4)

%Para x=4.4
f5=subs(f,0.4);
y5=y4+h*f5;
y5=vpa(y5,4)

%Para x=4.5
f6=subs(f,0.5);
y6=y5+h*f6;
y6=vpa(y6,4)

%Para x=4.6
f7=subs(f,4.6);
y7=y6+h*f7;
y7=vpa(y7,4)