П’ята лабораторна робота
x = [0:0.1:7];
f = exp(-x);
subplot(1, 2, 1)
plot(x, f)
g = sin(x);
subplot(1, 2, 2)
plot(x, g)
Створюю м файл
function f = myfun(x)
f=exp(-x).*sqrt((x.^2+1)/(x.^4+0.1));
>> x = [0:0.5:4];
>> y= myfun(x);
>> plot(x, y)
>> fplot('myfun', [0 4])
Створюю м-файл
function x = root2(a,b,c)
% розраховує (ax^2+bx+c)
% виволить попередження якщо результат комплексний
D = b^2-4*a*c;
if D<0
warning('результат комплексний')
end
x = (-b+sqrt(D))/4*a
x = (-b-sqrt(D))/4*a;
>> x = root2(1,4,4)
x =
-1
x =
-1
Створюю м-файл
function z = NOD(x, y)
% знаходить найбільший спільний дільник
while(x~=y)
if (x>y)
x=x-y;
else y=y-x;
end
end
z=y;
>> z=NOD(81, 36)
z =
9
Створюю м-файл
function q = prost(N)
%знаходить прості числа
i=1;
for n= 1:1:(150+10*N)
q(:,i)=n;
i=i+1;
end
p=2;
c=2;
while ((p^2)<=(150+10*N))
while (p<=(150+10*N))
q(:,p+c)=0;
p=p+c;
end
c=c+1;
p=c;
end
end
>> q=prost(13)
q =
Columns 1 through 9
1 2 3 0 5 0 7 0 0
Columns 10 through 18
0 11 0 13 0 0 0 17 0
Columns 19 through 27
19 0 0 0 23 0 0 0 0
Columns 28 through 36
0 29 0 31 0 0 0 0 0
Columns 37 through 45
37 0 0 0 41 0 43 0 0
Columns 46 through 54
0 47 0 0 0 0 0 53 0
Columns 55 through 63
0 0 0 0 59 0 61 0 0
Columns 64 through 72
0 0 0 67 0 0 0 71 0
Columns 73 through 81
73 0 0 0 0 0 79 0 0
Columns 82 through 90
0 83 0 0 0 0 0 89 0
Columns 91 through 99
0 0 0 0 0 0 97 0 0
Columns 100 through 108
0 101 0 103 0 0 0 107 0
Columns 109 through 117
109 0 0 0 113 0 0 0 0
Columns 118 through 126
0 0 0 0 0 0 0 0 0
Columns 127 through 135
127 0 0 0 131 0 0 0 0
Columns 136 through 144
0 137 0 139 0 0 0 0 0
Columns 145 through 153
0 0 0 0 149 0 151 0 0
Columns 154 through 162
0 0 0 157 0 0 0 0 0
Columns 163 through 171
163 0 0 0 167 0 0 0 0
Columns 172 through 180
0 173 0 0 0 0 0 179 0
Columns 181 through 189
181 0 0 0 0 0 0 0 0
Columns 190 through 198
0 191 0 193 0 0 0 197 0
Columns 199 through 207
199 0 0 0 0 0 0 0 0
Columns 208 through 216
0 0 0 211 0 0 0 0 0
Columns 217 through 225
0 0 0 0 0 0 223 0 0
Columns 226 through 234
0 227 0 229 0 0 0 233 0
Columns 235 through 243
0 0 0 0 239 0 241 0 0
Columns 244 through 252
0 0 0 0 0 0 0 251 0
Columns 253 through 261
0 0 0 0 257 0 0 0 0
Columns 262 through 270
0 263 0 0 0 0 0 269 0
Columns 271 through 279
271 0 0 0 0 0 277 0 0
Columns 280 through 288
0 0 0 0 0 0 0 0 0
Columns 289 through 294
0 0 0 0 0 0
Створюю м-файл
function s = sumaprost(N)
%знаходить суму простих чисел
s=0;
i=1;
for n= 1:1:(150+10*N)
q(:,i)=n;
i=i+1;
end
p=2;
c=2;
while ((p^2)<=(150+10*N))
while (p<=(150+10*N))
q(:,p+c)=0;
p=p+c;
end
c=c+1;
p=c;
end
s=sum(q);
end
>> s=sumaprost(13)
s =
7419