From 4d08c162d5ce4ac8b7eaeb9d985bfef62c886e99 Mon Sep 17 00:00:00 2001 From: Matt Argao Date: Sun, 31 Mar 2013 22:35:49 -0400 Subject: Added descriptions to each function. --- solution_stability.m | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) (limited to 'solution_stability.m') diff --git a/solution_stability.m b/solution_stability.m index 324621f..359335c 100644 --- a/solution_stability.m +++ b/solution_stability.m @@ -1,12 +1,13 @@ -function [ hist_h, hist_x ] = solution_stability( q_0, q_final, x_final, optics_positions, optics_focal_length, lambda, n_hist ) -%SOLUTION_STABILITY Summary of this function goes here -% Detailed explanation goes here +function [ hist_h, hist_x ] = solution_stability( q_0, q_final, x_final, optics_positions, optics_focal_length, lambda, n_hist, self_flag ) +%Plots histogram of n_hist number of sample points +%Area under plot is equal to the stability of the solution. + lens_displacement = 10^-3; x = lens_displacement*randn(n_hist,3); for i=1:n_hist xd= x(i,:); - fitness_array(i) = fitness( q_0, q_final, x_final, optics_positions + xd, optics_focal_length, lambda ); + fitness_array(i) = fitness( q_0, q_final, x_final, optics_positions + xd, optics_focal_length, lambda, self_flag ); end [hist_h, hist_x] = hist(fitness_array,100); -- cgit v1.2.3