%heat distribution 2D

clear all;
close all;

cnn_setenv;

mCNN.TemGroup='temlib_diff';

InitState=loadCNN('init_heat.bmp');

State=1.*InitState;
for repeat=1:100
    mCNN.INPUT1=1*zeros(size(State));
    mCNN.STATE=1*State;
    mCNN.Boundary=2;
    mCNN.TimeStep=0.01;
    mCNN.IterNum=3;

    loadtem('HEAT2D')
    runtem;
    
    State=1.*mCNN.STATE;

    subplot(1,2,1);
    surf(InitState);
    subplot(1,2,2);
    surf(mCNN.OUTPUT);
    
    drawnow
end