
program hesam;Here is a link to the complete (Delphi Pascal) unit that does the FEM Calculations: Below is the $\,40\times 30\,$ grid that has been used for sample calculations and a contour map of some results with $V_0=0$ , $V_1=1$ and random $a_1,a_2$.
Uses Laplace;
procedure test; var k : integer; begin Starten; {Initialize } for k := 0 to 9 do begin { FEM Calculations } Rekenen(Random,Random,0,1); { Store in 'results' file } Opschrijven(k); end; end;
begin test; end.
Sample output ('result' file) - can you see where it is? - :
x y V[x,y] 20 0 5.00000000000000E-0001 20 1 5.00000000000000E-0001 20 2 5.00000000000000E-0001 20 3 4.99999999999999E-0001 20 4 4.99999999999999E-0001 20 5 4.99999999999999E-0001 20 6 4.99999999999999E-0001 20 7 4.99999999999999E-0001 20 8 4.99999999999999E-0001 20 9 5.00000000000000E-0001 20 10 5.00000000000001E-0001 20 11 5.00000000000001E-0001 20 12 5.00000000000001E-0001 20 13 5.00000000000002E-0001 20 14 5.00000000000002E-0001 20 15 5.00000000000002E-0001 20 16 5.00000000000003E-0001 20 17 5.00000000000003E-0001 20 18 5.00000000000003E-0001 20 19 5.00000000000002E-0001 20 20 5.00000000000002E-0001 20 21 5.00000000000003E-0001 20 22 5.00000000000002E-0001 20 23 5.00000000000002E-0001 20 24 5.00000000000001E-0001 20 25 5.00000000000001E-0001 20 26 5.00000000000001E-0001 20 27 5.00000000000001E-0001 20 28 5.00000000000000E-0001 20 29 5.00000000000000E-0001 20 30 5.00000000000000E-0001