
#
# This file is automatically generated with 
# the System Biology Format Converter (http://sbfc.sourceforge.net/)
# from an SBML file.
#
# The conversion system has the following limitations:
#  - You may have to re order some reactions and Assignment Rules definition
#  - Delays are not taken into account
#  - You should change the lsode parameters (start, end, steps) to get better results
#

#
# The following line is there to be sure that Octave think that this file 
# is a script and not function file
#
1;

#
# Model name = Chatterjee2010_BloodCoagulation
#
# is urn:miriam:biomodels.db:MODEL1108260014
# isDescribedBy urn:miriam:pubmed:20941387
#
function z=pow(x,y),z=x^y;endfunction
function z=root(x,y),z=y^(1/x);endfunction
function z = piecewise(varargin)
	numArgs = nargin;
	result = 0;
	foundResult = 0;
	for k=1:2: numArgs-1
		if varargin{k+1} == 1
			result = varargin{k};
			foundResult = 1;
			break;
		endif
	end
	if foundResult == 0
		result = varargin{numArgs};
	endif
z = result;
endfunction

function z=function_function_1(k1,substrate1,substrate2,k2,epsilon,product), z=(k1*substrate1*substrate2-k2*product/epsilon);endfunction

function xdot=f(x,t)
# Compartment: id = compartment_1, name = compartment_1
# Warning compartment compartment_1 has no rule and is not constant
# Parameter:   id =  parameter_1, name = IIa*
	global_par_parameter_1=0.0;
# Parameter:   id =  parameter_2, name = f_IIa*
# Parameter:   id =  parameter_3, name = epsilon
# Parameter:   id =  parameter_4, name = epsilon_max
# Parameter:   id =  parameter_5, name = k
	global_par_parameter_5=0.0050;
# Parameter:   id =  parameter_6, name = epsilon_max_0
	global_par_parameter_6=0.01;
# assignmentRule: variable = parameter_2
	global_par_parameter_2=global_par_parameter_1^1.6123/(global_par_parameter_1^1.6123+(2.4279E-9)^1.6123);
# assignmentRule: variable = parameter_4
	global_par_parameter_4=global_par_parameter_6+(1-global_par_parameter_6)*global_par_parameter_2;
# rateRule: variable = parameter_3
global_par_parameter_3 = x(82);
# Reaction: id = reaction_1, name = R1

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_1_k1=3200000.0;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_1_k2=0.031;

	reaction_reaction_1=compartment_compartment_1*(reaction_reaction_1_k1*x(1)*x(3)-reaction_reaction_1_k2*x(2));
# Reaction: id = reaction_2, name = R2

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_2_k1=2.3E7;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_2_k2=3.1E-5;

	reaction_reaction_2=compartment_compartment_1*(reaction_reaction_2_k1*x(1)*x(5)-reaction_reaction_2_k2*x(4));
# Reaction: id = reaction_3, name = R3

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_3_k1=440000.0;

	reaction_reaction_3=compartment_compartment_1*reaction_reaction_3_k1*x(4)*x(3);
# Reaction: id = reaction_4, name = R4

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_4_k1=1.3E7;

	reaction_reaction_4=compartment_compartment_1*reaction_reaction_4_k1*x(6)*x(3);
# Reaction: id = reaction_5, name = R5

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_5_k1=23000.0;

	reaction_reaction_5=compartment_compartment_1*reaction_reaction_5_k1*x(7)*x(3);
# Reaction: id = reaction_6, name = R6

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_6_k1=2.5E7;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_6_k2=0.0105;

	reaction_reaction_6=compartment_compartment_1*(reaction_reaction_6_k1*x(4)*x(9)-reaction_reaction_6_k2*x(8));
# Reaction: id = reaction_7, name = R6b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_7_k1=6.0;

	reaction_reaction_7=compartment_compartment_1*reaction_reaction_7_k1*x(8);
# Reaction: id = reaction_8, name = R7

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_8_k1=2.2E7;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_8_k2=19.0;

	reaction_reaction_8=compartment_compartment_1*(reaction_reaction_8_k1*x(4)*x(6)-reaction_reaction_8_k2*x(10));
# Reaction: id = reaction_9, name = R8

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_9_k1=1.0E7;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_9_k2=2.4;

	reaction_reaction_9=compartment_compartment_1*(reaction_reaction_9_k1*x(4)*x(11)-reaction_reaction_9_k2*x(12));
# Reaction: id = reaction_10, name = R8b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_10_k1=1.8;

	reaction_reaction_10=compartment_compartment_1*reaction_reaction_10_k1*x(12);
# Reaction: id = reaction_11, name = R9

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_11_k1=7500.0;

	reaction_reaction_11=compartment_compartment_1*reaction_reaction_11_k1*x(14)*x(6);
