from math import exp def d_pre(A, alpha, h, t): d = A*((exp(-t/alpha))/(1+(t**h))) return d x = d_pre(0.88,0.23, 0.38, .5)

use this function to perform a full forward model.

Calculate the time series for one of your sets of parameter values from the brute force example from 0.01 to 0.51 seconds at 0.01 second intervals.

Plot this using Matplotlib with labels of time (s) on the abscissa and predicted response (unitless) on the ordinate.