%script file for solving induced.m as a differential equation
clear all;
global beta km n d2 inducer %declare parameters as global

beta = input('Maximum expression level given excess inducer '); %requests input for beta
km = input('Concentration of inducer required for half-maximal expression '); %requests input for km
n = input('Hill coefficient' ); %requests input for hill coefficient
d2 = input('degradation rate '); %requests input for protein degradation rate
inducer = input('Concentration of inducer '); %requests input for concentration of inducer used
time = input('timepoints '); %requests number of time points

[t,y] = ode45(@induced,[0 time],0);