Ecuación Diferencial Parcial Elíptica
description
Transcript of Ecuación Diferencial Parcial Elíptica
Ecuacin diferencial parcial elptica
h=1;k=1;>> x=[0:h:4]'
x =
0 1 2 3 4
>> y=[0:k:4]'
y =
0 1 2 3 4
>> iz=exp(y)-cos(y)
iz =
0 2.1780 7.8052 21.0755 55.2518
>> der=exp(y).*cos(4)-exp(4).*cos(y)
der =
-55.2518 -31.2763 17.8910 40.9230 0
>> ab=cos(x)-exp(x)
ab =
0 -2.1780 -7.8052 -21.0755 -55.2518
>> ar=exp(4).*cos(x)-exp(x).*cos(4)
ar =
55.2518 31.2763 -17.8910 -40.9230 0
>> S=solve('4*U1=2.1780+U2-2.1780+U4','4*U2=U1+U3-7.8052+U5','4*U3=U2-31.2763-21.0755+U6','4*U4=7.8052+U5+U1+U7','4*U5=U4+U6+U8+U2','4*U6=U5+17.8910+U9+U3','4*U7=21.0755+U8+31.2763+U4','4*U8=U7+U9-17.8910+U5','4*U9=U8+40.9230-40.9230+U6')
S =U1: [1x1 sym] U2: [1x1 sym] U3: [1x1 sym] U4: [1x1 sym] U5: [1x1 sym] U6: [1x1 sym] U7: [1x1 sym] U8: [1x1 sym] U9: [1x1 sym]
U1=vpa(S.U1,4) U1 = 0.0 >> U2=vpa(S.U2,4) U2 = -5.5106107142857142857142857142857 >> U2=vpa(S.U3,4) U2 = -14.237242857142857142857142857143 >> U3=vpa(S.U3,4) U3 = -14.237242857142857142857142857143 >> U4=vpa(S.U4,4) U4 = 5.5106107142857142857142857142857 >> U5=vpa(S.U5,4) U5 = 0.0 >> U6=vpa(S.U6,4) U6 = 0.91343928571428571428571428571429 >> U7=vpa(S.U7,4) U7 = 14.237242857142857142857142857143 >> U8=vpa(S.U8,4) U8 = -0.91343928571428571428571428571429 >> U9=vpa(S.U9,4) U9 = 0.0