hello everyone

I would like to model the growth phenomenon in the cardiac muscle according to the following equations below. the code I have written is the one where I calculate the growth tensor to form the elastic tensor.

But I have a quadrature error.

Has anyone ever modelled this phenomenon before? Please help.

where theta is the growth factor, lamda_e the elastic stretch and the lamda_h is the hydrostatic stretch.

here is the code function

def growthfunction (f_0, s_0, n_0, d,v):

t = np.linspace(0,1,4)

dt=0.9

N=len(t)

O_max=2

O_min=0

gam=0.1

lam=0.2

lamda_h=1

I = Identity(d)

F=I+v

O_0=1

iterration=0

err=1

F_g=O_0*I

for j in range(N):

```
F_gt=F_g.T
FF_g=inner(f_0,(F_g.T*F_g)*f_0)
lamda_g=sqrt(FF_g)
#lamda_e=(1/lamda_g)*(f_0*(F.T/O_0)*(F/O_0)*f_0)^0.5
O_1=O_0+dt*((1/gam)*(((O_max-O_0)/(O_max-O_min))**lam)*(lamda_g-lamda_h))
F_g=O_1*I
O_0=O_1
print ('\n growth steps ' ,j)
```

return F_g