# Reaction: id = reaction_12, name = R10

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_12_k1=2.0E7;

	reaction_reaction_12=compartment_compartment_1*reaction_reaction_12_k1*x(7)*x(15);
# Reaction: id = reaction_13, name = R11

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_13_k1=1.0E7;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_13_k2=1.0E-4;

	reaction_reaction_13=compartment_compartment_1*function_function_1(reaction_reaction_13_k1, x(16), x(13), reaction_reaction_13_k2, global_par_parameter_3, x(17));
# Reaction: id = reaction_14, name = R12

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_14_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_14_k2=1.0E-5;

	reaction_reaction_14=compartment_compartment_1*function_function_1(reaction_reaction_14_k1, x(17), x(9), reaction_reaction_14_k2, global_par_parameter_3, x(18));
# Reaction: id = reaction_15, name = R12b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_15_k1=8.2;

	reaction_reaction_15=compartment_compartment_1*reaction_reaction_15_k1*x(18);
# Reaction: id = reaction_16, name = R13

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_16_k1=22000.0;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_16_k2=6.0E-5;

	reaction_reaction_16=compartment_compartment_1*function_function_1(reaction_reaction_16_k1, x(19), x(20), reaction_reaction_16_k2, global_par_parameter_3, x(16));
# Reaction: id = reaction_17, name = R14

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_17_k1=0.0010;

	reaction_reaction_17=compartment_compartment_1*reaction_reaction_17_k1*x(18);
# Reaction: id = reaction_18, name = R15

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_18_k1=0.0010;

	reaction_reaction_18=compartment_compartment_1*reaction_reaction_18_k1*x(17);
# Reaction: id = reaction_19, name = R16

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_19_k1=2.0E7;

	reaction_reaction_19=compartment_compartment_1*reaction_reaction_19_k1*x(7)*x(21);
# Reaction: id = reaction_20, name = R17

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_20_k1=4.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_20_k2=0.0080;

	reaction_reaction_20=compartment_compartment_1*function_function_1(reaction_reaction_20_k1, x(6), x(22), reaction_reaction_20_k2, global_par_parameter_3, x(23));
# Reaction: id = reaction_21, name = R18

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_21_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_21_k2=2.06;

	reaction_reaction_21=compartment_compartment_1*function_function_1(reaction_reaction_21_k1, x(23), x(14), reaction_reaction_21_k2, global_par_parameter_3, x(79));
# Reaction: id = reaction_22, name = R18b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_22_k1=63.5;

	reaction_reaction_22=compartment_compartment_1*reaction_reaction_22_k1*x(79);
# Reaction: id = reaction_23, name = R19

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_23_k1=1.5E7;

	reaction_reaction_23=compartment_compartment_1*reaction_reaction_23_k1*x(23)*x(25);
# Reaction: id = reaction_24, name = R20

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_24_k1=900000.0;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_24_k2=3.6E-4;

	reaction_reaction_24=compartment_compartment_1*function_function_1(reaction_reaction_24_k1, x(6), x(26), reaction_reaction_24_k2, global_par_parameter_3, x(27));
# Reaction: id = reaction_25, name = R21

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_25_k1=3.2E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_25_k2=0.011;

	reaction_reaction_25=compartment_compartment_1*(reaction_reaction_25_k1*x(10)*x(26)-reaction_reaction_25_k2*x(28));
# Reaction: id = reaction_26, name = R22

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_26_k1=5.0E7;

	reaction_reaction_26=compartment_compartment_1*reaction_reaction_26_k1*x(4)*x(27);
# Reaction: id = reaction_27, name = R23

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_27_k1=1500.0;

	reaction_reaction_27=compartment_compartment_1*reaction_reaction_27_k1*x(6)*x(29);
# Reaction: id = reaction_28, name = R24

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_28_k1=7100.0;

	reaction_reaction_28=compartment_compartment_1*reaction_reaction_28_k1*x(25)*x(29);
# Reaction: id = reaction_29, name = R25

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_29_k1=490.0;

	reaction_reaction_29=compartment_compartment_1*reaction_reaction_29_k1*x(13)*x(29);
# Reaction: id = reaction_30, name = R26

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_30_k1=7100.0;

	reaction_reaction_30=compartment_compartment_1*reaction_reaction_30_k1*x(7)*x(29);
# Reaction: id = reaction_31, name = R27

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_31_k1=230.0;

	reaction_reaction_31=compartment_compartment_1*reaction_reaction_31_k1*x(4)*x(29);
# Reaction: id = reaction_32, name = R28

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_32_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_32_k2=6100.0;

	reaction_reaction_32=compartment_compartment_1*(reaction_reaction_32_k1*x(36)*x(7)-reaction_reaction_32_k2*x(37));
# Reaction: id = reaction_33, name = R28b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_33_k1=53.8;

	reaction_reaction_33=compartment_compartment_1*reaction_reaction_33_k1*x(37);
