Основы, программирования Fortran, C++, D Санкт-Петербургский государственный университет фак. Прикладной математики – процессов управления гр. 103 – прикладная математика и информатика 2013. 10. 12. Числовые объекты и операции. Занятие № 6 Храмушин Василий Николаевич – V. Khramushin@SPBU. ru каф. компьютерного моделирования и многопроцессорных систем
Шестое занятие 12 октября … Встроенные типы данных Все типы языка Си++: целые числа разной разрядности, вещественные числа, логические величины, перечисляемые значения, символы и их кодировка. Общая информация == Вычислительная математика: ++ комплексные числа и арифметика, Целые числа ++ векторы и тензоры аналитической Вещественные числа геометрии и базовые операции с Логические величины синтезированными вычислительными объектами Символы и байты Кодировка, многобайтовые символы Наборы перечисляемых значений Классный компьютер \ws 101 -operatorShared folderKhram
гр. 103 – Прикладная математика и информатика Власов Александрович + 2. Данилов Артем Дмитриевич + 3. Каюмова Диана Шамильевна + 4. Колесова Александровна -5. Лапин Егор Дмитриевич + 6. Леганькова Дарья Андреевна + 7. Минаева Ульяна Сергеевна + 8. Мозин Виталий Романович + 9. Пенкрат Андрей Олегович + 10. Глотова Екатерина Андреевна 11. Пономарева Юлия Константиновна + 12. Сатина Алена Юрьевна + + 13. Сафронов Вадим Викторович + 14. Симанков Сергей Сергеевич + . Смирнова Дарья Сергеевна 1. + + + + + -+ + + + + + ----+ + + -+ /
Вопросы от 5 октября. Функции: 1. Алгоритмическое и функциональное программирование (реентерабельность, рекурсия, демоны и драйверы…. ); 2. 3. 4. Виды процедур функций в алгоритмических языках программирования, особенности реализации (процедуры, функции, операторы) Соглашения о связях – правила передачи параметров в процедурных блоках программы … (по адресу, по значению? ). … как начало элементарного программирования …
Аппаратные реализации числовых данных Байт без знака Короткое целое число без знака Целое число без знака Длинное целое без знака Число одинарной точности двойной точности увеличенной точности Логическое значение char -128 – 127 unsigned char 0 – 255 short -32768 – 32767 unsigned short 0 – 65535 int -2147483648 – 2147483647 unsigned int 0 – 4294967295 long -2147483648 – 2147483647 unsigned long 0 – 4294967295 float ± 3. 4 e-38 – ± 3. 4 e+38 (7 цифр) double ± 1. 7 e-308 – ± 1. 7 e+308 (15 цифр) long double ± 1. 2 e-4932 – ± 1. 2 e+4932 bool true (истина), false (ложь) Данные в языке С могут именоваться с точностью до каждого бита. Выравнивание любых данных всегда производится на размер машинного слова (32 или 64 бит = 4 или 8 байт)
Форма записи числовых констант //! варианты представления целочисленных констант const int a = 0100; // восьмеричная константа 64 auto int b = 64; // десятичное целое в стеке static int c = 0 x 40; // шестнадцатеричное число register int d = 1<<6; // двоичная 1 в шестой позиции volatile float e=64. 0, f=6. 4 e 1; // плавающее число restrict *char g="+"Уникум". "; // ссылка без алиасов //! собственно числовые величины unsigned [int] h = 32768 u; // 0 х80000000 беззнаковое целое long p = 32768 l; // 0 х80000000 обычное целое x 32 unsigned long q = 65536 ul; // 0 х100 длинное беззнаковое long r = 2147483648 L ; // ? = -2147483648 unsigned long s = 4294967296 UL; // ? = 1 0 х0 float t = 64. 0 F; // по умолчанию все константы - double x = 64. 0; // предпочтительный тип числа long double y=64. 0 L; // х80 причуды сопроцессора IBMx 86
Логические, битовые поля и символы //! двухбайтовые целые числа short n = 32767; // старое целое word – х16 unsigned short j = 65535 u; // fixed – символ UTF-8 struct Time. And. Date // битовое поле (? ? всего 39 бит – 5 байт) { unsigned unsigned }; hours mins secs week. Day month year : 5; : 6; : 3; : 6; : 5; : 8; // часы от 0 до 24 // минуты // секунды от 0 до 60 // день недели // день месяца от 1 до 31 // месяц от 1 до 12 // год от 0 до 100 //! логические и символьные bool k = true; // true = 1, false = 0, – двоичная логика char l = 64; // signed char -128 – байт byte m = 64; // unsigned char 0 – 255 – байт // GCC – для русских программ очень важно: -funsigned-char
Алфавит и кодировки //! Кодировки шрифтов: DOS-866, Win-1251, UTF-8 char // unsigned char 0 – 255 – байт a b f n r t v ' " \ звонок возврат на один символ назад перевод страницы новая строка перевод каретки горизонтальная табуляция вертикальная табуляция апостроф двойные кавычки обратная дробная черта char zero = ' '; unsigned char bit = ' x. FF'; const char tab = ' 10'; wchar_t utf = L‘ab'; char str[] = "просто строка";
Перечисления и наименования
Базовые числовые объекты и операции
Оформление письменного задания (титульный лист) Правительство Российской Федерации Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Санкт-Петербургский государственный университет» Каюмова Диана Шамильевна Компиляция и интерпретация (достоинства в использовании и проблемы реализации формализованных и специальных языков программирования ) Самостоятельная работа по курсу «Основы программирования» Специальность: 010501. 65 «Прикладная математика и информатика» Научный руководитель: к. т. н. Храмушин В. Н. ______ Студентка гр. 103. Каюмова Д. Ш. ______ Санкт-Петербург 2013
Содержание документа ГОСТ 7. 32 -2001 «Отчет о НИР» , Минск, 2001 Структурными элементами отчета о НИР являются: - титульный лист; - список исполнителей; - реферат; - содержание; - нормативные ссылки; - определения; - обозначения и сокращения; - введение; - основная часть; - заключение; - список источников; - приложения. Общепринятая форма научного творчества: 1 стр. Титульный лист; 2 стр. Расширенное название, аннотация (1 абзац), ключевые слова 3 стр. Реферат, список обозначений. 4 стр. Введение, постановка задачи 5 стр. Основная часть 3 -4 с. Последняя стр. Заключение, список использованных источников, приложения. Не печатать! Выслать на V. Khramushin@SPBU. ru
14 сентября. Основы: Вопросы Концепция построения и особенности реализации языков программирования на примерах Fortran(57), Algol(58), В(69), C(72), C++(83), D(2001) и др; 2. Машинно-зависимые типы данных (стек, статика, куча), особенности их использования, скорость доступа и эффективность выполнения операций на линейном процессоре. 3. Компиляция и интерпретация, достоинства в использовании и проблемы реализации формализованных и специальных языков программирования (декларативные алгоритмы и функциональные языки 1. искусственного интеллекта) 21 сентября. Операции: 4. Ключевые слова и зарезервированные идентификаторы языков программирования (объявление данных и исполняемые операторы – Алгол – Фортран, С, С++, Ди); Математические и логические операции исторических языков программирования (математические и программные выражения) 6. Приоритеты и последовательность константных и исполняемых выражений (литералы времени компиляции и 5. интерпретация выражений на этапах сборки и исполнения программы). … только суть – что и зачем?
Вопросы 28 сентября. Операторы: 7. Структура типовой программы на алгоритмическом языке программирования (операторы описания данных и выполняемые команды); 8. Типовое построение среды компиляции и сборки программ на алгоритмических языках программирования (структура библиотек на примере проекта GCC. GNU) 9. Реализация простейшей интерактивной графики в среде GLFW и Glut для Open. GL (особенности исполнения программы, управляемой внешними прерываниями). … и то уже не просто … 5 октября. Функции: 10. Алгоритмическое и функциональное программирование (реентерабельность, рекурсия, демоны и драйверы…. ); 11. Виды процедур функций в алгоритмических языках программирования, особенности реализации (процедуры, функции, операторы) 12. Соглашения о связях – правила передачи параметров в процедурных блоках программы … (по адресу, по значению? ). 13. … начало элементарного программирования …
Вопросы 12 октября: «Данные» : 13. Особенности реализации числовых данных для типовой архитектуры арифметико-логического устройства и организации памяти компьютера (логические, целые, беззнаковые, плавающие, х16 -32 -64 числа); 14. Форма записи числовых констант с вариантами уточнения их аппаратного представления (восьмеричные, десятичные, шестнадцатеричные целые, битовые поля или плавающие числа с фиксированной точкой и экспонентой) 15. Алфавит, символьные кодировки и варианты именованных перечислений и макроопределений (шрифты Dos, Win и UTF, последовательности и макроопределения). 16. … далее объектно-ориентированное программирование …
Седьмое занятие 19 октября … Числовые объекты и операции Способы описания классов. Создание объектов. Обращение к атрибутам и методам объектов. == Вычислительная Понятие класса математика: Определение методов класса ++ векторы и тензоры аналитической Переопределение операций геометрии и базовые операции с Подписи методов синтезированными и необязательные аргументы вычислительными объектами Запись классов
Практика, 12 октября, 2 -3 пара Персональные компьютеры CBcodeblocks. exe – интегрированная среда разработки Min. GW – компиляторы Fortran, C, C++, (D) Min. GWbin – исполняемые модули (*. exe, *. dll) Min. GWinclude – прототипы стандартных пакетов (*. h) Min. GWlib – библиотеки программ Win 32 и Сlib (*. а) Классный компьютер login: student 101 e, password: student 101 e \ws 101 -operatorShared folderKhramGL
Элементарные числовые объекты логические конструкции этапа компиляции С++ Скаляры: индексы: int; числа: real; typedef double real; A = R + r·a a = a·r =(A-R)·r Векторы: точки: point struct{ real x, y, z; }; векторы: vector struct{ point a, v; }; Тензоры: матрицы: tensor struct{ point x, y, z; }; базисы: basis struct{ point a; tensor m; }; Алгебраические структуры: группы; кольца; поля - самостоятельно
Тензорная арифметика A = R + r·a a = a·r =(A-R)·r Вектор Тензор Вектор r = r – 1 – дуальность В исчислении флюксий Ньютона: Выражение для скорости в следующий момент времени – тензор v=v·r конвективных скоростей (живых сил). Уравнение движения: В численных методах механики необходимо вводить смешанные тензоры, и – тензорную математику
Соглашения о «связях» тензорной арифметики Прописными буквами обозначаются величины, измеренные относительно центра единой абсолютной системы координат; Строчными буквами – в локальном базисе, связанном с подвижной и деформируемой частицей сплошной среды. Основные обозначения: A = Ai – пространственная точка (вектор), измеренная и отмеченная в глобальной системе координат [м]; a = ak – векторный отсчет внутри локального базиса элементарной частицы жидкости [м-2]; r = rk = rik – тензор формы крупной частицы жидкости [м 3]; r = r j = r kj – тензор плотности [кг·м-3]; М = м ij = r · r – смешанный тензор, соотносящий внутреннее состояние частицы на абсолютную систему отсчета [кг]; v = v · r – тензор конвективных скоростей, характеризующий внутреннюю кинематику (живые силы) в локальном базисе расчетного континуума [сек]. (= длительности или частоты-1)
Вычислительные объекты и операции Объектно-ориентированные средства С++ #include <Math. h> // // Стандартный математический пакет - для скалярных числовых величин Новый вычислительный объект - пространственная точка class Point { real x, y, z; public: Point(): x(0), y(0), z(0){} // конструкторы Point( real _x, real _y=0, real _z=0 ) : x( _x ), y( _y ), z( _z ){} Point operator - (){ x=-x, y=-y, z=-z; return *this; } Point operator +=( Point P ); // Point operator -=( Point P ); // определение вычислительных Point operator *=( real W ); // операций и преобразований Point operator /=( real W ); // friend Point operator + ( Point P, Point Q ){ return P+=Q; } friend Point operator - ( Point P, Point Q ){ return P-=Q; } friend real abs( Point P ) { return sqrt( P. x*P. x + P. x*P. x ) }. . . }; inline Point: : operator +=( Point P ) { x+=P. x, y=P. y, z=P. z; return *this; }
Вычислительные объекты и операции определение функций и переопределение операций С++ #include ". . MathPoint. h" // собственный пакет для векторов // // Сложение векторов // Point: : operator +=( Point P ) { x+=P. x, y=P. y, z=P. z; return *this; } // // удлинение модуля вектора // Point: : operator += (real W ) { W *= W/sqr( *this )+1. 0; x*=W; y*=W; z*=W; } // // возведение в квадрат // inline double sqr( real w ){ return w*w; } inline double sqr( Point P ) { return P. x*P. x + P. y*P. y + P. z*P. z; } // // ошибочное переопределение, которое компилятор не заметит Point: : operator+( Point P ) { x+=P. x, y+=P. y, z+=P. z; return *this; }
Point. h 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 #ifndef __Point_HXX #define __Point_HXX #include <Math. h> // // Представление точки в двумерном пространстве (Point) // и в дискретном поле графического монитора (point) // struct Point; // struct point // Точка с дискретным отсчетом { // int x, y; operator Point(); point& operator+=( const point& p ){ x+=p. x, y+=p. y; return *this; point& operator-=( const point& p ){ x-=p. x, y-=p. y; return *this; friend bool operator==( const point& p, const point& q ) { return p. x==q. x && p. y==q. y; friend bool operator!=( const point& p, const point& q ) { return p. x!=q. x || p. y!=q. y; friend point operator+( point p, const point& q ){ return p+=q; } friend pint operator-( point p, const point& q ){ return p-=q; } }; } }
Point. h #define Ret( X, Y ){ Point _P_={ float( X ), float( Y ) }; return _P_; } #define Refer const Point& 24 // 25 struct Point // Представление точки на плоскости 26 { // 27 float x, y; operator point(); 28 29 friend Point operator+( Point P, Refer Q ){ return P+=Q; } 30 friend Point operator-( Point P, Refer Q ){ return P-=Q; } 31 friend Point operator*( Point P, Refer Q ){ return P*=Q; } 32 friend Point operator/( Point P, Refer Q ){ return P/=Q; } 33 friend Point operator+( Point P, float D ){ P. x+=D; return P; } 34 friend Point operator-( Point P, float D ){ P. x-=D; return P; } 35 friend Point operator*( Point P, float D ){ return P*=D; } 36 friend Point operator*( float D, Point P ){ return P*=D; } 37 friend Point operator/( Point P, float D ){ return P/=D; } 38 39 friend bool operator==(Refer A, Refer B){return A. x==B. x&&A. y==B. y; } 40 friend bool operator!=(Refer A, Refer B){return A. x!=B. x||A. y!=B. y; }
Point. h Point& Point& Point& set( float X ){ x=X, y=0. 0; return *this; } set( float X, float Y ){ x=X, y=Y; return *this; } operator -()Ret( -x, -y ) operator ~()Ret( y, x ) operator =( float X ){ x=X, y=0. 0; return *this; } operator+=( Refer P ){ x+=P. x, y+=P. y; return *this; } operator-=( Refer P ){ x-=P. x, y-=P. y; return *this; } operator*=( float A ){ x*=A, y*=A; return *this; } operator/=( float A ){ x/=A, y/=A; return *this; } operator*=( Refer B ){ set( x*B. x-y*B. y, x*B. y+y*B. x ); return *this; } operator/=( Refer B ){ return *this *= conj( B )/norm( B ); } friend Point conj( Refer P ) Ret( P. x, -P. y ) friend float norm( Refer P ){ return P. x*P. x+P. y*P. y; } friend Point sqr( Refer P ) Ret( P. x*P. x-P. y*P. y, 2*P. x*P. y ) friend float abs( Refer P ){ return hypot( P. y, P. x ); } friend Point dir( Point P ){ return P/=abs( P ); } friend float arg( Refer P ){ return atan 2( P. y, P. x ); } friend Point exp( Refer P ); //friend Point log( Refer P )Ret( log( abs( P ) ), arg( P ) ) //friend Point pow( Refer P, Refer pw ){ return exp( log( P )*pw ); } //friend Point sqrt( Refer P ){ static Point W={ 0. 5, 0 }; return pow( P, W ); } }; inline Point set( float X )Ret( X, 0 ) inline Point set( float X, float Y )Ret( X, Y ) inline float sqr( float R ){ return R*R; } inline float abs( float A ){ return fabs( A ); } inline float norm( float X, float Y ){ return X*X+Y*Y; } inline Point polar( float angle )Ret( cos( angle ), sin( angle ) ) inline Point polar( float mag, float angle )Ret( mag*cos(angle), mag*sin(angle) ) inline Point exp( Refer P ){ return polar( exp( P. x ), P. y ); } //{ P. x=exp( P. x ); return P. set( P. x*cos( P. y ), P. x*sin( P. y ) ); } #undef Refer #undef Ret #endif
Вычислительные объекты и операции Определения и операции с тензорами в С++ (пакет Tv_Graphics. h) // // BASIS. hxx Ver 1. 0 90. 11. 20 // // Процедуры построения альтернативного геометрического Базиса // и базиса по четырехугольной деформированной площадке // // ¦ x. X y. X ¦ ¦ Xx Yx ¦ Y y Def // M = ¦ ¦ m = ¦ ¦ T A --¬ // ¦ x. Y y. Y ¦ ¦ Xy Yy ¦ ¦ ¦ A ¬a -¦ // ¦ ¦ / v // -1 conjg( M ) ¦ y. Y -y. X ¦ ¦ ¦ / // m = M = ----- = ¦ ¦ / ¦M¦ ¦ ¦ / // det( M ) ¦-x. Y x. X ¦ ¦ ¦ r // ¦ +------>x // A = Re + a · M ¦ Re // ¦ // a = ( A-Re ) / M L=========>X // #ifndef _Basis_Hxx #define _Basis_Hxx #include ". . Type. h" + беглый обзор пакета 2 -D процедур
Вычислительные объекты и операции Определения и операции с тензорами в С++ // // // // // // // B A S I S - Определение локального базиса внутри глобального Re y, x size angle - координаты центра локального базиса определяющие векторы длины осей ортогонального базиса ориентация меридиана (оси Y) в градусах Axis( Re, y, x ) - Явное определение линейного базиса с помощью граничных векторов { y, x } Decart( Re, x ) - Здесь вектор { x } задает не более чем направление, что исключает (смысл и) геометрические размерения для определяемого таким образом базиса Euler( Re, angle, size ) - 'Географический базис заданный направлением и размерами прямоугольной области Здесь под направлением понимается угол в °(гр) от Северного Меридиана, отсчитываемый по ходу Часовой стрелки Area Прямое определение по опорным точкам c 2 ( Aim - Результирующая прямоугольная площадка LuЇ Ru Crs - Тройка угловых точек, просматриваемых как: ¦Crs ) - на выходе площадь под угловыми точками Ld+---+Rd Ld - перенос при сохранении деформации Aim c 1 Ld, Rd - устанавливает ортогональную систему координат Ld, Rd, Lu - это еще линейный базис Ld, Rd, Lu, Ru – билинейная аппроксимация локального поля
Вычислительные объекты и операции Определения и операции с тензорами в С++ // struct Base // Линейные преобразования в единичном базисе { // Point y, x, // Вид локального базиса в глобальном Y, X; // Отображение главных осей Base(); // float Axis( Point Y, Point X ); // Построение физического базиса Point Goint( Point A ); //{ A. x*X. x+A. y*Y. x, Point Goext( Point a ); //{ a. x*x. x+a. y*y. x, A. y*Y. y+A. x*X. y; } a. y*y. y+a. x*x. y; } }; // struct Basis: Base // Линейные преобразования с глобальным отсчетом { // Point R; // Глобальный Центр отсчета Basis(); // float Axis ( Point Re, Point Y, Point X ); float Euler ( Point Re, float Angle, Point size ); float Decart( Point Re, Point x ); Point Goint( Point A ); //{ return Base: : Goint( A-R ); } Point Goext( Point a ); //{ return Base: : Goext( a )+R; } }; //
Вычислительные объекты и операции Определения и операции с тензорами в С++ // struct Basis_Area: Basis // { // Basis S; // Field B; // float ky, kx; // Basis_Area(); float Point Не. Линейные - с внутренними размерениями Встречный базис в противоположном углу Локальный Центр начала отсчета с размерами Деформация до произвольного четырехугольника // luЇ // ¦Crs[3] Area( Field Aim, Point Ld, Point Rd ); // ld+---¬rd Area( Field Aim, Point Ld, Point Rd, Point Lu ); Area( Field Aim, Point Ld, Point Rd, Point Lu, Point Ru ); // Goint( Point& A ); // Переход в локальный базис Aim Goext( Point& a ); // Возврат в деформированный базис Crs // }; // // Поиск точки пересечения двух прямых, даже если таковые не сходятся // int Cross( Point L 11, Point L 12, Point L 21, Point L 22, Point& Cx ); // #endif
Вычислительные объекты и операции Процедуры тензорно-векторных преобразований в С++ float Base: : Axis( Point _Y, Point _X ) { y=_Y; // x=_X; // Axis - Определение локального базиса float det=y. y*x. x - y. x*x. y; // Y. y=x. x/det; Y. x=-y. x/det; // Обращение матрицы переходов X. x=y. y/det; X. y=-x. y/det; return det; } static Point _P; Point Base: : Goint( Point A ) { _P. x=A. x*X. x+A. y*Y. x; _P. y=A. y*Y. y+A. x*X. y; return _P; } Point Base: : Goext( Point a ) { _P. x=a. x*x. x+a. y*y. x; _P. y=a. y*y. y+a. x*x. y; return _P; } // // Тензорное преобразование и переносом начала осей координат // Point Basis: : Goint( Point A ){ return Base: : Goint( A-R ); } Point Basis: : Goext( Point a ){ return Base: : Goext( a )+R; } float Basis: : Axis( Point Re, Point _Y, Point _X ) { R=Re; return Base: : Axis( _Y, _X ); }
Вычислительные объекты и операции Процедуры тензорно-векторных преобразований в С++ // Здесь добавляется возможность деформации базиса, ( 93. 11. 23 ) // что не желательно нагружать на Basis. // Area // ( Aim - Результирующая прямоугольная площадка LuЇ // Corners - Тройка угловых точки просматриваемых как: ¦Crs[3] // ) - на выходе площадь под угловыми точками Ld+---¬Rd // Aim // Area Lu Ru // ( Aim - Результирующая прямоугольная площадка 3 -----¬ 4 // Corners - Четыре угловых точки просматриваемых как: ¦ Cs ¦ // ) - на выходе площадь под угловыми точками 1 L-----2 // Ld Rd float Basis_Area: : Area( Field Aim, Point Ld, Point Rd, Point Lu ) { Point P; P. x=Rd. x; P. y=Lu. y; return Area( Aim, Ld, Rd, Lu, P ); } float Basis_Area: : Area( Field Aim, Point Ld, Point Rd ) { return Area( Aim, Ld, Rd, Ld+(Rd-Ld)*set(0, 1)*Aim. Ly/Aim. Lx ); } float Basis_Area: : Area( Field Aim, Point Ld, Point Rd, Point Lu, Point Ru ) { float A=Basis: : Axis( Ld, Lu-Ld, Rd-Ld ); B=Aim; Point D=Basis: : Goint( Ru ); kx=D. x-1; ky=D. y-1; return ( A+S. Axis( Ru, Rd-Ru, Lu-Ru ) )/2. 0; }
Вычислительные объекты и операции Процедуры тензорно-векторных преобразований в С++ Point Basis_Area: : Goint( Point& P ) // Переход в локальный базис Aim { // Point A=S. Goint( P ), a=Basis: : Goint( P ); A. x=B. Jx+B. Lx*( a. x/( A. x+a. x ) ); A. y=B. Jy+B. Ly*( a. y/( A. y+a. y ) ); return A; } // Point Basis_Area: : Goext( Point& p ) // Возврат в центральный базис { // It's Right !!! Point A, a; A. x=( p. x-B. Jx )/B. Lx; // A. y=( p. y-B. Jy )/B. Ly; a. x=A. x*( 1+kx*A. y ); a. y=A. y*( 1+ky*A. x ); return Basis: : Goext( a ); } // // Поиск точки пересечения двух прямых линий, // даже если таковые не сходятся // int Cross( Point L 11, Point L 12, Point L 21, Point L 22, Point& Cx ) { Base Det; if( abs(Det. Axis( L 22 -L 21, L 12 -L 11 ))<1. e-12 )return -2; // -- линии L 21=L 22=Det. Goint( L 22 -L 11 ); // почти параллельны L 21. y=0; Cx=Det. Goext( L 21 )+L 11; if( L 22. x<0. 0 || L 22. y<0. 0 )return -1; if( L 22. x>1. 0 || L 22. y>1. 0 )return 1; return 0; }
Скачать презентацию Основы программирования Fortran C D Санкт-Петербургский государственный университет
Bpl 2013-10-12 - Данные.ppt