Cómo resolver un problema de valor límite utilizando el método de diferencia finita

Estoy resolviendo un problema dado para h = 0.25. La solución se adjunta en imágenes.

¡¡¡DISFRUTAR!!!

1

2

3

CODIGO MATLAB

a = [- 4 2,25 0; 1.75 -4 2.25; 0 1.75 -4]

a =

-4.0000 2.2500 0

1.7500 -4.0000 2.2500

0 1.7500 -4.0000

>> b = -. 125 * [1 1 1] ‘

b =

-0,1250

-0,1250

-0,1250

>> y = a \ b

y =

0.1005

0,1231

0,0851

>> x = [0 .25 .5 .75 1]

x =

0 0.2500 0.5000 0.7500 1.0000

>> y = [0 0.1005 0.1231 0.0851 0]

y =

0 0.1005 0.1231 0.0851 0

>> trama (x, y)

>> rejilla en

Resultados