1 36456
2 33761
3 38309
4 36638
5 36111
6 32272
7 27107
8 32067
9 26357
10 34666
11 30034
12 30354
13 27336
14 21941
15 26401
16 18164
17 26762
18 26991
19 26834
20 24589
21 19174
22 23886
23 24236
24 23924
25 22350
26 18574
27 20333
28 16072
29 20529
30 21957
31 19046
32 19522
33 19460
34 18938
35 18967
36 18591
37 18381
38 18252
39 18008
40 18105
41 18020
42 17478
43 17191
44 17297
45 16348
46 15877
47 15338
48 15034
49 14678
50 14376
51 14128
52 13959
53 13840
54 13818
55 13714
56 13588
57 13437
58 12905
59 13535
60 13354
61 13094
62 12850
63 12776
64 12529
65 12724
66 11794
67 11602
68 11459
69 11496
70 11515
71 11477
72 11422
73 10986
74 11047
75 11066
76 11054
77 11107
78 11230
79 11276
80 11831
81 12085
82 12331
83 12699
84 12900
85 13198
86 13780
87 14260
88 14638
89 15050
90 15241
91 15495
92 15682
93 15753
94 15785
95 16037
96 16311
97 16745
98 17185
99 17596
100 18371
101 19296
102 20228
103 21147
104 22270
105 23567
106 25138
107 26994
108 29335
109 31629
110 34295
111 37223
112 39537
113 42161
114 44673
115 47444
116 50497
117 53185
118 56213
119 58430
120 59287
121 61958
122 65145
123 68966
124 73304
125 78388
126 84160
127 93028
128 100766
129 107977
130 115974
131 124484
132 133929
133 143114
134 153113
135 163586
136 175723
137 188437
138 203913
139 218937
140 232675
141 248257
142 264852
143 281380
144 297279
145 309865
146 321233
147 330155
148 339945
149 349038
150 356787
151 364870
152 371115
153 373320
154 378505
155 381353
156 386099
157 390029
158 390727
159 392331
160 391812
161 388058
162 383159
163 376017
164 342794
165 354315
166 341022
167 328934
168 319438
169 307822
170 295473
171 283507
172 273656
173 263676
174 255247
175 245652
176 237330
177 228090
178 217638
179 205750
180 194948
181 185028
182 175174
183 142207
184 153274
185 145609
186 137948
187 130692
188 123236
189 117368
190 111601
191 105864
192 100068
193 94942
194 90090
195 85775
196 82090
197 77722
198 73923
199 69721
200 66320
201 63190
202 60615
203 58140
204 56512
205 54658
206 53118
207 51404
208 50151
209 49229
210 48427
211 47754
212 46939
213 46242
214 45588
215 44614
216 43704
217 43269
218 42963
219 42548
220 42080
221 41862
222 41643
223 41286
224 40827
225 40306
226 39743
227 39110
228 38461
229 38431
230 38527
231 38328
232 38587
233 38577
234 37975
235 38173
236 38031
237 38009
238 38209
239 38330
240 38541
241 39942
242 40227
time = data(:,1);
A_data= data(:,2);
beta0 = 0.5;
alpha0 = 0.004;
gamma0 = 0.1;
upsilon0 = 0.13;
epsilon0= 0.07;
lamda0= 0.1;
sigma0= 0.07;
kappa0= 0.03;
nu0 = 0.0001;
xi0 = 0.0002;
lb =[0,0,0,0,0,0,0,0,0,0]; ub = [1,1,1,1,1,1,1,1,1,1];
B0 = [beta0; alpha0; gamma0; upsilon0; epsilon0; lamda0; sigma0; kappa0; nu0; xi0 ];
options = optimoptions(@lsqnonlin,'Algorithm','trust-region-reflective');
[B,resnorm,RESIDUAL,exitflag,OUTPUT,LAMBDA,Jacobian] = lsqnonlin(@Kinetics,B0,time,A_data,lb,ub,options );
disp(B)
function Av = Kinetics(B,time);
plot(time,A_data,time,Av)
function A = Kinetics(B, t)
x0 = [1217378052,100,3,2,1,1,1,1];
[T,Av] = ode45(@DifEq, t, x0);
function dA = DifEq(t, x)
N = 1390000000;
pi = 70000;
zeta = 0.1;
eta = 0.2;
theta = 0.3;
iota = 0.3;
delta = 0.1;
rho = 0.5;
mu = 0.0000425;
beta = B(1);
alpha =B(2);
gamma = B(3);
upsilon =B(4);
epsilon = B(5);
lamda = B(6);
sigma = B(7);
kappa = B(8);
nu = B(9);
xi = B(10);
xdot = zeros(8,1);
xdot(1) = pi -beta*(zeta*x(3)+eta*x(4)+theta*x(5)+iota*x(6))*(x(1)/N) -mu*x(1);
xdot(2) = beta*(zeta*x(3)+eta*x(4)+theta*x(5)+iota*x(6))*(x(1)/N) -(delta+mu)*x(2);
xdot(3) = rho*delta*x(2)-(lamda+gamma+nu+mu)*x(3);
xdot(4) = (1-rho)*delta*x(2)-(sigma+kappa+mu)*x(4);
xdot(5) = lamda*x(3) + sigma*x(4)-(alpha+upsilon+mu)*x(5);
xdot(6) = alpha*x(6) + kappa*x(4)- (epsilon+xi+mu)*x(6);
xdot(7) = gamma*x(3) + upsilon*x(5) + epsilon*x(6);
xdot(8) = nu*x(3) + xi*x(6);
dA = xdot;
end
A = Av(:,1);
end
i got an error;
Error using lsqnonlin (line 190)
Invalid datatype. Options argument must be created with OPTIMOPTIONS.
Error in error_14bmodel (line 19)
[B,resnorm,RESIDUAL,exitflag,OUTPUT,LAMBDA,Jacobian] = lsqnonlin(@Kinetics,B0,time,A_data,lb,ub,options );
i am fresher of matlab .i requesting you sir please help anyone how to resolve this error
if possible please refer some covid 19 or malaria,dengue hepities b models sir
thank you sir