Programa "Elementos finitos"

Funciona para resolver 

(laplaciano)U=F(x,y)
U=0 en la frontera

es con elementos triangulares; se meten los puntos (Primero los puntos
interiores en los que las funciones de forma Wi=1 en el punto i;
la funcin W1 vale 1 en el punto 1; la funcin W2 vale 1 en el punto 2, etc.)

Primero pide el nmero de puntos (total) y sus coordenadas.

Luego pide el nmero de elementos finitos triangulares y los nmeros de los
puntos que son los vrtices de cada tringulo.

A continuacin, pide el nmero de funciones de forma (W) que se van a utilizar
Los puntos que definen las W son los primeros, as que no es necesario dar
ms datos.

Luego el programa da el valor de las derivadas segn x y segn y de cada
W en cada conjunto K, seguido por las superficies de cada conjunto K.

Luego da los elementos de la matriz, gradientes de los productos de las 
funciones de forma, es decir, $W1W2 significa "Integral de W1 por W2"
Dado que la matriz es simtrica ($WiWj=$WjWi), solo se dan los elementos del 
tringulo superior de la matriz.

Despus se pide la funcin F, que puede ser:
-F constante en cada elemento K, pero distinto entre los diferentes K
 Entonces se pide el valor de F en cada K.
-F constante en todo el recinto.
 Simplemente se pide el valor de F.
-F funcin cualquiera
 Para usar esta opcin hay que haber definido previamente en las lneas
 5 a 18 la funcin en la forma 5 F=F(x,y), por ejemplo, si nos dicen que 
 F vale X^2-3*X*Y^4, escribiramos:
 5 F=X^2-3*X*Y^4
 Si la funcin vara de conjunto a conjunto, por ejemplo F=XY en K1, F=3*X^2 en
 K2, F=XY^3 en K5, F=0 en el resto, haramos
 5 IF R+1=1 THEN F=X*Y:RETURN
 6 IF R+1=2 THEN F=3*X^2:RETURN
 7 IF R+1=5 THEN F=X*Y^:RETURN
 8 F=0:RETURN

Despus, el programa calcula los coeficientes de cada W y los da; eso es la
solucin aproximada del problema.

A continuacin, puede calcularse el valor de la funcin en cualquier punto
indicando en que conjunto K est y diciendo sus coordenadas.

(El programa necesita tener la subrutina ECUACION en el rea P2!)

Longitud del programa en la calculadora: 2993 Bytes.
Longitud de la subrutina ECUACION: 497 Bytes.