Coding Examples

def morphogen_gradient(x):
    return 100 * (1 - x**2)

function dydt = geneReg(t, y)
    dydt = zeros(2,1);
    dydt(1) = y(2) - y(1);
    dydt(2) = -y(1) + y(2)^2;
end

library(ggplot2)
ggplot(data, aes(x=stage, y=expression)) + geom_line()