Solução sadiku - 15.30 fem

%Questão 30 (FEM)
%SOLUÇÃO PARA EQUAÇÃO DE LAPLACE POR ELEMENTOS FINITOS
%ND= Nº DE NÓS
%NE= Nº DE ELEMENTOS
%NP= Nº DE NÓS FIXOS (ONDE O POTENCIAL É PREESTABELECIDO)
%NDP(I)= Nº DO NÓ COM POTENCIAL PREESTABELECIDO, I=1,2...NP
%VAL(I)= VALOR DO POTENCIAL PREESTABELECIDO NÓ NDP(I)
%NL(I,J)= LISTA DOS NÓS PARA CADA ELEMENTO I, ONDE J=1,2,3 REFERE-SE AO Nº DO NÓ LOCAL
%CE(I,J)= MATRIZ DE COEFICIENTES DO ELEMENTO
%C(I,J)= MATRIZ DE RIGIDEZ GLOBAL
%B(I)= MATRIZ DO LADO DIREITO NO SISTEMA DE EQUAÇÕES SIMULTÂNEAS
%X(I), Y(I)= COORDENADAS GLOBAIS DO NÓ I
%XL(J), YL(J)= COORDENADAS LOCAIS DO NÓ J=1,2,3
% V(I)= POTENCIAL NO NÓ I
% MATRIZES P(I) E Q(I) ESTÃO DEFINIDAS NA EQUAÇÃO 15.62a
%***PRIMEIRA ETAPA***
NE=50;
ND=36;
NP=20;
NL= [1 8 7
1 2 8
2 9 8
2 3 9
3 10 9
3 4 10
4 11 10
4 5 11
5 12 11
5 6 12
7 14 13
7 8 14
8 15 14
8 9 15
9 16 15
9 10 16
10 17 16
10 11 17
11 18 17
11 12 18
13 20 19
13 14 20
14 21 20
14 15 21
15 22 21
15 16 22
16 23 22
16 17 23
17 24 23
17 18 24
19 26 25
19 20 26
20 27 26
20 21 27
21 28 27
21 22 28
22 29 28
22 23 29
23 30 29
23 24 30
25 32 31
25 26 32
26 33 32
26 27 33
27 34 33
27 28 34
28 35 34
28 29 35
29 36 35
29 30 36];
X= [0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0];
Y= [0.0 0.0 0.0 0.0 0.0 0.0 0.2 0.2 0.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.4 0.4 0.6 0.6 0.6 0.6 0.6 0.6 0.8 0.8 0.8 0.8 0.8 0.8 1.0 1.0 1.0 1.0 1.0 1.0];
NDP= [1 2 3 4 5 6 12 18 24 30 36 35 34 33 32 31 25 19 13 7];
VAL= [ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 58.8 95.1 95.1 58.8 0.0 0.0 0.0 0.0 0.0];
%***SEGUNDA ETAPA – CÁLCULO DA MATRIZ DE COEFICIENTES PARA CADA ELEMENTO E //CONEXÃO GLOBAL***
B=zeros(ND,1);
C=zeros(ND,ND); for I=1:NE
%DETERMINAÇÃO DAS COORDENADAS LOCAIS XL(J), YL(J) PARA O ELEMENTO I
K= NL(I, [1:3]);
XL= X(K);
YL= Y(K);
P=zeros(3,1);
Q=zeros(3,1);
P(1)= YL(2) - YL(3);
P(2) = YL(3) -