МОСКОВСКИЙГОСУДАРСТВЕННЫЙИНСТИТУТЭЛЕКТРОННОЙТЕХНИКИ
(ТЕХНИЧЕСКИЙУНИВЕРСИТЕТ)
Расчет площади сложной фигурыс помощью метода
имитацеонногомоделирования.
ЛогвиненкоВ.
Москва. 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); } } §