Hola Bruna! Está perfecto el desarrollo, se resuelve así.
- Tenes 5 nodos --> 5 incognitas.
- Las condiciones de borde te dan información sobre y0 y y4 --> 2 datos.
- 5 incognitas - 2 datos = 3 ecuaciones --> te queda un sistema de 3x3.
Lo único, resolviste mal el sistema porque te comiste el + y3 en la ec de la 2da fila. Debería dar y1=0, y2=0, y3= 0.25
En wolfram podes chequear la solucion del sistema: https://www.wolframalpha.com/input?i=%7B%7B-2%2C1%2C0%7D%2C%7B1%2C-2%2C1%7D%2C%7B0%2C1%2C-4%2F3%7D%7D*%7B%7By1%7D%2C%7By2%7D%2C%7By3%7D%7D+%3D+%7B%7B0*0.5%5E2%7D%2C%7B1*0.5%5E2%7D%2C%7B0*0.5%5E2-1%2F3%7D%7D
(para y1 e y2 da un numero muy chico en lugar de 0 por error de redondeo, si lo haces a mano da 0, pero para darse una idea está bueno)
Saludos!