# Reaction: id = reaction_34, name = R29

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_34_k1=5.0E-4;

	reaction_reaction_34=compartment_compartment_1*reaction_reaction_34_k1*x(39);
# Reaction: id = reaction_35, name = R30

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_35_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_35_k2=750.0;

	reaction_reaction_35=compartment_compartment_1*function_function_1(reaction_reaction_35_k1, x(40), x(39), reaction_reaction_35_k2, global_par_parameter_3, x(41));
# Reaction: id = reaction_36, name = R30b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_36_k1=0.033;

	reaction_reaction_36=compartment_compartment_1*reaction_reaction_36_k1*x(41);
# Reaction: id = reaction_37, name = R31

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_37_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_37_k2=3600.0;

	reaction_reaction_37=compartment_compartment_1*function_function_1(reaction_reaction_37_k1, x(40), x(42), reaction_reaction_37_k2, global_par_parameter_3, x(43));
# Reaction: id = reaction_38, name = R31b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_38_k1=40.0;

	reaction_reaction_38=compartment_compartment_1*reaction_reaction_38_k1*x(43);
# Reaction: id = reaction_39, name = R32

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_39_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_39_k2=45.3;

	reaction_reaction_39=compartment_compartment_1*function_function_1(reaction_reaction_39_k1, x(39), x(44), reaction_reaction_39_k2, global_par_parameter_3, x(45));
# Reaction: id = reaction_40, name = R32b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_40_k1=5.7;

	reaction_reaction_40=compartment_compartment_1*reaction_reaction_40_k1*x(45);
# Reaction: id = reaction_41, name = R33

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_41_k1=27000.0;

	reaction_reaction_41=compartment_compartment_1*reaction_reaction_41_k1*x(42)*x(44);
# Reaction: id = reaction_42, name = R34

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_42_k1=0.011;

	reaction_reaction_42=compartment_compartment_1*reaction_reaction_42_k1*x(44);
# Reaction: id = reaction_43, name = R35

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_43_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_43_k2=2.4;

	reaction_reaction_43=compartment_compartment_1*(reaction_reaction_43_k1*x(40)*x(47)-reaction_reaction_43_k2*x(48));
# Reaction: id = reaction_44, name = R36

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_44_k1=3600.0;

	reaction_reaction_44=compartment_compartment_1*reaction_reaction_44_k1*x(40)*x(49);
# Reaction: id = reaction_45, name = R37

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_45_k1=21.6;

	reaction_reaction_45=compartment_compartment_1*reaction_reaction_45_k1*x(40)*x(29);
# Reaction: id = reaction_46, name = R38

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_46_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_46_k2=5.0;

	reaction_reaction_46=compartment_compartment_1*(reaction_reaction_46_k1*x(51)*x(7)-reaction_reaction_46_k2*x(52));
# Reaction: id = reaction_47, name = R38b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_47_k1=1.3E-4;

	reaction_reaction_47=compartment_compartment_1*reaction_reaction_47_k1*x(52);
# Reaction: id = reaction_48, name = R39

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_48_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_48_k2=200.0;

	reaction_reaction_48=compartment_compartment_1*function_function_1(reaction_reaction_48_k1, x(40), x(51), reaction_reaction_48_k2, global_par_parameter_3, x(54));
# Reaction: id = reaction_49, name = R39b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_49_k1=5.7E-4;

	reaction_reaction_49=compartment_compartment_1*reaction_reaction_49_k1*x(54);
# Reaction: id = reaction_50, name = R40

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_50_k1=3190000.0;

	reaction_reaction_50=compartment_compartment_1*reaction_reaction_50_k1*x(53)*x(51);
# Reaction: id = reaction_51, name = R41

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_51_k1=320.0;

	reaction_reaction_51=compartment_compartment_1*reaction_reaction_51_k1*x(53)*x(29);
# Reaction: id = reaction_52, name = R42

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_52_k1=1800.0;

	reaction_reaction_52=compartment_compartment_1*reaction_reaction_52_k1*x(53)*x(49);
# Reaction: id = reaction_53, name = R43

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_53_k1=100.0;

	reaction_reaction_53=compartment_compartment_1*reaction_reaction_53_k1*x(53)*x(57);
# Reaction: id = reaction_54, name = R44

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_54_k1=4300.0;

	reaction_reaction_54=compartment_compartment_1*reaction_reaction_54_k1*x(53)*x(59);
# Reaction: id = reaction_55, name = R45

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_55_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_55_k2=41.0;

	reaction_reaction_55=compartment_compartment_1*function_function_1(reaction_reaction_55_k1, x(53), x(11), reaction_reaction_55_k2, global_par_parameter_3, x(60));
