%Script for simulation of isomerpar.m
    clear
    clc
    close all
    
%%Load Data

    load  Data_ex5_E2012.mat


%% Define model parameters
    umaxl = 0.0887;
    ksl = 0.0194 ;
    kil = 6.502 ;
    umaxz = 0.2016;
    ksz = 4.993;
    kd = 0.0094;
    yl = 0.4106;
    yz= 0.6418;
    alpha=0.0525;
    betha=0.0026;

    par = [umaxl ksl kil umaxz ksz kd yl yz alpha betha];
%% Define initial conditions
    lac_cont = 55; 
    xyl_cont = 42.5;
    bio_cont = 1;
    pro_cont = 0.5;

    x0=[lac_cont xyl_cont bio_cont pro_cont];
%% Define simulation start/end time
    t0=0;
    tfin=100;

%% Simulation options were set in our_simulink.mdl;
    %Solve the model
    sim('isomerpar_simulink_mSfunc')

%% Plot the results
    figure
subplot (2,2,1)
hold on
plot(t,y(:,1),'r','linewidth',1.5)
plot(t1,L_data,'.','linewidth',.5)
xlabel('')
ylabel('Lactose [g/L]')
axis([0 100 0 60])
grid on
hold off

subplot (2,2,2)
hold on
plot(t,y(:,2),'r','linewidth',1.5)
plot(t1,Z_data,'.','linewidth',.5)
hold off
xlabel('')
ylabel('Xylose [g/L]')
axis([0 100 0 50])
grid on

subplot (2,2,3)
hold on
plot(t,y(:,3),'r','linewidth',1.5)
hold off
xlabel('Time (h)')
ylabel('Biomasse [g/L]')
axis([0 100 0 50])
grid on

subplot (2,2,4)
hold on
plot(t,y(:,4),'r','linewidth',1.5)
plot(t1,P_data,'.','linewidth',.5)
hold off
xlabel('Time')
ylabel('Protein [g/L]')
axis([0 100 0 8])
grid on

%% Save

saveas(gcf,'E5','jpg')