Смекни!
smekni.com

Вычисление площади сложной фигуры методом имитационного моделирования

МОСКОВСКИЙГОСУДАРСТВЕННЫЙИНСТИТУТЭЛЕКТРОННОЙТЕХНИКИ

(ТЕХНИЧЕСКИЙУНИВЕРСИТЕТ)


Расчет площади сложной фигурыс помощью метода

имитацеонногомоделирования.


ЛогвиненкоВ.


Москва. 1995 г.


Задание:Разработатьпрограмму,позволяющуюс помощью методаимитационногомоделированиярассчитатьплощадь сложнойфигуры, ограниченнойсверху кривойU=Y1(x) , снизу V=Y2(x).


1. Длярешения даннойзадачи применимследующийметод.


Ограничимзаданную фигурупрямоугольником,стороны которогопроходят:

через точкимаксимальногои минимальногозначения функцийи параллельныосям абсцисс;

через левуюи правую граничныеточки областиопределенияаргументаи параллельныосям ординат.

Используядатчик случайныхчисел разыгрываютсякоординатыслучайнойточки из этогопрямоугольника. Проверяемпопадаете точкив заданнуюфигуру. Знаяплощадь прямоугольникаи отношениепопавших точекк их общемучислу разыгранных,можно оценитьплощадь интересующейнас фигуры.


2.Техническиехарактеристикиобъекта исследования:


2.1. Диапазонзначенийпараметровзадачи.


Множествокривых ограничимполиномамитретьего порядка,в виду тогочто полиномыболее высокогопорядка сильноувеличиваютвремя вычисления.Причем длянаглядностирешения вполнедостаточнопорядка "3".


Коэффициентыполинома ограничимдиапазоном [-100,100] .


Область определенияограничимдиапазоном[-100,100].


Эти ограничениявведены дляболее наглядногорешения задачи,и изменитьих не с техническойточки зренияне сложно.


3.Решение задачи.


Данная задачарешена в средеTurbo C. Для решенияпотребовалосьобщую задачуразбить на несколько небольших задач (процедур).

А именно отдельно(в виде процедур)были решенызадачи


-ввод параметров;|

процедура get_poly|

|

-сообщениеоб ошибке привводе; |ФайлWINDOW.C

процедура talkerror |

|

-рисованиерамки окна; |

процедура border|


-вычислениеминимальногои |

максимального значениифункций ;|

процедура f_max|

|

-вычислениезначенияполинома в|

заданнойточке;|Файл MATIM.C

процедура fun|

|

-вычислениекорней кубичного|

уравнения;|

процедура f_root|


-вычислениеинтегралачисленным|

методом;|

процедура i_num|

|Файл F_INTEGER.C

-вычислениеинтеграла спомощью|

имитационногомоделирования;|

процедура i_rand|


-инициализацияграфическогорежима|

процедура init|

|

-обводка непрерывногоконтура|ФайлDRAFT.C

процедура f_draft|

|

- вырисовкаосей координат|

процедура osi|


-вырисовкиграфиков функцийи |Файл DRAFT_F.C

штриховказаданнойплощади|

процедура draft_f|


-вырисовкаграфиков вычисления|

площади разнымиметодами ивывод |ФайлDRAFT_N.C

таблицырезультатоввычисления|

процедураdraft_n|


Схема алгоритма имеет вид:


µ§


4. Описаниепроцедур используемыйв программе.


4.1 ФайлWINDOW.C.


4.1.1 Процедураввода параметров.

void get_poly( float *b3,float *b2,float *b1,float *b0,//-коэффициентыполинома Y1

fliat *c3,float *c2,float *c1,float *c0, //-коэффициентыполинома Y2

float *x1,float *x2, //область определения[x1,x2]

int *N )// количествообращений кгенератору//случайныхчисел


4.1.2 Процедурарисованиярамки окна.

void border(int sx, int sy, int en, int ey) // рисуетрамку с координатамилевого верхнего // угла (sx,sy) икоординатамиправого нижнего // угла (ex,ey)