# Reaction: id = reaction_56, name = R45b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_56_k1=7.7;

	reaction_reaction_56=compartment_compartment_1*reaction_reaction_56_k1*x(60);
# Reaction: id = reaction_57, name = R46

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_57_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_57_k2=0.64;

	reaction_reaction_57=compartment_compartment_1*function_function_1(reaction_reaction_57_k1, x(13), x(9), reaction_reaction_57_k2, global_par_parameter_3, x(61));
# Reaction: id = reaction_58, name = R46b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_58_k1=7.0E-4;

	reaction_reaction_58=compartment_compartment_1*reaction_reaction_58_k1*x(61);
# Reaction: id = reaction_59, name = R47

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_59_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_59_k2=2.1;

	reaction_reaction_59=compartment_compartment_1*function_function_1(reaction_reaction_59_k1, x(6), x(15), reaction_reaction_59_k2, global_par_parameter_3, x(62));
# Reaction: id = reaction_60, name = R47b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_60_k1=0.023;

	reaction_reaction_60=compartment_compartment_1*reaction_reaction_60_k1*x(62);
# Reaction: id = reaction_61, name = R48

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_61_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_61_k2=0.9;

	reaction_reaction_61=compartment_compartment_1*(reaction_reaction_61_k1*x(5)*x(11)-reaction_reaction_61_k2*x(63));
# Reaction: id = reaction_62, name = R48b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_62_k1=3.6E-5;

	reaction_reaction_62=compartment_compartment_1*reaction_reaction_62_k1*x(63);
# Reaction: id = reaction_63, name = R49

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_63_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_63_k2=210.0;

	reaction_reaction_63=compartment_compartment_1*function_function_1(reaction_reaction_63_k1, x(5), x(9), reaction_reaction_63_k2, global_par_parameter_3, x(64));
# Reaction: id = reaction_64, name = R49b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_64_k1=1.6E-6;

	reaction_reaction_64=compartment_compartment_1*reaction_reaction_64_k1*x(64);
# Reaction: id = reaction_65, name = R50

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_65_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_65_k2=636.0;

	reaction_reaction_65=compartment_compartment_1*(reaction_reaction_65_k1*x(65)*x(7)-reaction_reaction_65_k2*x(66));
# Reaction: id = reaction_66, name = R50b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_66_k1=84.0;

	reaction_reaction_66=compartment_compartment_1*reaction_reaction_66_k1*x(66);
# Reaction: id = reaction_67, name = R51

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_67_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_67_k2=742.6;

	reaction_reaction_67=compartment_compartment_1*(reaction_reaction_67_k1*x(68)*x(7)-reaction_reaction_67_k2*x(69));
# Reaction: id = reaction_68, name = R51b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_68_k1=7.4;

	reaction_reaction_68=compartment_compartment_1*reaction_reaction_68_k1*x(69);
# Reaction: id = reaction_69, name = R52

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_69_k1=1000000.0;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_69_k2=0.064;

	reaction_reaction_69=compartment_compartment_1*(reaction_reaction_69_k1*x(68)^2-reaction_reaction_69_k2*x(72));
# Reaction: id = reaction_70, name = R53

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_70_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_70_k2=701.0;

	reaction_reaction_70=compartment_compartment_1*(reaction_reaction_70_k1*x(72)*x(7)-reaction_reaction_70_k2*x(73));
# Reaction: id = reaction_71, name = R53b

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_71_k1=49.0;

	reaction_reaction_71=compartment_compartment_1*reaction_reaction_71_k1*x(73);
# Reaction: id = reaction_72, name = R54

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_72_k1=1.0E8;
	# Local Parameter:   id =  k2, name = k2
	reaction_reaction_72_k2=1000.0;

	reaction_reaction_72=compartment_compartment_1*(reaction_reaction_72_k1*x(71)*x(7)-reaction_reaction_72_k2*x(75));
# Reaction: id = reaction_73, name = R55

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_73_k1=16000.0;

	reaction_reaction_73=compartment_compartment_1*reaction_reaction_73_k1*x(73)*x(29);
# Reaction: id = reaction_74, name = R56

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_74_k1=16000.0;

	reaction_reaction_74=compartment_compartment_1*reaction_reaction_74_k1*x(69)*x(29);
# Reaction: id = reaction_75, name = R57

	# Local Parameter:   id =  k1, name = k1
	reaction_reaction_75_k1=10000.0;

	reaction_reaction_75=compartment_compartment_1*reaction_reaction_75_k1*x(75)*x(29);

#Event: id=event_1
	event_event_1=sin(time*100) > 0;

	if(event_event_1) 
		global_par_parameter_1=piecewise(x(7), global_par_parameter_1 < x(7), global_par_parameter_1);
	endif

