Microscope [?]

clear all; close all; clc;

f1 = 3.5;
f2 = 6;
% f1 = 3.5;

%     figure()

s1 = [10:1:15];
M = zeros(1,length(s1));

L1 = f2+(s1(1)*f1)/(f1-s1(1));  %Boundary limit
L2 = f2+(s1(end)*f1)/(f1-s1(end));  %Boundary limit

Lmin = min(L1,L2);
Lmax = max(L1,L2);


RangeL = [Lmin*1000:100:Lmax*1000];

%LRange = [615:100:1756];

for j = RangeL
    L = j/1000;
    for i = 1:length(s1)
        num = L*f2*(f1-s1(i))-s1(i)*f1*f2;
        den = s1(i)*(f1*(f2-L+s1(i))+s1(i)*(L-f2));
    %     s1
        M(i) = -num/den;
    end
    figure()
    plot(s1,M);
    title(['F1 = 3.5cm, F2 = 6cm, Length = ' num2str(L),'cm']);
    ylabel('Magnification, times');
    xlabel('Distance to object, cm');
    grid on
%     hold on

end
%figure()

% hold on
% hold off
% legend(['L=' num2str(L)]);
% legend('L = 0.615','L = 0.715','L = 0.815','L = 0.915','L = 1.015','L = 1.115','L = 1.215','L = 1.315','L = 1.415','L = 1.515','L = 1.615','L = 1.715')


%
% legend('3,5','4','4.5','5');