Page 103 - 363_
P. 103
104
Для цього слід активізувати розділ "File" головного меню командного
вікна MatLAB, обрати у підменю, що з'явиться на екрані, розділ "New", а потім
- команду "M-file". На екрані з'явиться вікно текстового редактора. В ньому
слід набрати наступний текст:
function у = F1(x, d)
% Процедура, що обчислює значення функції
% у = (d^3)*ctg(x)*sqrt(sin(x)^4-cos(x)^4).
% Звернення у = F1(x, d).
у = (d^3)*cot(x).*sqrt(sin(x).^4-cos(x).^4);
Після цього необхідно зберегти цей текст під ім'ям “Fl.m”. Необхідний
М-файл утворено.
Відтепер можна користуватися цією функцією при розрахунках. Так,
якщо ввести команду
>> у = F1(1, 0.1),
то одержимо результат
y = 4.1421e-004.
Слід зауважити, що аналогічно можна одержати одразу вектор усіх
значень означеної функції при різних значеннях аргументу, якщо останні
зібрати у деякий вектор. Так, якщо сформувати вектор:
>> zet= 0 : 0.3 :1.8;
і звернутися до цієї ж процедури
>> my = Fl(zet,l),
то одержимо:
Warning: Divide by zero
my =
Columns 1 through 4
NaN+ Infi 0+2.9369i 0+0.8799i 0.3783
Columns 5 through 7
0.3339 0.0706 -0.2209
Примітки.
1. Можливість використання сформованої процедури як для
окремих чисел, так і для векторів і матриць зумовлено застосуванням у
запису відповідного М-файлу замість звичайних знаків арифметичних дій
їх аналогів з попередньою крапкою.
2. Щоб уникнути виведення на екран небажаних проміжних
результатів, необхідно у тексті процедури усі обчислювальні оператори
завершувати символом ";".