#Event: id=event_2
	event_event_2=time > 0;

	if(event_event_2) 
		compartment_compartment_1=compartment_compartment_1*5;
	endif
	xdot=zeros(82,1);
	# rateRule: variable = parameter_3
	xdot(82) = global_par_parameter_5*(global_par_parameter_4-global_par_parameter_3);

	# Species:   id = species_1, name = TF, affected by kineticLaw

	xdot(1) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_1) + (-1.0 * reaction_reaction_2));

	# Species:   id = species_2, name = TF_VII, affected by kineticLaw

	xdot(2) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_1));

	# Species:   id = species_3, name = VII, affected by kineticLaw

	xdot(3) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_1) + (-1.0 * reaction_reaction_3) + (-1.0 * reaction_reaction_4) + (-1.0 * reaction_reaction_5));

	# Species:   id = species_4, name = TF_VIIa, affected by kineticLaw

	xdot(4) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_2) + (-1.0 * reaction_reaction_3) + ( 1.0 * reaction_reaction_3) + (-1.0 * reaction_reaction_6) + (-1.0 * reaction_reaction_8) + (-1.0 * reaction_reaction_9) + ( 1.0 * reaction_reaction_10) + (-1.0 * reaction_reaction_26) + (-1.0 * reaction_reaction_31));

	# Species:   id = species_5, name = VIIa, affected by kineticLaw

	xdot(5) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_2) + ( 1.0 * reaction_reaction_3) + ( 1.0 * reaction_reaction_4) + ( 1.0 * reaction_reaction_5) + (-1.0 * reaction_reaction_61) + ( 1.0 * reaction_reaction_62) + (-1.0 * reaction_reaction_63) + ( 1.0 * reaction_reaction_64));

	# Species:   id = species_6, name = Xa, affected by kineticLaw

	xdot(6) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_4) + ( 1.0 * reaction_reaction_4) + (-1.0 * reaction_reaction_8) + (-1.0 * reaction_reaction_11) + ( 1.0 * reaction_reaction_11) + ( 1.0 * reaction_reaction_15) + (-1.0 * reaction_reaction_20) + (-1.0 * reaction_reaction_24) + (-1.0 * reaction_reaction_27) + ( 1.0 * reaction_reaction_58) + (-1.0 * reaction_reaction_59) + ( 1.0 * reaction_reaction_60) + ( 1.0 * reaction_reaction_64));

	# Species:   id = species_7, name = IIa, affected by kineticLaw

	xdot(7) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_5) + ( 1.0 * reaction_reaction_5) + ( 1.0 * reaction_reaction_11) + (-1.0 * reaction_reaction_12) + ( 1.0 * reaction_reaction_12) + (-1.0 * reaction_reaction_19) + ( 1.0 * reaction_reaction_19) + ( 1.0 * reaction_reaction_23) + (-1.0 * reaction_reaction_30) + (-1.0 * reaction_reaction_32) + ( 1.0 * reaction_reaction_33) + (-1.0 * reaction_reaction_46) + ( 1.0 * reaction_reaction_47) + (-1.0 * reaction_reaction_65) + ( 1.0 * reaction_reaction_66) + (-1.0 * reaction_reaction_67) + ( 1.0 * reaction_reaction_68) + (-1.0 * reaction_reaction_70) + ( 1.0 * reaction_reaction_71) + (-1.0 * reaction_reaction_72));

	# Species:   id = species_8, name = TF_VIIa_X, affected by kineticLaw

	xdot(8) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_6) + (-1.0 * reaction_reaction_7));

	# Species:   id = species_9, name = X, affected by kineticLaw

	xdot(9) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_6) + (-1.0 * reaction_reaction_14) + ( 1.0 * reaction_reaction_17) + (-1.0 * reaction_reaction_57) + (-1.0 * reaction_reaction_63));

	# Species:   id = species_10, name = TF_VIIa_Xa, affected by kineticLaw

	xdot(10) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_7) + ( 1.0 * reaction_reaction_8) + (-1.0 * reaction_reaction_25));

	# Species:   id = species_11, name = IX, affected by kineticLaw

	xdot(11) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_9) + (-1.0 * reaction_reaction_55) + (-1.0 * reaction_reaction_61));

	# Species:   id = species_12, name = TF_VIIa_IX, affected by kineticLaw

	xdot(12) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_9) + (-1.0 * reaction_reaction_10));

	# Species:   id = species_13, name = IXa, affected by kineticLaw

	xdot(13) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_10) + (-1.0 * reaction_reaction_13) + ( 1.0 * reaction_reaction_17) + ( 1.0 * reaction_reaction_18) + (-1.0 * reaction_reaction_29) + ( 1.0 * reaction_reaction_56) + (-1.0 * reaction_reaction_57) + ( 1.0 * reaction_reaction_58) + ( 1.0 * reaction_reaction_62));

	# Species:   id = species_14, name = II, affected by kineticLaw

	xdot(14) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_11) + (-1.0 * reaction_reaction_21));

	# Species:   id = species_15, name = VIII, affected by kineticLaw

	xdot(15) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_12) + (-1.0 * reaction_reaction_59));

	# Species:   id = species_16, name = VIIIa, affected by kineticLaw

	xdot(16) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_12) + (-1.0 * reaction_reaction_13) + ( 1.0 * reaction_reaction_16) + ( 1.0 * reaction_reaction_60));

	# Species:   id = species_17, name = IXa_VIIIa, affected by kineticLaw

	xdot(17) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_13) + (-1.0 * reaction_reaction_14) + ( 1.0 * reaction_reaction_15) + (-1.0 * reaction_reaction_18));

	# Species:   id = species_18, name = IXa_VIIIa_X, affected by kineticLaw

	xdot(18) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_14) + (-1.0 * reaction_reaction_15) + (-1.0 * reaction_reaction_17));

	# Species:   id = species_19, name = VIIIa1_L, affected by kineticLaw

	xdot(19) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_16) + ( 1.0 * reaction_reaction_17) + ( 1.0 * reaction_reaction_18));

	# Species:   id = species_20, name = VIIIa2, affected by kineticLaw

	xdot(20) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_16) + ( 1.0 * reaction_reaction_17) + ( 1.0 * reaction_reaction_18));

	# Species:   id = species_21, name = V, affected by kineticLaw

	xdot(21) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_19));

	# Species:   id = species_22, name = Va, affected by kineticLaw

	xdot(22) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_19) + (-1.0 * reaction_reaction_20));

	# Species:   id = species_23, name = Xa_Va, affected by kineticLaw

	xdot(23) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_20) + (-1.0 * reaction_reaction_21) + ( 1.0 * reaction_reaction_22) + (-1.0 * reaction_reaction_23) + ( 1.0 * reaction_reaction_23));

	# Species:   id = species_24, name = Xa_Va_IIa
