Interpolacion De Una Matriz

CODIGO DEL PROGRAMA

program Interpolacion_Lagrange_2D
integer n,m,i,j
dimension f(0:30,0:10)
dimension x(0:10)
dimension y(0:10)
data n,m/4,5/
data(x(i),i=0,4)/0.0,0.5,1.0,1.5,2.0/
data (y(j),j=0,5)/0.0,0.5,1.0,1.5,2.0,2.5/
data((f(i,j),j=0,5),i=0,4)/0.0,5.0,10.00,15.00,20.00,25.00,5.00,7.51,10.05,12.70,15.67,20.00,10.00,10.00,10.00,10.00,10.00,10.00,15.00,12.51,9.95,7.32,4.33,0.0,20.00,15.00,10.00,5.00,0.00,-5.00/

print *
print *, 'Interpolacion de Lagrange en 2D'print *
print *, ' Tabla de valores '
print *
print *, 0.00 0.5 1.0 1.5 2.0 2.5 '
print *,' x\y________________________________________________'
print *,' '
print *, ' 0.0 0.00 5.00 10.00 15.00 20.00 25.00 '
print *, ' 0.5 5.00 7.51 10.05 12.70 15.67 20.00 '
print *, ' 1.0 10.00 10.0010.00 10.00 10.00 10.00 '
print *, ' 1.5 15.00 12.51 9.95 7.32 4.33 0.00 '
print *, ' 2.0 20.00 15.00 10.00 5.00 0.00 -5.00 '
print *

10 print *,'Ingrese x:'read *, val1
print *,'Ingrese y:'
read *, val2

if (val1.lt.x(0) .or. val1.gt.x(4)) then
print *, 'Cuidado!!: X esta en el rango de extrapolacion'print *, 'Desea continuar? 0-continuar/1-Salir'
print*
read *, resp
print *
if(resp.eq.0) go to 10
print *
stop 'fin del programa'
end if



if (val2.lt.y(0) .or. val2.gt.y(5)) then
print *, 'Cuidado!!: y esta en elrango de extrapolacion'
print *
print *, 'Desea continuar? 0-continuar/1-Salir'
print*
read *, resp
print *...