Код программы в среде MATLAB

Блок-схема правой части

Блок-схема производной правой части

Блок-схема для решения уравнения методом Ньютона

Блок-схема для главной функции

Код программы в среде Visual C++

#include

#include

#include

using namespace std;

double f1(double x)

{ double y;

y=sin(x)+x-0.5;

return y;

}

double f2(double x)

{ double y;

y=cos(x)+1;

return y;

}

double niu(double eps, double x, int& t)

{ double y;

t=1;

y=x-f1(x)/f2(x);

while(fabs(y-x)>=eps)

{t++;

x=y;

y=x-f1(x)/f2(x);

}

cout<<"t="<

return y;

}

void main()

{ double eps, x, y;

int t;

cout<<"eps=";

cin>>eps;

cout<

cout<<"x Код программы в среде MATLAB=";

cin>>x;

cout<

y=niu(eps,x,t);

cout<<"y="<

}

Результат работы программы

Код программы в среде MATLAB

1. Файл f_niu.m

function z=f_niu(x)

z=sin(x)+x-0.5;

2. Файл f_niu1.m

function z=f_niu1(x)

z=cos(x)+1;

3. Файл f_niu2.m

function z=f_niu2(x)

z=-sin(x);

4. Файл niuton.m

function niuton(f,f1,f2,a,b,eps,x)

y=x-feval(f,x)/feval(f1,x)

t=1;

while abs(y-x)>=eps

x=y;

y=x-feval(f,x)/feval(f1,x);

t=t+1;

end;

x=y

t

Результат работы программы

>> niuton('f_niu','f_niu1','f_niu2',0.2,0.3,0.0000000001,0.25)

y Код программы в среде MATLAB =

0.2513

x =

0.2513

t =


documentaqcggub.html
documentaqcgoej.html
documentaqcgvor.html
documentaqchcyz.html
documentaqchkjh.html
Документ Код программы в среде MATLAB