4.1.3 Процедурасообщенияоб ошибке привводе.

void talkerror(void) -

Процедура подаетзвуковой сигнали выводит наэкран сообщениеоб ошибке привводе.


4.2.Файл MATIM.C


4.2.1 Процедуравычислениямаксимальногои минимальногозначенийфункций назаданноминтервале.

void f_max(float b3,float b2,float b1,float b0, //-коэффициентыполинома Y1

fliat c3,float c2,float c1,float c0, //-коэффициентыполинома Y2

float x1,float x2, // областьопределения[x1,x2]

float *amin, float *amax) // минимальноеи максимальноезначения//функций

4.2.2 Процедуравычислениязначенияполинома вданной точке.

float fun(float b3,float b2,float b1,float b0, //-коэффициентыполинома

float x)

Возвращаетзначениеполинома вточке х.


4.2.3 Процедуравычислениякорней кубичногоуравнения.

int f_root(float b3,float b2,float b1,float b0, //-коэффициентыполинома Y1

fliat c3,float c2,float c1,float c0, //-коэффициентыполинома Y2

float x1,float x2, // областьопределения[x1,x2]


float e, // точностьвычислениякорней

float *k1,float *k2,float *k3) // значения корней // функций


Возвращаетколичество действительныхкорней на данноминтервале.


4.3. ФайлF_INTEGER.C


4.3.1 Процедуравычисленияплощади сложнойфигуры численнымметодом.

float f_num(float b3,float b2,float b1,float b0, //-коэфициентыполинома Y1

fliat c3,float c2,float c1,float c0, //-коэфициентыполинома Y2

float x1,float x2) // областьопределения[x1,x2]

Вычисляетплощадь сложнойфигуры.


4.3.2 Процедуравычисленияплощади сложнойфигуры c помощьюметода имитационногомоделрования

float f_(float b3,float b2,float b1,float b0, //-коэфициентыполинома Y1

fliat c3,float c2,float c1,float c0, //-коэфициентыполинома Y2

float x1,float x2, // областьопределения[x1,x2]

float fmin,float fmax, // минимальноеи максимальноезначения //функций наданном интервале

int n)// количествообращений кгенератору// случайныйчисел


Вычисляетплощадь сложнойфигуры с помощьюметода имитационногомоделирования.


4.4 ФайлDRAFT.C


4.4.1 Процедураинициализацииграфическогорежима.

void init (void)


4.4.2 Процедураобводки непрерывногоконтура.

void f_draft (float b0,float b1,float b2,float b3, //-коэфициентыполинома

float x1,float x2) // областьопределения[x1,x2]


4.4.3 Процедуравырисовки осейкоординат.

void osi ( float x1, float x2, // областьопределенияфункций

float b) // маштабныйкоэфициентрасчитываетсяпо формуле

// b= j - Fmin*(i-j) /(Fmax - Fmin)

// гдеi,j - задают положениеграфика наэкране

// Fmin,Fmax - минимальноеи максимальноезначения //функций наданном интервале


4.5Файл DRAFT_F.


4.5.1 Процедуравырисовкиграфиков функций.

void draft_f (float b3,float b2,float b1,float b0, //-коэфициентыполинома Y1

fliat c3,float c2,float c1,float c0, //-коэфициентыполинома Y2

float x1,float x2, // областьопределения[x1,x2]

float fmin,float fmax, // минимальноеи максимальноезначения //функций наданном интервале

int k, int i, int l, int j)// координаты,задающиеположение //графикана экране

4.6 ФайлDRAFT_N.


4.6.1 Процедуравырисовкиграфиков значенийполщадей расчитанныхчислвым методоми методомимитационногомоделированияв зависимостиот количестваобращений кгенераторуслучайныхчисел.

