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

РАiЕТ ПЛОЩАДИ СЛОЖНОЙ ФИГУРЫ МЕТОДОМ ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ

Ва

Разработка программы, позволяющей с помощью метода имитационного моделирования рассчитать площадь сложной фигуры, ограниченной сверху кривой 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

ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа ВаВа 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 fmin,float fmax, ВаВаВаВа ВаВаВаВаВа // минимальное и максимальное значения ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа ВаВаВа ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа //функций на данном интервале

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

Ва

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. Borland International, Inc. 1988.

8.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);

ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа Ва }

}

Ва

Ва Ва Ва Ва Ва Ва Ва Ва Ва Ва Ва Ва Ва

Вместе с этим смотрят:

Расчет размерных цепей
Расчетно-графическая работа по специальным главам математики
Расширения полей. Формальное присоединение элементов
Решение балансовых уравнений