# Warning species is not changed by either rules or reactions

	xdot(24) = ;

	# Species:   id = species_25, name = mIIa, affected by kineticLaw

	xdot(25) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_22) + (-1.0 * reaction_reaction_23) + (-1.0 * reaction_reaction_28));

	# Species:   id = species_26, name = TFPI, affected by kineticLaw

	xdot(26) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_24) + (-1.0 * reaction_reaction_25));

	# Species:   id = species_27, name = Xa_TFPI, affected by kineticLaw

	xdot(27) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_24) + (-1.0 * reaction_reaction_26));

	# Species:   id = species_28, name = TF_VIIa_Xa_TFPI, affected by kineticLaw

	xdot(28) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_25) + ( 1.0 * reaction_reaction_26));

	# Species:   id = species_29, name = ATIII, affected by kineticLaw

	xdot(29) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_27) + (-1.0 * reaction_reaction_28) + (-1.0 * reaction_reaction_29) + (-1.0 * reaction_reaction_30) + (-1.0 * reaction_reaction_31) + (-1.0 * reaction_reaction_45) + (-1.0 * reaction_reaction_51) + (-1.0 * reaction_reaction_73) + (-1.0 * reaction_reaction_74) + (-1.0 * reaction_reaction_75));

	# Species:   id = species_30, name = Xa_ATIII, affected by kineticLaw

	xdot(30) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_27));

	# Species:   id = species_31, name = mIIa_ATIII, affected by kineticLaw

	xdot(31) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_28));

	# Species:   id = species_32, name = IXa_ATIII, affected by kineticLaw

	xdot(32) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_29));

	# Species:   id = species_33, name = IIa_ATIII, affected by kineticLaw

	xdot(33) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_30));

	# Species:   id = species_34, name = TF_VIIa_ATIII, affected by kineticLaw

	xdot(34) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_31));

	# Species:   id = species_35, name = AMC, affected by kineticLaw

	xdot(35) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_33));

	# Species:   id = species_36, name = Boc_VPR_MCA, affected by kineticLaw

	xdot(36) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_32));

	# Species:   id = species_37, name = Boc_VPR_MCA_IIa, affected by kineticLaw

	xdot(37) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_32) + (-1.0 * reaction_reaction_33));

	# Species:   id = species_38, name = Boc_VPR, affected by kineticLaw

	xdot(38) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_33));

	# Species:   id = species_39, name = XII, affected by kineticLaw

	xdot(39) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_34) + (-1.0 * reaction_reaction_35) + (-1.0 * reaction_reaction_39));

	# Species:   id = species_40, name = XIIa, affected by kineticLaw

	xdot(40) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_34) + (-1.0 * reaction_reaction_35) + ( 2.0 * reaction_reaction_36) + (-1.0 * reaction_reaction_37) + ( 1.0 * reaction_reaction_38) + ( 1.0 * reaction_reaction_40) + (-1.0 * reaction_reaction_43) + (-1.0 * reaction_reaction_44) + (-1.0 * reaction_reaction_45) + (-1.0 * reaction_reaction_48) + ( 1.0 * reaction_reaction_49));

	# Species:   id = species_41, name = XIIa_XII, affected by kineticLaw

	xdot(41) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_35) + (-1.0 * reaction_reaction_36));

	# Species:   id = species_42, name = PK, affected by kineticLaw

	xdot(42) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_37) + (-1.0 * reaction_reaction_41));

	# Species:   id = species_43, name = XIIa_PK, affected by kineticLaw

	xdot(43) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_37) + (-1.0 * reaction_reaction_38));

	# Species:   id = species_44, name = K, affected by kineticLaw

	xdot(44) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_38) + (-1.0 * reaction_reaction_39) + ( 1.0 * reaction_reaction_40) + (-1.0 * reaction_reaction_41) + ( 2.0 * reaction_reaction_41) + (-1.0 * reaction_reaction_42));

	# Species:   id = species_45, name = XII_K, affected by kineticLaw

	xdot(45) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_39) + (-1.0 * reaction_reaction_40));

	# Species:   id = species_46, name = K.inhibited, affected by kineticLaw

	xdot(46) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_42));

	# Species:   id = species_47, name = CTI, affected by kineticLaw

	xdot(47) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_43));

	# Species:   id = species_48, name = XIIa_CTI, affected by kineticLaw

	xdot(48) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_43));

	# Species:   id = species_49, name = C1inh, affected by kineticLaw

	xdot(49) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_44) + (-1.0 * reaction_reaction_52));

	# Species:   id = species_50, name = XIIa_ATIII, affected by kineticLaw

	xdot(50) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_45));

	# Species:   id = species_51, name = XI, affected by kineticLaw

	xdot(51) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_46) + (-1.0 * reaction_reaction_48) + (-1.0 * reaction_reaction_50));

	# Species:   id = species_52, name = XI_IIa, affected by kineticLaw

	xdot(52) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_46) + (-1.0 * reaction_reaction_47));

	# Species:   id = species_53, name = XIa, affected by kineticLaw

	xdot(53) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_47) + ( 1.0 * reaction_reaction_49) + (-1.0 * reaction_reaction_50) + ( 2.0 * reaction_reaction_50) + (-1.0 * reaction_reaction_51) + (-1.0 * reaction_reaction_52) + (-1.0 * reaction_reaction_53) + (-1.0 * reaction_reaction_54) + (-1.0 * reaction_reaction_55) + ( 1.0 * reaction_reaction_56));

	# Species:   id = species_54, name = XIIa_XI, affected by kineticLaw

	xdot(54) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_48) + (-1.0 * reaction_reaction_49));

	# Species:   id = species_55, name = XIa_C1inh, affected by kineticLaw

	xdot(55) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_52));

	# Species:   id = species_56, name = XIa_alpha1AT, affected by kineticLaw

	xdot(56) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_53));

	# Species:   id = species_57, name = alpha1AT, affected by kineticLaw

	xdot(57) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_53));

	# Species:   id = species_58, name = XIa_alpha2AP, affected by kineticLaw

	xdot(58) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_54));

	# Species:   id = species_59, name = alpha2AP, affected by kineticLaw

	xdot(59) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_54));

	# Species:   id = species_60, name = XIa_IX, affected by kineticLaw

	xdot(60) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_55) + (-1.0 * reaction_reaction_56));

	# Species:   id = species_61, name = IXa_X, affected by kineticLaw

	xdot(61) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_57) + (-1.0 * reaction_reaction_58));

	# Species:   id = species_62, name = Xa_VIII, affected by kineticLaw

	xdot(62) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_59) + (-1.0 * reaction_reaction_60));

	# Species:   id = species_63, name = VIIa_IX, affected by kineticLaw

	xdot(63) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_61) + (-1.0 * reaction_reaction_62));

	# Species:   id = species_64, name = VIIa_X, affected by kineticLaw

	xdot(64) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_63) + (-1.0 * reaction_reaction_64));

	# Species:   id = species_65, name = Fbg, affected by kineticLaw

	xdot(65) = (1/(compartment_compartment_1))*((-1.0 * reaction_reaction_65));

	# Species:   id = species_66, name = Fbg_IIa, affected by kineticLaw

	xdot(66) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_65) + (-1.0 * reaction_reaction_66));

	# Species:   id = species_67, name = FPA, affected by kineticLaw

	xdot(67) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_66));

	# Species:   id = species_68, name = Fbn1, affected by kineticLaw

	xdot(68) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_66) + (-1.0 * reaction_reaction_67) + (-2.0 * reaction_reaction_69));

	# Species:   id = species_69, name = Fbn1_IIa, affected by kineticLaw

	xdot(69) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_67) + (-1.0 * reaction_reaction_68) + (-1.0 * reaction_reaction_74));

	# Species:   id = species_70, name = FPB, affected by kineticLaw

	xdot(70) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_68) + ( 1.0 * reaction_reaction_71));

	# Species:   id = species_71, name = Fbn2, affected by kineticLaw

	xdot(71) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_68) + (-1.0 * reaction_reaction_72));

	# Species:   id = species_72, name = Fbn1x2, affected by kineticLaw

	xdot(72) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_69) + (-1.0 * reaction_reaction_70));

	# Species:   id = species_73, name = Fbn1x2_IIa, affected by kineticLaw

	xdot(73) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_70) + (-1.0 * reaction_reaction_71) + (-1.0 * reaction_reaction_73));

	# Species:   id = species_74, name = Fbn2x2, affected by kineticLaw

	xdot(74) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_71));

	# Species:   id = species_75, name = Fbn2_IIa, affected by kineticLaw

	xdot(75) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_72) + (-1.0 * reaction_reaction_75));

	# Species:   id = species_76, name = Fbn1x2_IIa_ATIII, affected by kineticLaw

	xdot(76) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_73));

	# Species:   id = species_77, name = Fbn1_IIa_ATIII, affected by kineticLaw

	xdot(77) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_74));

	# Species:   id = species_78, name = Fbn2_IIa_ATIII, affected by kineticLaw

	xdot(78) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_75));

	# Species:   id = species_79, name = Xa_Va_II, affected by kineticLaw

	xdot(79) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_21) + (-1.0 * reaction_reaction_22));

	# Species:   id = species_80, name = XIIa_C1inh, affected by kineticLaw

	xdot(80) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_44));

	# Species:   id = species_81, name = XIa_ATIII, affected by kineticLaw

	xdot(81) = (1/(compartment_compartment_1))*(( 1.0 * reaction_reaction_51));

