O’zbekiston respublikasi oliy ta’lim vazirligiga qarashli muhammad al-xorazimiy nomidagi toshkent axborot texnologiyalari universiteti irisov Lazizbek laboratoriya ishi-1

O’ZBEKISTON RESPUBLIKASI OLIY TA’LIM VAZIRLIGIGA QARASHLI MUHAMMAD AL-XORAZIMIY NOMIDAGI TOSHKENT AXBOROT TEXNOLOGIYALARI UNIVERSITETI

Irisov Lazizbek

LABORATORIYA ISHI-1

Fan nomi: Algaritmlarni loyihalash

Variant 10.
Berilgan integral qiymatini to‘g‘ri to‘rtburchaklar, trapetsiyalar, Simpson usullarida ε>0 aniqlikda hisoblang. Aniqlikka erishganlik sharti sifatida |S2n-Sn|< ε tengsizlikdan foydalaning (boshlang‘ich n=10 deb olish mumkin).

Berilgan tenglamaning taqribiy yechimini ε>0 aniqlikda urinmalar (Nyuton) va vatarlar usullarida hisoblang. Aniqlikka erishganlik sharti sifatida |xn+1-xn|<ε tengsizlikdan foydalaning. Boshlang‘ich x0ϵ(a;b) olinadi.

Nyuton (Urinmalar) usuli.

#include

#include

using namespace std;

double func(double );

double birinchiHosilaFunc(double );

void newton_raphson(double ,double );

int main()

{
double x2,e;
cout<<"Boshlangich taxminni kiriting: ";
cin>>x2;
cout<<"Istalgan aniqlikni kiriting e: ";
cin>>e;
newton_raphson(x2,e);
return 0;
}
double func(double x)
{
double res= (pow(x,7)-pow(x,2)+4);
return res;
}
double birinchiHosilaFunc(double x)
{
return (7*pow(x,6) - 2*x);
}
void newton_raphson(double x2, double e)
{
double x1;
do
{
x1=x2;
x2 = x1 - func(x1)/ birinchiHosilaFunc (x1);
}while(fabs(x2-x1) >= e);
cout<<"\nTenglam ildizi: "<<(x2);
}
Vatarlar usuli

#include

#include

using namespace std;

float f(float x) {
float f = pow(x,7)-pow(x,2)+4;
return f;
}
void secant(float x1, float x2, float E) {
float n = 0, xm, x0, c;
if (f(x1) * f(x2) < 0) {
do {
x0 = (x1 * f(x2) - x2 * f(x1)) / (f(x2) - f(x1));
c = f(x1) * f(x0);
x1 = x2;
x2 = x0;
n++;
if (c == 0)
break;
xm = (x1 * f(x2) - x2 * f(x1)) / (f(x2) - f(x1));
} while (fabs(xm - x0) >= E);
cout << "tenglama ildizi=" << (x0) << endl;
cout << "qadamlar soni = " << n << endl;
} else
cout << "javobi topilmadi";
}
int main() {
float x1 = 0, x2 = 1, E = 0.0001;
secant(x1, x2, E);
return 0;
}