й довідник.- СПб .: Питер, 2002. - 528 с.
Додаток 1
Протокол програми (в М-файлі)
format long
% метод трапеції=1=1=1=3=100=abs (xn-x0)/n=int (x0, a, b); i=1: n - 1
x (i)=x0 + i * h;
s=s + 2 * int (x (i), a, b); ( метод трапецій )=(s + int (xn, a, b)) * h/2
disp ( аналітичні рішення )=intt (xn, a, b) -intt (x0, a, b)
% метод Симсона
j=0;=int (x0, a, b); j lt; n - 1=j + 1;=s + 4 * int (x0 + j * h, a, b); j == n - 1
continue=j + 1;=s + 2 * int (x0 + j * h, a, b); ( метод Симсона )=h * (s + int (xn, a , b))/3
t=[0.06943184 0.33000984 0.66999052 0.93056816];=[0.17392742 0.32607258 0.32607258 0.17392742];=0; i=1: 4
s=s + w (i) * int (x0 + (xn-x0) * t (i), a, b);
end ( квадратурні формули Гаусса )
s=(xn-x0) * s
tch=[0.083751 0.312730 0.5 0.68727 0.916249];=0; i=1: 5
s=s + 0.2 * int (x0 + (xn-x0) * tch (i), a, b); ( метод Чебишева )=(xn-x0) * s
function y=int (x, a, b)=1/((a * x + b) ^ 2); y=intt (x, a, b)=- 1/(a * (a * x + b));
Додаток 2
Протокол програми (в М-файлі) long=1
% кількість кроків=1000;
% обнуленіемассів=zeros (1, n);
y=zeros (1, n);
% н.у
x (1)=0; (1)=0.1;
% крок рішення
h=(xend-x (1))/n; i=2: n
k1=odddde (x (i - 1), y (i - 1));
k2=odddde (x (i - 1) + h/2, y (i - 1) + h * k1/2);
k3=odddde (x (i - 1) + h/2, y (i - 1) + h * k2/2);
k4=odddde (x (i - 1) + h, y (i - 1) + h * k3);
y (i)=y (i - 1) + h * (k1 + 2 * k2 + 2 * k3 + k4)/6;
x (i)=x (i - 1) + h;
% plot (x, y)=0.1;
[xm, ym]=ode45 (odddde, [0: 0.001: 1], y0); (x, y, xm, ym)
Додаток 3
Протокол програми (в М-файлі)
% Введення рядками=[2 - 4 3 1;
- +1 5 - 7 - 3;
10 - 2 4 4;
- 1 1 - 1 - 1];=[7;- 24; 34;- 6];=a b
% Введення стовпцями=[2 - 1 10 - 1;
- 4 +5 - 2 січня;
3 - 7 +4 - 1;
1 - 3 +4 - 1];=[7 - 24 34 - 6];=b/a=b * a ^ - 1=b * inv (a)
Додаток 4
Протокол програми (в М-файлі)
% початкове значення кореня=1;
% лічильник інтерація
n=0; abs (ur (x)/dur (x)) gt; 1e - 7
x=x-ur (x)/dur (x);
n=n + 1;
if n gt;=1000
disp ( рішення немає )
break
end(x)(n)y=ur(x)=x^4-5*x^3+8*x^2-5*x+1;y=dur(x)long=1;=7;(1)=x-h;(1)=ur(x(1));i=1:n;(i+1)=x(1)+i*h;(i+1)=ur(x(i+1));(x)(y)i=1:n- 1 (i)=y (i + 1) -y (i); i=1: n - 2y (i)=dy (i + 1) -dy (i); i=1: n - 3y (i )=d2y (i + 1) -d2y (i); i=1: n - 4y (i)=d3y (i + 1) -d3y (i);=(dy (1) + 0.5 * d2y (1)-d3y (1)/6 + d4y (1)/12)/h;