%Transient conduction

clear all; clc;
C1 = 1.1191;
zeta1 = 0.8603;
Fo = 0:0.1:3;
theta_star_o = C1*exp(-zeta1^2*Fo);

%Plot 
figure(1); %Create figure #1
clf; %Clear current figure
h=plot(Fo,theta_star_o,'-k'); 
set(gca,'Box','on','LineWidth',2,...
'FontName','Helvetica',...
'FontSize',14);
%Set axis properties
set(h,'LineWidth',2);
%Set linewidth for curve
xlabel('t');
ylabel('\theta(0,t)');
%title('Transient conduction: Plane wall with convection');
axis square; %Make axis box square

%Read and plot ANSYS results
ANSYS = dlmread('thetao_vs_t_dt_01.txt');
hold on;
h2=plot(ANSYS(:,1),ANSYS(:,2),'--r'); 
set(h2,'LineWidth',2);
%legend('One term series solution','Numerical solution (ANSYS)');

ANSYS2 = dlmread('thetao_vs_t.txt');
hold on;
h2=plot(ANSYS2(:,1),ANSYS2(:,2),':b'); 
set(h2,'LineWidth',2);
legend('One term series solution','\Deltat = 0.01',...
    'ANSYS, Default \Deltat');