endfunction

#Initial conditions vector
x0=zeros(82,1);
x0(82) = 0.01;
x0(1) = 5.0E-12;
x0(2) = 0.0;
x0(3) = 1.0E-8;
x0(4) = 0.0;
x0(5) = 1.0E-10;
x0(6) = 0.0;
x0(7) = 0.0;
x0(8) = 0.0;
x0(9) = 1.6E-7;
x0(10) = 0.0;
x0(11) = 9.0E-8;
x0(12) = 0.0;
x0(13) = 0.0;
x0(14) = 1.4E-6;
x0(15) = 7.0E-10;
x0(16) = 0.0;
x0(17) = 0.0;
x0(18) = 0.0;
x0(19) = 0.0;
x0(20) = 0.0;
x0(21) = 2.0E-8;
x0(22) = 0.0;
x0(23) = 0.0;
x0(24) = 0.0;
x0(25) = 0.0;
x0(26) = 2.5E-9;
x0(27) = 0.0;
x0(28) = 0.0;
x0(29) = 3.4E-6;
x0(30) = 0.0;
x0(31) = 0.0;
x0(32) = 0.0;
x0(33) = 0.0;
x0(34) = 0.0;
x0(35) = 0.0;
x0(36) = 1.0E-5;
x0(37) = 0.0;
x0(38) = 0.0;
x0(39) = 3.4E-7;
x0(40) = 0.0;
x0(41) = 0.0;
x0(42) = 4.5E-7;
x0(43) = 0.0;
x0(44) = 0.0;
x0(45) = 0.0;
x0(46) = 0.0;
x0(47) = 4.2E-6;
x0(48) = 0.0;
x0(49) = 2.5E-6;
x0(50) = 0.0;
x0(51) = 3.1E-8;
x0(52) = 0.0;
x0(53) = 0.0;
x0(54) = 0.0;
x0(55) = 0.0;
x0(56) = 0.0;
x0(57) = 4.5E-5;
x0(58) = 0.0;
x0(59) = 1.0E-6;
x0(60) = 0.0;
x0(61) = 0.0;
x0(62) = 0.0;
x0(63) = 0.0;
x0(64) = 0.0;
x0(65) = 9.0E-6;
x0(66) = 0.0;
x0(67) = 0.0;
x0(68) = 0.0;
x0(69) = 0.0;
x0(70) = 0.0;
x0(71) = 0.0;
x0(72) = 0.0;
x0(73) = 0.0;
x0(74) = 0.0;
x0(75) = 0.0;
x0(76) = 0.0;
x0(77) = 0.0;
x0(78) = 0.0;
x0(79) = 0.0;
x0(80) = 0.0;
x0(81) = 0.0;


#Creating linespace
t=linspace(0,90,100);

#Solving equations
x=lsode("f",x0,t);

#ploting the results
plot(t,x);
