#include "PDEHeat.h"

/*
 * Constructor for the heat equation PDE solver. Parameters are the same names
 * as are documented in PDEHeat.h.
 * 
 * PDE's constructor is called implicitly in the initialisation list. Remember
 * to set initial and boundary conditions here.
 */
PDEHeat::PDEHeat(int N, double dt, double a, double b, double theta, double D) :
  PDE(N, dt, a, b, theta)
{
  // You need to fill in this function.
  
}

/*
 * Sets the operators this->le (differential operator at time n) and this->li
 * (differential operator at time n+1). Remember to leave the first and last
 * rows of these matrices set to zero.
 */
void PDEHeat::setOperator() {
  // You need to fill in this function.
  
}

/*
 * Sets the boundary condition for this PDE. As a test, you might want to set
 * these to zero and use the initial condition defined in the comment for
 * setIC().
 */
void PDEHeat::setBC() {
  // You need to fill in this function.
  
}

/*
 * Sets the initial condition for this PDE. As a test you might want to use a
 * function like (x-a)*(b-x)
 */
void PDEHeat::setIC() {
  // You need to fill in this function.
  
}