ecuaciones diferenciales parciales
DESCRIPTION
Programación en Octave de problema de la Ecuación de CalorTRANSCRIPT
-
Alumno: Snchez Prez Erik Hernn Materia: Ecuaciones Diferenciales Parciales II
Proyecto Parte 2 Resultado de coordenadas Cdigo: fprintf('Resultado de coordenadas\n');
fprintf('----------------------------------------\n');
fprintf(' RADIO x y \n');
fprintf('----------------------------------------\n');
for i=20:2:26
for j=0:9
x=i*cos(j*(pi/9));
y=i*sin(j*(pi/9));
fprintf('%6.0f %15.6f %15.6f \n',i,x,y);
end
end
for i=22:2:26
x=i*cos(pi/2);
y=i*sin(pi/2);
fprintf('%6.0f %15.6f %15.6f \n',i,x,y);
end
Matriz de conectividad
Cdigo: n=52;
nd=3;
da=9;
dr=3;
c=zeros(n,nd);
k=0;
j=1;
for i=1:((da-1)*(dr))/2
if j==4
k=k+1;
j=1;
endif
-
Alumno: Snchez Prez Erik Hernn Materia: Ecuaciones Diferenciales Parciales II
c(((2*i)-1),1)=i+k;
c(((2*i)-1),2)=i+1+k;
c(((2*i)-1),3)=i+4+k;
j=j+1;
endfor
k=0;
j=1;
for i=1: ((da-1)*(dr))/2
if j==4
k=k+1;
j=1;
endif
c(((2*i)),1)=i+1+k;
c(((2*i)),2)=i+4+k;
c(((2*i)),3)=i+5+k;
j=j+1;
endfor
k=0;
j=1;
for i=1: 3
if j==2
k=k+1;
j=1;
endif
c((((dr)*((da-1))))+((2*i)-1)+k,1)=(((dr)*((da-1)))-7)+i-1;
c((((dr)*((da-1))))+((2*i)-1)+k,2)=(((dr)*((da-1)))-7)+i;
c((((dr)*((da-1))))+((2*i)-1)+k,3)=(((dr)*((da-1)))-7)+i+3;
j=j+1;
endfor
k=0;
j=1;
for i=1: 3
-
Alumno: Snchez Prez Erik Hernn Materia: Ecuaciones Diferenciales Parciales II
if j==2
k=k+1;
j=1;
endif
c((((dr)*((da-1))))+((2*i)+1)+k,1)=(((dr)*((da-1)))-7)+i+3;
c((((dr)*((da-1))))+((2*i)+1)+k,2)=(((dr)*((da-1)))-7)+i+6;
c((((dr)*((da-1))))+((2*i)+1)+k,3)=(((dr)*((da-1)))-7)+i+7;
j=j+1;
endfor
k=0;
j=1;
for i=1: 3
if j==2
k=k+1;
j=1;
endif
c((((dr)*((da-1))))+((2*i))+k,1)=(((dr)*((da-1)))-7)+i;
c((((dr)*((da-1))))+((2*i))+k,2)=(((dr)*((da-1)))-7)+i+7;
c((((dr)*((da-1))))+((2*i))+k,3)=(((dr)*((da-1)))-7)+i+3;
j=j+1;
endfor
k=0;
j=1;
for i=1: ((da-1)*(dr))/2
if j==4
k=k+1;
j=1;
endif
c((((dr)*((da-1)))+9)+((2*i)),1)=(((dr)*((da-1))))+i+k-1;
c((((dr)*((da-1)))+9)+((2*i)),2)=(((dr)*((da-1))))+i+k;
c((((dr)*((da-1)))+9)+((2*i)),3)=(((dr)*((da-1))))+i+4+k;
j=j+1;
-
Alumno: Snchez Prez Erik Hernn Materia: Ecuaciones Diferenciales Parciales II
endfor
k=0;
j=1;
for i=1: ((da-1)*(dr))/2
if j==4
k=k+1;
j=1;
endif
c((((dr)*((da-1)))+9)+((2*i-1)),1)=(((dr)*((da-1))))+i-1+k;
c((((dr)*((da-1)))+9)+((2*i-1)),2)=(((dr)*((da-1))))+i+3+k;
c((((dr)*((da-1)))+9)+((2*i-1)),3)=(((dr)*((da-1))))+i+4+k;
j=j+1;
endfor
c