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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Файл WINDOW. C

– ввод параметров;

Процедура get_poly

– сообщение об ошибке при вводе;

Процедура talkerror

– рисование рамки окна;

Процедура border

Файл MATIM. C

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

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

Процедура f_max

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

Заданной точке;

Процедура fun

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

Уравнения;

Процедура f_root

Файл F_INTEGER. C

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

Методом;

Процедура i_num

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

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

Процедура i_rand

Файл DRAFT. C

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

Процедура init

– обводка непрерывного контура

Процедура 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 Список использованной литературы.

1. Язык программирования Си для персонального компьютера. С. О. Бочков, Д. М. Субботин.

2. С++ . Описание языка программирования. Бьярн Страустрап.

3. TURBO C. User’s Guide. Borland International, Inc. 1988.

.4. TURBO C. Reference Guide. Borland International, Inc. 1988.

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

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

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

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

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

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

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

Файл sq. c

/*

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

*/

#include <stdio. h>

#include <graphics. h>

#include <conio. h>

#include <math. h>

#include <stdlib. h>

#include “matim. c”

#include “window. c”

#include “f_integr. c”

#include “draft. c”

#include “draft_f. c”

#include “draft_e. c”

Int k=20,i=15,l=270,j=140;

Void main(void)

{

Float b0,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<minc)?minb:minc;

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 <stdio. h>

#include <conio. h>

#include <math. h>

#include <stdlib. h>

#include <dos. h>

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

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

Void f_max(float a3,float a2,float a1,float a0,float x1,float x2,float *amin, float *amax)

{

Float dx, 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)

Fmin=Fx;

} while ( x<x2 );

*amin=Fmin;

*amax=Fmax;

}

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

Int f_root(float a0,float a1,float a2,float a3,float x1,float x2,float e, float *k1,float *k2,float *k3)

{ float ku1,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<0)

{ 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<x2);

*k1=ku1;

*k2=ku2;

*k3=ku3;

Return c;

}

Float fun(float a3,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 <stdio. h>

#include <conio. h>

#include <math. h>

#include <stdlib. h>

#include <dos. h>

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

Void border(int sx, int sy, int ex, int ey){

Int i;

For (i=sx+1;i<ex;i++){ //гоpизонтальные линии

Gotoxy(i, sy);

Putch(205);

Gotoxy(i, ey);

Putch(205);

}

For (i=sy+1;i<ey;i++){ //веpтикальные линии

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)

{

Float b3,b2,b1,b0,c3,c2,c1,c0,x1,x2;

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

Char bc0[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<1) { talkerror(); goto V; }

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 <stdio. h>

#include <math. h>

#include <time. h>

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

Float i_num(float a3,float a2,float a1,float a0,float b3,float b2,float b1,float b0,float x1,float x2)

{

Float xt, sx=0,f1,f2,e=0.01;

Xt=x1;

While (xt<x2)

{

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(float a3,float a2,float a1,float a0,float b3,float b2,float b1,float b0,float x1,float x2,float fmin, float fmax, int n)

{

Float s, sn=0,f1,f2,min, max, x, y;

Int i;

Time_t t;

Srand((unsigned) time (&;t));

//randomize();

For(i=1;i<n;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<f2)?f1:f2;

If (y>=min) {

If (y<=max)

Sn++;

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

}

}

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

Return s;

}

Файл draft. c

/*

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

*/

#include <stdio. h>

#include <graphics. h>

#include <conio. h>

#include <math. h>

#include <stdlib. h>

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(float a0,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<l+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,float x2,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, float a3,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 <stdio. h>

#include <graphics. h>

#include <conio. h>

#include <math. h>

#include <stdlib. h>

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(float a0,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<l+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,float x2,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, float a3,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 <stdio. h>

#include <graphics. h>

#include <conio. h>

#include <math. h>

#include <stdlib. h>

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

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

Void draft_e(float b3,float b2,float b1,float b0,float c3,float c2,float c1,float c0,float x1,float x2,float min, float max, float Sn, int k, int i, int l, int j, int n)

{

Float dx, 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<11;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);

}

}


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