clear all close all cf=500; % W*s/K=J/K kie=5; % W/K; a=-kie/cf; b1=1/cf; b2=kie/cf; c=1; d=0; ex_oven_scheme