%% ZADANIE 7)
% Wczytaj dane zapisane w pliku V0.mat i wyœwietl na rysunku 
% siatkê elementów skoñczonych oraz wartoœci potencja³u V0 
% dla punktów okreœlonych wspó³rzêdnymi zapisanymi w danych x i y

clear;
close all;
clc

load V0;		  % plik pobierz ze strony: ks.zut.edu.pl/mm/V0.mat

t(4,:) = t(1,:);
figure(1), hold on,
for w = 1:size(t,2),
  plot(p(1,t(1:4,w)),p(2,t(1:4,w))),
  text(mean(p(1,t(1:3,w))),mean(p(2,t(1:3,w))),num2str(w)),
end
axis equal tight;

figure(2),
surf(x,y,reshape(V0,5,5));
axis equal tight;
colorbar

%% ZADANIE 8)
% Wykorzystaj funkcje mesSPR i mesV do wykonania
% aproksymacji dla innych wspó³rzêdnych

x=0:3:30;  y=0:3:30;  VV(1:length(x),1:length(y))=NaN;
for k = 1:length(y),
  for w = 1:length(x),
    for e = 1:size(t,2),
      if mesSPR(p(1,t(1:3,e)),p(2,t(1:3,e)),x(w),y(k))
        VV(w,k)=mesV(p(1,t(1:3,e)),p(2,t(1:3,e)),x(w),y(k),V0(t(1:3,e)));
      end
    end
  end
end
figure(3), surf(x,y,VV);
axis equal tight;
colorbar
