An insulin pump system for continuous monitoring and delivery of insulin to a subject, wherein the insulin pump system comprises a controller, wherein the controller is adapted to perform the following steps: obtaining values associated with blood levels of insulin for a subject; mapping the data using a transfer function; generating a linear difference model comprising a plurality of states; obtaining a defined value zone for the insulin in the subject; calculating a next administration dose and/or time of delivery of the insulin based on a predicted insulin value using the linear difference model and the defined zone; delivering the insulin to the subject based upon the calculated next administration; wherein the obtaining values comprises: obtaining insulin and meal data values and continuous glucose monitoring data values for the subject; mapping the insulin administration and meal data values using transfer functions: I map s = 1 τ 1 s + 1 τ 2 s + 1 I s M map s = 1 τ 3 s τ 4 s + 1 M s to obtain new states I map and M map, wherein τ are time points of measurement; wherein the generating a linear difference model comprises generating a linear difference model comprising: G k + 1 = α 1 G k + ⋯ + α p G k − p + β 11 I map, k − d 1 − 1 + ⋯ + β 1 q 1 I map, k − d 1 − q 1 + ⋯ + β 21 M map, k − d 2 − 1 + ⋯ + β 2 q 2 M map, k − d 2 − q 2 I map, k + 1 = γ 1 I map, k + γ 2 I map, k − 1 + γ 3 I k + γ 4 I k − 1 M map, k + 1 = δ 1 M map, k + δ 2 M map, k − 1 + δ 3 M k + δ 4 M k − 1, wherein G represents glucose blood concentration, α and β are model coefficients, d 1 and d 2 are insulin and meal time delays, respectively, and p, q 1, and q 2 are the orders of glucose, insulin, and meal, respectively; and y and δ represent weighting factors for insulin and meal respectively after being absorbed into the blood; and wherein obtaining a defined value zone comprises obtaining a predefined glycemic zone for the subject. obtaining values associated with blood levels of insulin for a subject; mapping the data using a transfer function; generating a linear difference model comprising a plurality of states; obtaining a defined value zone for the insulin in the subject; calculating a next administration dose and/or time of delivery of the insulin based on a predicted insulin value using the linear difference model and the defined zone; delivering the insulin to the subject based upon the calculated next administration; wherein the obtaining values comprises: obtaining insulin and meal data values and continuous glucose monitoring data values for the subject; mapping the insulin administration and meal data values using transfer functions: I map s = 1 τ 1 s + 1 τ 2 s + 1 I s M map s = 1 τ 3 s τ 4 s + 1 M s to obtain new states I map and M map, wherein τ are time points of measurement; wherein the generating a linear difference model comprises generating a linear difference model comprising: G k + 1 = α 1 G k + ⋯ + α p G k − p + β 11 I map, k − d 1 − 1 + ⋯ + β 1 q 1 I map, k − d 1 − q 1 + ⋯ + β 21 M map, k − d 2 − 1 + ⋯ + β 2 q 2 M map, k − d 2 − q 2 I map, k + 1 = γ 1 I map, k + γ 2 I map, k − 1 + γ 3 I k + γ 4 I k − 1 M map, k + 1 = δ 1 M map, k + δ 2 M map, k − 1 + δ 3 M k + δ 4 M k − 1, wherein G represents glucose blood concentration, α and β are model coefficients, d 1 and d 2 are insulin and meal time delays, respectively, and p, q 1, and q 2 are the orders of glucose, insulin, and meal, respectively; and y and δ represent weighting factors for insulin and meal respectively after being absorbed into the blood; and wherein obtaining a defined value zone comprises obtaining a predefined glycemic zone for the subject. obtaining insulin and meal data values and continuous glucose monitoring data values for the subject; mapping the insulin administration and meal data values using transfer functions: I map s = 1 τ 1 s + 1 τ 2 s + 1 I s M map s = 1 τ 3 s τ 4 s + 1 M s to obtain new states I map and M map, wherein τ are time points of measurement; wherein the generating a linear difference model comprises generating a linear difference model comprising: G k + 1 = α 1 G k + ⋯ + α p G k − p + β 11 I map, k − d 1 − 1 + ⋯ + β 1 q 1 I map, k − d 1 − q 1 + ⋯ + β 21 M map, k − d 2 − 1 + ⋯ + β 2 q 2 M map, k − d 2 − q 2 I map, k + 1 = γ 1 I map, k + γ 2 I map, k − 1 + γ 3 I k + γ 4 I k − 1 M map, k + 1 = δ 1 M map, k + δ 2 M map, k − 1 + δ 3 M k + δ 4 M k − 1, wherein G represents glucose blood concentration, α and β are model coefficients, d 1 and d 2 are insulin and meal time delays, respectively, and p, q 1, and q 2 are the orders of glucose, insulin, and meal, respectively; and y and δ represent weighting factors for insulin and meal respectively after being absorbed into the blood; and wherein obtaining a defined value zone comprises obtaining a predefined glycemic zone for the subject.