Subroutine Fcn(n, t, Y, Yprime)

Implicit none

Integer i, j, n

Real*8 Cp, k, L, dx, Rho, t, Y1, Y(n), Yprime(n)

Y1 = 50

k = 37.7d0

Rho = 7822.d0

Cp = 444.d0

L = 0.10d0

! n is the number of nodes in the slab (x).

dx = (L - 0.d0)/dble(n - 1)

! dT/dt = Yprime = Derivative of the temperature.

! If the program starts with Debug, this condition is necessary.

Yprime(1) = 0.0

! If the program starts with Debug, this condition is not necessary.

Yprime(n) = 0.0

Do i = 2, n-1

! Symmetry condition at x = 0 m.

Y(1) = Y(2)

! Evaluate the derivative of the temperature,

!.

! is the temperature of the surface at.

End Do

End subroutine