clc; clear;
format long;
x=0.01:0.005:0.06;
E=10^(-6);
f=@(x)sqrt(2*x+0.4).*atan(cos(3*x+1))
u=@(x)sqrt(2*x+0.4);
v=@(x)atan(x);
w=@(x)cos(3*x+1);
u(x)
v(x)
w(x)
E1=E/3/0.3640
E2=E/3/0.52
E3=E/3

%косинус
y=zeros(length(x),1);
for l=1:length(x);
    pmax=2;
    df=E1*2;% для входа в цикл
     while df > E1;
         p=pmax;  %p-1 так как нумерация - с нуля
         df=1/factorial(abs(2*(p-1)))*(abs((3*x(l)+1)))^(abs(2*(p-1)));
         if df > E1;
             pmax=pmax+1;
         end
     end
    arg=0;
    pmax
	for j=1:pmax; %вычисляем косинус
             arg=arg+(-1)^(abs(j-1))/factorial(2*(j-1))*(abs((3*x(l)+1)))^(abs(2*(j-1)));
    end
    y(l)=arg;
end
y %значения косинуса
(w(x))'
y-(w(x))'

%арктангенс
a=zeros(length(y),1);
for l=1:length(y);
     df=E2*2;% для входа в цикл
     nmax=2;
     while df > E2;
         n=nmax;
         df=1/abs(2*n-1)*(abs((y(l))))^(abs(2*n-1));
         if df > E2;
             nmax=nmax+1;
         end
     end
     arg=0;
     nmax
     for j=1:nmax; %вычисляем арктангенс
             arg=arg+(-1)^(abs(j-1))/abs(2*j-1)*(abs(y(l)))^(abs(2*j-1));
     end
     a(l)=arg;
end
a  %значения арктангенса
(v(y))
a-(v(y))

%корень
z=zeros(length(x),1);%для значений корня
for l=1:length(x);
    f1=0; f2=0;
    df=2*E3;
    kmax=2;
     while df>E3;
         k=kmax;
         f1=100;
         for j=1:(k-1);
             f1=0.5*(f1+(2*x(l)+0.4)/f1);
         end
         f2=0.5*(f1+(2*x(l)+0.4)/f1);
         df=abs(f2-f1);
         if df > E3;
             kmax=kmax+1;
         end
     end
     kmax
     z(l)=f2;
end
z
(u(x))'
z-(u(x))'
znach=zeros(length(x),1);
for l=1:length(x);
     znach(l)=a(l)*z(l);
end;
znach
(f(x))'
znach-(f(x))'