mn lab 7

function matrix = galerkin(panels)

 [rows,cols,numpanels] = size(panels);
 matrix = zeros(numpanels, numpanels);

 M = 9;
 [tempx tempy] = meshgrid(1/(2*M):1/M:1,1/(2*M):1/M:1);
 evalgrid = [reshape(tempx,M*M,1) reshape(tempy,M*M,1) zeros(M*M,1)];

for i=1:numpanels
    numverts = panels(1,1,i);
    panel_i = panels(2:numverts+1,:,i);

    for i = 1 : 3
        if (i == 3)
            evalpoints(:, i) = 0;
        else
            evalpoints(:, i) = evalgrid(:, i)*(1/sqrt(numpanels))+ panel_i(1, i);
        end
    end
    for j=1:numpanels
        numverts = panels(1,1,j);
        panel = panels(2:numverts+1,:,j);
        [area, collocpt, Z, fss] = calcp(panel, evalpoints);
        matrix(i,j) = sum(fss*(area/(M*M)));
    end
end