Rozmowy o pięknym sporcie, jakim jest z pewnością piłka nożna
1. Napisać funkcję generującą wartości losowe z rozkładu gamma τkxk−1 exp(−xτ)Γ(k), dla parametrów k = 9,τ = 0.5 (2 pkt.)
funkcja.m
function [ gamma ] = funkcja( x,t,k )
gamma=(t.^k).*(x.^(k-1)).*(exp(-x*t))
end
while liczba <= 1000
x=rand()*zakres;
y=rand();
if y < tau^k*x^(k-1)*exp(-x*tau)/gamma(k)
wyniki(liczba)=x;
liczba=liczba+1;
end
end
2. Narysować krzywą gęstości (1 pkt.)
krzywa.m
x=(1:60);
t=0.5;
k=9;
funkcja(x,t,k);
plot(x,ans);
3. Stworzyć histogram z wylosowanych 1000 próbek (2 pkt.)
* Stworzyć animację z losowań do histogramu (2 pkt.)
histogram.m
clear
for i=1:1000
x=rand(1,1000);
t=0.5;
k=9;
funkcja(x,t,k);
hist(ans);
drawnow
pause(0.5);
end
http://www.mathworks.com/help/stats/makedist.html
::::::::::::::::::::::::::::::::::::::::::::::GLINAFUN::::::::::::::::::::::::::::::::::::::::::::::::::
Offline
circle.m
function o = circle(x,y,r)
hold on
th = 0i/50:2*pi;
xunit = r * cos(th) + x;
yunit = r * sin(th) + y;
o = plot(xunit, yunit);
funkcja1.m
function [x y] = funkcja1(r,R,h,t)
x = (R-r)*cos(t)+h*cos((R-r)*t/r);
y = (R-r)*sin(t)-h*sin((R-r)*t/r);
end
funkcja2.m
function [x y] = funkcja2(r,R,h,t)
x = (R+r)*cos(t)-h*cos((R+r)*t/r);
y = (R+r)*sin(t)-h*sin((R+r)*t/r);
end
krzywa.m
clear
R = 5;
r = 3;
h = 5;
t = (1:0.1:500);
[x1 y1] = funkcja2(r,R,h,t);
plot(x1,y1)
circle(0,0,R)
circle(R-r,0,r)
hold on
x=(R-r:0.01:R-r+h);
plot(x, 0);
zad3.m
clear
prompt = 'podaj R';
R = input(prompt);
prompt = 'podaj r';
r = input(prompt);
prompt = 'podaj h';
h = input(prompt);
t = (1:0.1:500);
[x1, y1] = funkcja1(r,R,h,t);
plot(x1,y1)
circle(0,0,R)
circle(R-r,0,r)
hold on
x=(R-r:0.01:R-r+h);
plot(x, 0);
zad4.m
clear;
r=2;
h=3;
R=10;
d=linspace(0,2*pi,500);
%duze kolo
XR=R*cos(d);
YR=R*sin(d);
plot(XR,YR,'red');
f=0;
for t = linspace(0,2*pi,500);
f=f+1;
hold on;
%ta epichondria czy cos
[x(f) y(f)]= funkcja1(r,R,h,t);
%male kolo
xr=r*cos(d)+(R-r)*cos(t);
yr=r*sin(d)+(R-r)*sin(t);
%promien rysowania z malego kola
prx=(R-r)*cos(t);
pry=(R-r)*sin(t);
X=linspace(x(f),prx,1000);
Y=linspace(y(f),pry,1000);
%rysowanie
h3=plot(x,y);
h1=plot(xr,yr,'black');
h2=plot(X,Y);
%pauza i usuwanie
pause(0.0002);
delete(h1);
delete(h2);
delete(h3);
end
%na koniec wyswietlenie wszystkiego
h1=plot(xr,yr,'black');
h2=plot(X,Y);
h3=plot(x,y);
::::::::::::::::::::::::::::::::::::::::::::::GLINAFUN::::::::::::::::::::::::::::::::::::::::::::::::::
Offline
tygrysy.m
clear
pierwszy=1;
x=rand(1,30)*100;
y=rand(1,30)*100;
plot(x,y,'o')
hold on
for i=1:1:30
if y(pierwszy)>y(i)
pierwszy=i;
end
end;
buffor = [x(1) y(1)];
x(1)=x(pierwszy);
y(1)=y(pierwszy);
x(pierwszy)=buffor(1);
y(pierwszy)=buffor(2);
for i=1:1:30
d(i)=abs(x(i)-x(1))+abs(y(i)-y(1))
end;
for i=2:1:30
if x(i)<x(1)
alfa(i)=2-(y(i)-y(1))/d(i);
else
alfa(i)=(y(i)-y(1))/d(i);
end;
end;
n=2;
numerlisty=[1 3 3 3 2 2 2 2 2 2 2 2 2 2 8 9 9 2 2 2 2 2 2 2 2 2 2 2 2 2 ];
while n<31
for i=2:1:30
if alfa(i)< alfa(numerlisty(n))
numerlisty(n)=i;
end
end
alfa(numerlisty(n))=1000;
n=n+1;
end
nowyx=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
nowyy=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
%%%
iks=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
ygrek=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
for i=1:1:30
iks(i)=x(i);
ygrek(i)=0;
plot(iks,ygrek,'g')
end
%%
for i=1:1:30
nowyx(i)=x(numerlisty(i));
nowyy(i)=y(numerlisty(i));
end
plot(nowyx,nowyy,'g')
x=nowyx;
y=nowyy;
stosx =[x(3) x(2) x(1)];
stosy = [y(3) y(2) y(1)];
for i=4:1:30
while (wyznacznik(stosx(2),stosy(2),stosx(1),stosy(1),x(i),y(i)))>0 %&&length(stosx)>2
stosx=pop(stosx);
stosy=pop(stosy);
end
stosx=push(stosx,x(i));
stosy=push(stosy,y(i));
end
stosx=push(stosx,stosx(end));
stosy=push(stosy,stosy(end));
plot(stosx, stosy, 'r')
pop.m
function o = pop(stos)
stos(1) = [];
o =stos;
push.m
function o = push(stos, punkt)
dlugosc = length (stos);
n=dlugosc+1;
while n>1
stos(n)=stos(n-1);
n=n-1;
end
stos(1)=punkt;
o =stos;
wyznacznik.m
function o = wyznacznik(p1x, p1y, p2x, p2y, p3x, p3y)
o =p1y*p2x+p2y*p3x+ p3y*p1x- p3y* p2x- p1y*p3x- p2y*p1x;
::::::::::::::::::::::::::::::::::::::::::::::GLINAFUN::::::::::::::::::::::::::::::::::::::::::::::::::
Offline
::::::::::::::::::::::::::::::::::::::::::::::GLINAFUN::::::::::::::::::::::::::::::::::::::::::::::::::
Offline
clear;
i=1; %x, y, g, p, d, l
mapa=zeros(20,20,5);
odwiedzone=zeros(20,20);
polozeniex=1;
polozeniey=1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(20,20,5)=1;
ostatniruch=[];
while any(1==any(0==odwiedzone))
a=rand*4;
if a==4
ruch=0
else
ruch=floor(a) %wartosc calkowita
end
if ruch==0 && polozeniey>1 && odwiedzone(polozeniex,polozeniey-1)==0 %gora
mapa(polozeniex,polozeniey,1)=1;
polozeniey=polozeniey-1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,3)=1;
ostatniruch=push(ostatniruch,1);
elseif ruch==1 && polozeniex<20 && odwiedzone(polozeniex+1,polozeniey)==0 %prawo
mapa(polozeniex,polozeniey,2)=1;
polozeniex=polozeniex+1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,4)=1;
ostatniruch=push(ostatniruch,2);
elseif ruch==2 && polozeniey<20 && odwiedzone(polozeniex,polozeniey+1)==0 %dol
mapa(polozeniex,polozeniey,3)=1;
polozeniey=polozeniey+1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,1)=1;
ostatniruch=push(ostatniruch,3);
elseif ruch==3 && polozeniex>1 && odwiedzone(polozeniex-1,polozeniey)==0 %lewo
mapa(polozeniex,polozeniey,4)=1;
polozeniex=polozeniex-1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,2)=1;
ostatniruch=push(ostatniruch,4);
%%%%%%%%%%%%%%%%%%%%%%%%%
elseif polozeniey>1 && odwiedzone(polozeniex,polozeniey-1)==0
mapa(polozeniex,polozeniey,1)=1;
polozeniey=polozeniey-1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,3)=1;
ostatniruch=push(ostatniruch,1)
elseif polozeniex<20 && odwiedzone(polozeniex+1,polozeniey)==0
mapa(polozeniex,polozeniey,2)=1;
polozeniex=polozeniex+1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,4)=1;
ostatniruch=push(ostatniruch,2);
elseif polozeniey<20 && odwiedzone(polozeniex,polozeniey+1)==0
mapa(polozeniex,polozeniey,3)=1;
polozeniey=polozeniey+1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,1)=1;
ostatniruch=push(ostatniruch,3);
elseif polozeniex>1 && odwiedzone(polozeniex-1,polozeniey)==0
mapa(polozeniex,polozeniey,4)=1;
polozeniex=polozeniex-1;
odwiedzone(polozeniex,polozeniey)=1;
mapa(polozeniex,polozeniey,5)=1;
mapa(polozeniex,polozeniey,2)=1;
ostatniruch=push(ostatniruch,4);
elseif ostatniruch(1)==1
ostatniruch=pop(ostatniruch);
polozeniey=polozeniey+1;
elseif ostatniruch(1)==2
ostatniruch=pop(ostatniruch);
polozeniex=polozeniex-1;
elseif ostatniruch(1)==3
ostatniruch=pop(ostatniruch);
polozeniey=polozeniey-1;
elseif ostatniruch(1)==4
ostatniruch=pop(ostatniruch);
polozeniex=polozeniex+1;
end
end
rectangle('Position',[0,0,25,25],'FaceColor',[0 0 0])
hold on
x0=0;
y0=20;
for j=1:1:20
for i=1:1:20
if mapa(i,j,3)==1
rectangle('Position',[i-0.45,j-0.45,0.9,1],'FaceColor',[1 1 1],'EdgeColor','w')
end
if mapa(i,j,2)==1
rectangle('Position',[i-0.45,j-0.45,1,0.9],'FaceColor',[1 1 1],'EdgeColor','w')
end
if mapa(i,j,1)==1
rectangle('Position',[i-0.45,j-0.55,0.9,1],'FaceColor',[1 1 1],'EdgeColor','w')
end
if mapa(i,j,4)==1
rectangle('Position',[i-0.55,j-0.45,1,0.9],'FaceColor',[1 1 1],'EdgeColor','w')
end
end
end
function o = push(stos, punkt)
dlugosc = length (stos);
n=dlugosc+1;
while n>1
stos(n)=stos(n-1);
n=n-1;
end
stos(1)=punkt;
o =stos;
function o = pop(stos)
stos(1) = [];
o =stos;
::::::::::::::::::::::::::::::::::::::::::::::GLINAFUN::::::::::::::::::::::::::::::::::::::::::::::::::
Offline
clear;
i=1;
mapa=zeros(20,20,4);
odwiedzone=zeros(20,20);
wspolrzednaex=2;
wspolrzednay=20;
odwiedzone(wspolrzednay,wspolrzednaex)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
drogax=[wspolrzednaex];
drogay=[wspolrzednay];
ostatniruch=[];
pktkoncowy=[20 1];
znaleziony=0;
while any(1==any(0==odwiedzone))
a=rand*3;
if a==3
ruch=0
else
ruch=floor(a) %wartosc calkowita
end
if ruch==0 && mod(wspolrzednaex,2)==0 && mod(wspolrzednay,2)~=0 && wspolrzednay<20 && odwiedzone(wspolrzednaex,wspolrzednay+1)==0 %dol
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay+1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif ruch==0 && mod(wspolrzednaex,2)==0 && mod(wspolrzednay,2)==0 && wspolrzednay>1 && odwiedzone(wspolrzednaex,wspolrzednay-1)==0 %gora
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay-1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif ruch==0 && mod(wspolrzednaex,2)~=0 && mod(wspolrzednay,2)~=0 && wspolrzednay>1 && odwiedzone(wspolrzednaex,wspolrzednay-1)==0 %gora
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay-1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif ruch==0 && mod(wspolrzednaex,2)~=0 && mod(wspolrzednay,2)==0 && wspolrzednay<20 && odwiedzone(wspolrzednaex,wspolrzednay+1)==0 %dol
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay+1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif ruch==1 && wspolrzednaex<20 && odwiedzone(wspolrzednaex+1,wspolrzednay)==0 %prawo
mapa(wspolrzednaex,wspolrzednay,2)=1;
wspolrzednaex=wspolrzednaex+1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,1)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
ostatniruch=push(ostatniruch,2);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif ruch==2 && wspolrzednaex>1 && odwiedzone(wspolrzednaex-1,wspolrzednay)==0 %lewo
mapa(wspolrzednaex,wspolrzednay,1)=1;
wspolrzednaex=wspolrzednaex-1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,2)=1;
ostatniruch=push(ostatniruch,1);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
%%%%%%%%%%%%%%%%%%%%%%%%%
elseif mod(wspolrzednaex,2)==0 && mod(wspolrzednay,2)~=0 && wspolrzednay<20 && odwiedzone(wspolrzednaex,wspolrzednay+1)==0 %dol
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay+1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif mod(wspolrzednaex,2)==0 && mod(wspolrzednay,2)==0 && wspolrzednay>1 && odwiedzone(wspolrzednaex,wspolrzednay-1)==0 %gora
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay-1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif mod(wspolrzednaex,2)~=0 && mod(wspolrzednay,2)~=0 && wspolrzednay>1 && odwiedzone(wspolrzednaex,wspolrzednay-1)==0 %gora
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay-1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif mod(wspolrzednaex,2)~=0 && mod(wspolrzednay,2)==0 && wspolrzednay<20 && odwiedzone(wspolrzednaex,wspolrzednay+1)==0 %dol
mapa(wspolrzednaex,wspolrzednay,3)=1;
wspolrzednay=wspolrzednay+1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,3)=1;
ostatniruch=push(ostatniruch,3);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif wspolrzednaex<20 && odwiedzone(wspolrzednaex+1,wspolrzednay)==0 %prawo
mapa(wspolrzednaex,wspolrzednay,2)=1;
wspolrzednaex=wspolrzednaex+1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,1)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
ostatniruch=push(ostatniruch,2);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif wspolrzednaex>1 && odwiedzone(wspolrzednaex-1,wspolrzednay)==0 %lewo
mapa(wspolrzednaex,wspolrzednay,1)=1;
wspolrzednaex=wspolrzednaex-1;
odwiedzone(wspolrzednaex,wspolrzednay)=1;
mapa(wspolrzednaex,wspolrzednay,4)=1;
mapa(wspolrzednaex,wspolrzednay,2)=1;
ostatniruch=push(ostatniruch,1);
if znaleziony==0
drogax=push(drogax,wspolrzednaex);
drogay=push(drogay,wspolrzednay);
end
elseif ostatniruch(1)==1
ostatniruch=pop(ostatniruch);
wspolrzednaex=wspolrzednaex+1;
if znaleziony==0
drogax=pop(drogax);
drogay=pop(drogay);
end
elseif ostatniruch(1)==2
ostatniruch=pop(ostatniruch);
wspolrzednaex=wspolrzednaex-1;
if znaleziony==0
drogax=pop(drogax);
drogay=pop(drogay);
end
elseif ostatniruch(1)==3 && mod(wspolrzednaex,2)==0 && mod(wspolrzednay,2)~=0
ostatniruch=pop(ostatniruch);
wspolrzednay=wspolrzednay+1;
if znaleziony==0
drogax=pop(drogax);
drogay=pop(drogay);
end
elseif ostatniruch(1)==3 && mod(wspolrzednaex,2)==0 && mod(wspolrzednay,2)==0
ostatniruch=pop(ostatniruch);
wspolrzednay=wspolrzednay-1;
if znaleziony==0
drogax=pop(drogax);
drogay=pop(drogay);
end
elseif ostatniruch(1)==3 && mod(wspolrzednaex,2)~=0 && mod(wspolrzednay,2)~=0
ostatniruch=pop(ostatniruch);
wspolrzednay=wspolrzednay-1;
if znaleziony==0
drogax=pop(drogax);
drogay=pop(drogay);
end
elseif ostatniruch(1)==3 && mod(wspolrzednaex,2)~=0 && mod(wspolrzednay,2)==0
ostatniruch=pop(ostatniruch);
wspolrzednay=wspolrzednay+1;
if znaleziony==0
drogax=pop(drogax);
drogay=pop(drogay);
end
end
if wspolrzednay==pktkoncowy(2) && wspolrzednaex==pktkoncowy(1)
znaleziony=1;
end
end
rectangle('Position',[0,0,25,25],'FaceColor',[1 1 1])
hold on
for j=1:1:20
for i=1:1:20
if mod(i+j,2)==0
x=i;
y=j;
if mapa(i,j,1)==0
x2 = [x x+1 ];
y2 = [y y+1 ];
plot(x2,y2,'green')
end
if mapa(i,j,2)==0
x2 = [x+1 x+2];
y2 = [y+1 y];
plot(x2,y2,'green')
end
if mapa(i,j,3)==0
x2 = [x x+2 ];
y2 = [y y ];
patch(x2,y2,'w','EdgeColor','green')
end
end
if mod(i+j,2)~=0
x=i;
y=j;
if mapa(i,j,1)==0
x2 = [x x+1 ];
y2 = [y+1 y ];
plot(x2,y2,'green')
end
if mapa(i,j,2)==0
x2 = [x+1 x+2];
y2 = [y y+1];
plot(x2,y2,'green')
end
if mapa(i,j,3)==0
x2 = [x x+2 ];
y2 = [y+1 y+1 ];
patch(x2,y2,'w','EdgeColor','green')
end
end
end
end
drogax=drogax+1;
drogay=drogay+0.5;
plot(drogax, drogay,'r')
clc;
clear all;
close all;
load data_map.mat;
x = data_map(:,1,1);
y = data_map(:,2,1);
z = data_map(:,3,1);
for i=1:length(x)
x1(i) = round(x(i)*1000)/1000;
y1(i) = round(y(i)*1000)/1000;
z1(i) = round(z(i)*1000)/1000;
end
[X,Y] = meshgrid(58.807:0.001:59.227, 55.28:0.001:55.7);
for i=1:421
for j=1:421
Z(i,j) = 0;
end
end
for n=1:2327
k(n) = round((x1(n) - 58.811)/0.001)+1;
l(n) = round((y1(n) - 55.28)/0.001)+1;
Z(l(n),k(n)) = z(n);
end
Z1=Z;
for a=1:100
for i=2:420
for j=2:420
if Z(i,j)==0
Z1(i,j)= ( Z1(i-1,j) + Z1(i+1,j) + Z1(i,j-1) + Z1(i,j+1) )/4;
end
end
end
end
Z3=zeros(421,421);
for i=1:420
for j=1:420
Z3(i,j)=Z1(i,j);
end
mesh(X,Y,Z3);
pause (0.002);
end
mesh(X,Y,Z1)
shading interp
%Z=interp2(x1,y1,z1,X,Y);
%figure(2)
%mesh(X,Y,Z)
%surf(Z,'EdgeColor','None');
%shading interp
%figure(3)
%surf(X,Y,Z)
%shading interp
%Z=interp2(x1,y1,z1,X,Y);
%figure(2)
%mesh(X,Y,Z)
%surf(Z,'EdgeColor','None');
%shading interp
%figure(3)
%surf(X,Y,Z)
%shading interp
::::::::::::::::::::::::::::::::::::::::::::::GLINAFUN::::::::::::::::::::::::::::::::::::::::::::::::::
Offline