x=-2:0.001:0.00001;
x2=0:0.01:2;
y1=sin(1./x);
y2=0;
plot(x,y1,x2,y2);

-----------------------------

more off
suma = 0;
n = 0;
do
  wyraz = 2.^n ./ factorial(n);
  suma = suma + wyraz;
  n = n + 1;
until wyraz < 0.01
n = n - 1
suma

-------------------------------

n=1;
wyraz=0;
suma=1;
do
  wyraz = 1./(n.*sqrt(n));
  suma = suma + wyraz;
  n = n+1;
until wyraz < 0.01
n = n-1
suma
wyraz

---------------------------------

more off
z = magic(7);
[n,s]=size(z);
sumka = 0;
if (rem(n,2) == 0)
  for (ii = 1 : n)
    sumka = sumka + z(ii,ii) + z(ii,n-ii+1);
  endfor
endif
if (rem(n,2) ~= 0)
  for (ii = 1 : n)
    sumka = sumka + z(ii,ii) + z(ii,n-ii+1);
  endfor
  x = ceil(n./2);
  sumka = sumka - z(x,x);
endif
sumka

----------------------------------

x1=linspace(2,10,12);
y3=sin(4*x1)./x1;
q=trapz(x1,y3)

----------------------------------

%zad 1
A=[];
suma=0;
for ii=1:10
  suma=suma+ii;
  A(ii)=suma;
  endfor

  %zad 2

 N=[1 2 -1 3
   -3 1 0 1
   2 1 1 0
   1 -1 1 1];
   W=[7
      0
      7
      1];
    wynik=N\W;

 %zad 3
 tol=3;
 xp=1;
 xk=2;
 Q=quad('calka',xp,xk,tol);
 Q
 %zad 4
 n=5;
 v=silnia(n);
 v

 function [v]=silnia(n)
  sil=1;
  for ii=1:1:n
    sil=sil*ii;
    endfor
    v=sil;

------------------------

more off
A = [1 1 1; 1 1 1];
B = [2; 2];
C = [3 3 3 3];
D = [A B
     C];

D(2,:)=[] %usuwanie

-------------------------

function [x] = fibonaczi(n)
if(n==0)
x=0;
endif
if(n==1)
x=1;
endif
if(n>1)
x=fibonaczi(n-1)+fibonaczi(n-2);
endif
endfunction

-------------------------------------

z=magic(4)
[m n] = size(z);
s = 0;
for ii=1:m-1
  for jj=m:-1:ii+1
    s=s+z(ii,jj);
  endfor
endfor
s