clear all close all clc %% tf = 1:3000; t0 = 0; MuMax = 0.5; Mu = MuMax .* tf; Lambda = Mu.^3./tf; %% F = exp(2.*Lambda ./ Mu) .* normcdf(-sqrt(Lambda./(tf-t0)).*((tf-t0)./Mu + 1)); plot(tf, F)