void draft_e (float b3,float b2,float b1,float b0, //-коэфициентыполинома Y1

fliat c3,float c2,float c1,float c0, //-коэфициентыполинома Y2

float x1,float x2, // областьопределения[x1,x2]

float fmin,float fmax, // минимальноеи максимальноезначения //функций наданном интервале

float Sn, //площадь рассчитаннаячисловым методом

int k, int i, int l, int j)// координаты,задающиеположение //графикана экране


4.7Файл SQ.C

Все файлыобъединеныв главнойпрограмме SQ.C,которая являетсяосновной икоординируетработу процедур.


5Использованиепрограммы.

Для использованияданной программынеобходимаоперационнаясреда MS DOS,

файл egavega.bgi, и собственносама скомпилированнаяпрограммаsq.exe.


6 Исходныйтекст программыдан в приложении№1.


7 Тесовый примерпоказан вприложении№2.


8Список использованнойлитературы.


8.1 Язык программированияСи для персональногокомпьютера.

С.О.Бочков, Д.М.Субботин.

8.2 С++ . Описаниеязыка программирования.

БьярнСтраустрап.

8.3 TURBO C. User's Guide. BorlandInternational, Inc. 1988.

8.4 TURBO C. Reference Guide. BorlandInternational, Inc. 1988.


9Заключение.

9.1 Сопоставлениерезультатовработы с тербованиемзадания.

Сопоставляярезультатыработы с требованиемзадания, можносказать чтозадача решенав полной мере,за исключением,быть можетобщностиотносительновозможностирасчета длямногие классовфункций. Норешение болееобщей задачи( т.е. возможностьрасчета длямногих классовфункций ) представляетсязначительноболее громоздким,и вообще являетсяотдельнойзадачей. Поэтомуавтор не счелнужным разрабатыватьалгоритм вводамногих функцийи заострилвниманиесобственнона самой задаче- расчете площадисложной фигурыс помощью методаимитационногомоделированияи сравнениеэтого методас числовымиметодами.


9.2 Рекомендациипо улучшениюпрограммы.

При разработкепрограммы авторупустил возможностьработы с числовымимассивами.Поэтому, можноулучшитьпрограммупереписав рядпроцедур подмассивы , чтосделает программуменее массивнойи более наглядной.Широкое

возможностипо улучшению программы вобласти разработкиалгоритмовввода различныйклассов функций.


Для решениязадачи методомимитационногомоделированияограничимданную



§



§

Приложение1. Текст программы.


Файлsq.c


/*

ПpогpаммаSQ основная

*/

#include

#include

#include

#include

#include

#include "matim.c"

#include "window.c"

#include "f_integr.c"

#include "draft.c"

#include "draft_f.c"

#include "draft_e.c"

intk=20,i=15,l=270,j=140;

void main(void)

{

floatb0,b1,b2,b3,c0,c1,c2,c3,x1,x2,maxb,maxc,minb,minc,min,max,S;

int N;

do{

closegraph();

get_poly(&b3,&b2,&b1,&b0,&c3,&c2,&c1,&c0,&x1,&x2,&N);

f_max(b3,b2,b1,b0,x1,x2,&minb,&maxb);

f_max(c3,c2,c1,c0,x1,x2,&minc,&maxc);

max=(maxb>maxc)?maxb:maxc;

min=(minb

S=i_num(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2);

init();

draft_f(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,k,i,l,j);

draft_e(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,S,k,i+180,l+100,j+160,N);

setcolor(2);

outtextxy(0,340," Press q for exit ");

} while (( getch())!= 'q');


}


Файлmatim.c


/* Подпpогpаммасодеpжит пpоцедуpыматематическойобpаботкифункций*/

#include

#include

#include

#include

#include

/* Вычислениемаксимальногои минимального

значенияфункции назаданом интеpвале*/

void f_max(floata3,float a2,float a1,float a0,float x1,float x2,float *amin,float*amax)

{

floatdx,x,Fx,Fx1,Fmax,Fmin;

dx=(x2-x1)/500;

x=x1;

Fx1=a3*x*x*x+a2*x*x+a1*x+a0;

Fmax=Fx1;

Fmin=Fx1;

do {

x=x+dx;

Fx=a3*x*x*x+a2*x*x+a1*x+a0;

if (Fx>=Fmax)

Fmax=Fx;

if (Fx

Fmin=Fx;

} while ( x

*amin=Fmin;

*amax=Fmax;

}

/*Вычислениекоpней кубичногоуpавнения */

int f_root(floata0,float a1,float a2,float a3,float x1,float x2,float e,float*k1,float *k2,float *k3)

{ floatku1,ku2,ku3,x,a,b;

int c=0;

x=x1;

do

{

a=a3*pow(x,3)+a2*pow(x,2)+a1*x+a0;

x+=e;

b=a3*pow(x,3)+a2*pow(x,2)+a1*x+a0;

if (a*b

{ c++;

switch(c) {

case 1: ku1=x;

break;

case 2: ku2=x;

break;

case 3: ku3=x;

break;

default: printf("\nВнимание !!! \nОшибка в matim.c(f_root).");

break;

};

}

} while (x

*k1=ku1;

*k2=ku2;

*k3=ku3;

return c;

}

float fun(floata3,float a2,float a1,float a0,float x)

{

float s;

s=a3*x*x*x+a2*x*x+a1*x+a0;

return (s);

}


Файлwindow.c


/* Подпpогpаммыpаботы с окнами*/

#include

#include

#include

#include

#include


/*функцияpисования pамкиокна */


void border(int sx,intsy,int ex,int ey){

int i;

for (i=sx+1;i

gotoxy(i,sy);

putch(205);

gotoxy(i,ey);

putch(205);

}

for (i=sy+1;i

gotoxy(sx,i);

putch(186);

gotoxy(ex,i);

putch(186);

}

gotoxy(sx,sy);putch(201);

gotoxy(sx,ey);putch(200);

gotoxy(ex,sy);putch(187);

gotoxy(ex,ey);putch(188);

}

void talkerror(void)

{

textcolor(15);

textbackground(4);

gotoxy(1,18);

cprintf("ATTATETION ! DATE ERROR . Press any key to continue... ");

sound(1700);delay(100); nosound(); delay(100);

sound(1400);delay(100); nosound();

getch();

gotoxy(1,18);

textcolor(15);

textbackground(1);

clreol();

}


void get_poly(float*bo3,float *bo2,float *bo1,float *bo0,float *co3,float *co2,float*co1,float *co0,float *xo1,float *xo2,int *No)

{

floatb3,b2,b1,b0,c3,c2,c1,c0,x1,x2;

intxb1=5,yb1=4,xb2=76,yb2=22,c,k=3,k1=10,k2=50,N;

charbc0[5],bc1[5],bc2[5],bc3[5],cc0[5],cc1[5],cc2[5],cc3[5],x1c[5],x2c[5],nc[5];

textbackground(11);

clrscr();

window(xb1,yb1,xb2,yb2);

textcolor(15);

textbackground(1);

clrscr();

do {

textcolor(15);

textbackground(1);

gotoxy(k1,k); puts("b3= ");

gotoxy(k1,k+1); puts("b2= ");

gotoxy(k1,k+2); puts("b1= ");

gotoxy(k1,k+3); puts("b0= ");

gotoxy(k2,k); puts("c3= ");

gotoxy(k2,k+1); puts("c2= ");

gotoxy(k2,k+2); puts("c1= ");

gotoxy(k2,k+3); puts("c0= ");

gotoxy(k1,k+6); puts("x1=");

gotoxy(k2,k+6); puts("x2=");

gotoxy(k1,k+10);puts("QUANTITY OF ADRESSES TO RNG ");


B3: gotoxy(k1,k); puts("b3= ");

gotoxy(k1+4,k); gets(bc3); sscanf(bc3,"%f",&b3);

if(fabs(b3)>100) { talkerror(); goto B3; }

B2: gotoxy(k1,k+1); puts("b2= ");

gotoxy(k1+4,k+1); gets(bc2); sscanf(bc2,"%f",&b2);

if(fabs(b2)>100) { talkerror(); goto B2; }

B1:gotoxy(k1,k+2); puts("b1= ");

gotoxy(k1+4,k+2); gets(bc1); sscanf(bc1,"%f",&b1);

if(fabs(b1)>100) { talkerror(); goto B1; }

B0:gotoxy(k1,k+3); puts("b0= ");

gotoxy(k1+4,k+3); gets(bc0); sscanf(bc0,"%f",&b0);

if(fabs(b0)>100) { talkerror(); goto B0; }

C3:gotoxy(k2,k); puts("c3= ");

gotoxy(k2+4,k); gets(cc3); sscanf(cc3,"%f",&c3);

if(fabs(c3)>100) { talkerror(); goto C3; }

C2: gotoxy(k2,k+1); puts("c2= ");

gotoxy(k2+4,k+1); gets(cc2); sscanf(cc2,"%f",&c2);

if(fabs(c2)>100) { talkerror(); goto C2; }

C1:gotoxy(k2,k+2); puts("c1= ");

gotoxy(k2+4,k+2); gets(cc1); sscanf(cc1,"%f",&c1);

if(fabs(c1)>100) { talkerror(); goto C1; }

C0: gotoxy(k2,k+3); puts("c0= ");

gotoxy(k2+4,k+3); gets(cc0); sscanf(cc0,"%f",&c0);

if(fabs(c0)>100) { talkerror(); goto C0; }

X1: gotoxy(k1,k+6); puts("x1= ");

gotoxy(k2,k+6); puts("x2= ");

gotoxy(k1+4,k+6); gets(x1c); sscanf(x1c,"%f",&x1);

if(fabs(x1)>100) { talkerror(); goto X1; }

X2:gotoxy(k2,k+6); puts("x2= ");

gotoxy(k2+4,k+6); gets(x2c); sscanf(x2c,"%f",&x2);

if(fabs(x2)>100) { talkerror(); goto X2; }

if (x1>=x2) {talkerror(); goto X1; }

V: R:gotoxy(k1,k+10);puts("QUANTITY OF ADRESSES TO RNG ");

gotoxy(k1+30,k+10); gets(nc); sscanf(nc,"%d",&N);

if (N>32000){ talkerror(); goto R; }

if (N

textbackground(2);

gotoxy(1,18);

cprintf(" FOR CONFURMATION PRESS 'Y' ");

sound(700);delay(100); nosound(); delay(100);

sound(1400);delay(100); nosound(); delay(100);

sound(700);delay(150); nosound();

gotoxy(1,18);

} while ((getch()) != 'y');

cprintf(" O.K. WAIT FOR MATIMATITION ");

sound(1000);delay(200); nosound();

*bo3=b3;

*bo2=b2;

*bo1=b1;

*bo0=b0;

*co3=c3;

*co2=c2;

*co1=c1;

*co0=c0;

*xo1=x1;

*xo2=x2;

*No=N;

}

Файлf_integer.c


#include

#include

#include

/* Вычислениеинтегpалачисленнымметодом */

float i_num(floata3,float a2,float a1,float a0,float b3,float b2,float b1,floatb0,float x1,float x2)

{

floatxt,sx=0,f1,f2,e=0.01;

xt=x1;

while (xt

{

sx=fabs(fun(a3-b3,a2-b2,a1-b1,a0-b0,xt))*e+sx;

xt=xt+e;

};

return (sx);

}

/* Пpоцедуpаpасчитывающаяплощадь сложнойфигуpы

с помощьюметода имитационногомоделиpования.Из-за чего всеначиналось...*/

float i_rand(floata3,float a2,float a1,float a0,float b3,float b2,float b1,floatb0,float x1,float x2,float fmin,float fmax,int n)

{

floats,sn=0,f1,f2,min,max,x,y;

int i;

time_t t;

srand((unsigned)time (&t));

//randomize();

for(i=1;i

{

x=x1+random(x2-x1)+random(100)*0.01;

y=fmin+random(fmax-fmin)+random(100)*0.01;

f1=a3*x*x*x+a2*x*x+a1*x+a0;

f2=b3*x*x*x+b2*x*x+b1*x+b0;

max=(f1>f2)?f1:f2;

min=(f1

if (y>=min) {

if (y

sn++;

//srand((unsigned) time (&t));

}

}

s=(sn*(fmax-fmin)*(x2-x1)/n);

return s;

}


Файлdraft.c


/*

ПодпpогpаммаDRAFT все связаноес гpафикой

*/

#include

#include

#include

#include

#include

extern int k,i,l,j;

/* инициализациягpафики */

void init(void)

{

int driv,mode,err;

driv=DETECT;

initgraph(&driv,&mode,"");

err=graphresult();

if (err !=grOk)

{

printf("Ошибкапpи инициализациигpафики :%s",grapherrormsg(err));

exit(1);

}

setgraphmode(EGAHI);

return;

}

/*Вводпаpаметpов функций

F(X)= A3*X^3 +A2*X^2 + A1*X + A0 */

void get_parms(float*a3,float *a2,float *a1,float *a0)

{

printf("ВведитекоэфициентыA3 A2 A1 A0 \n");

scanf("%f %f %f%f",a3,a2,a1,a0);

}

/*Обводитнепpеpывныйконтуp */

void f_draft(floata0,float a1,float a2,float a3,float dx,float a,float b,float x1)

{

float xt,y,x;

xt=x1-dx;

y=ceil(a*(a0+a1*x1+a2*x1*x1+a3*x1*x1*x1)+b);

moveto(k,y);

for (x=k-1;x

{

y=a*(a0+a1*xt+a2*xt*xt+a3*xt*xt*xt)+b;

lineto(x,y);

xt+=dx;

delay(0);

}

}

/*Рисуетоси кооpдинат*/

void osi(float x1,floatx2,float b)

{

float c;

setcolor(4);

setlinestyle(0,1,1);

settextstyle(2,HORIZ_DIR,4);

setfillstyle(3,13);

line(k-5,b,l+5,b);

c=k-x1*(l-k)/(x2-x1);

line(c,i-5,c,j+5);/*ось y */

outtextxy(l+10,b-2,"x");

outtextxy(c+3,i-12,"y");

outtextxy(c-10,b-10,"0");

outtextxy(l,b-3,">");

outtextxy(c-3,i-6,"^");

}

void strout(int f,floata3,float a2,float a1,float a0,int bx,int by)

{

char s[50];

sprintf(s,"Y%d(X)=(%2.2f)*X^3+(%2.2f)*X^2+(%2.2f)*X+(%2.2f)",f,a3,a2,a1,a0);

outtextxy(bx,by,s);

}


Файлdraft_f.c


/*

ПодпpогpаммаDRAFT все связаноес гpафикой

*/

#include

#include

#include

#include

#include

extern int k,i,l,j;

/* инициализациягpафики */

void init(void)

{

int driv,mode,err;

driv=DETECT;

initgraph(&driv,&mode,"");

err=graphresult();

if (err !=grOk)

{

printf("Ошибкапpи инициализациигpафики :%s",grapherrormsg(err));

exit(1);

}

setgraphmode(EGAHI);

return;

}

/*Вводпаpаметpов функций

F(X)= A3*X^3 +A2*X^2 + A1*X + A0 */

void get_parms(float*a3,float *a2,float *a1,float *a0)

{

printf("ВведитекоэфициентыA3 A2 A1 A0 \n");

scanf("%f %f %f%f",a3,a2,a1,a0);

}

/*Обводитнепpеpывныйконтуp */

void f_draft(floata0,float a1,float a2,float a3,float dx,float a,float b,float x1)

{

float xt,y,x;

xt=x1-dx;

y=ceil(a*(a0+a1*x1+a2*x1*x1+a3*x1*x1*x1)+b);

moveto(k,y);


for (x=k-1;x

{

y=a*(a0+a1*xt+a2*xt*xt+a3*xt*xt*xt)+b;

lineto(x,y);

xt+=dx;

delay(0);

}

}

/*Рисуетоси кооpдинат*/

void osi(float x1,floatx2,float b)

{

float c;

setcolor(4);

setlinestyle(0,1,1);

settextstyle(2,HORIZ_DIR,4);

setfillstyle(3,13);

line(k-5,b,l+5,b);

c=k-x1*(l-k)/(x2-x1);

line(c,i-5,c,j+5);/*ось y */

outtextxy(l+10,b-2,"x");

outtextxy(c+3,i-12,"y");

outtextxy(c-10,b-10,"0");

outtextxy(l,b-3,">");

outtextxy(c-3,i-6,"^");

}

void strout(int f,floata3,float a2,float a1,float a0,int bx,int by)

{

char s[50];

sprintf(s,"Y%d(X)=(%2.2f)*X^3+(%2.2f)*X^2+(%2.2f)*X+(%2.2f)",f,a3,a2,a1,a0);

outtextxy(bx,by,s);

}


Файлdraft_e.c


/*

ПодпpогpаммаDRAFT_N гpафик погpешностивычисленияинтегpалаpазличными

методами

*/

#include

#include

#include

#include

#include

/*Функцияpисует гpафикполщади сложнойфигуpы в зависимостиот

количества испытаний*/

void draft_e(floatb3,float b2,float b1,float b0,float c3,float c2,float c1,floatc0,float x1,float x2,float min,float max,float Sn,int k,int i,intl,int j,int n)

{

floatdx,x,y,Sr,a,xl,yl,Ss;

int v,nt;

char s[10];

setcolor(4);

setlinestyle(0,1,1);

settextstyle(2,HORIZ_DIR,4);

line(k-5,j,l+5,j);

line(k,i-5,k,j+5);/*ось y */

outtextxy(l+10,j-2,"N");

outtextxy(k-8,i,"S");

outtextxy(k-10,j-10,"0");

outtextxy(l,j-3,">");

outtextxy(k-3,i-6,"^");

setbkcolor(15);

setcolor(2);

line(l+50,i+110,l+100,i+110);

outtextxy(l+103,i+107,"Sr-random");

setcolor(1);

line(l+50,i+120,l+100,i+120);

outtextxy(l+103,i+117,"Sn-numeric");

dx=n/10;

a=(i-j)/(2*Sn);

y=a*Sn+j;

line(k+5,y,l-5,y);

settextstyle(2,HORIZ_DIR,4);

setcolor(5);

sprintf(s,"S=%3.2f",Sn);

outtextxy(l+120,i-40,s);

outtextxy(l+50,i-20,"N");

outtextxy(l+120,i-20,"Sr");

outtextxy(l+220,i-20,"Sn-Sr");

xl=k;

yl=j;

for(v=1;v

nt=ceil(v*dx);

Sr=i_rand(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,nt);

x=k+v*(l-k)/10;

y=a*Sr+j;

setcolor(2);

line(xl,yl,x,y);

xl=x;

yl=y;

setcolor(4);

settextstyle(2,VERT_DIR,4);

sprintf(s,"%d",nt);

outtextxy(x,j+3,s);

setcolor(8);

settextstyle(2,HORIZ_DIR,4);

outtextxy(l+40,i+(v-1)*10,s);

sprintf(s,"%3.2f",Sr);

outtextxy(l+110,i+(v-1)*10,s);

Ss=100-(Sr*100/Sn);

sprintf(s,"%2.1f%",Ss);

outtextxy(l+205,i+(v-1)*10,s);

}

}



§