Программа по расчету сдачи в кассе

% Напишите программу, которая определяет сдачу в кассе самообслуживания
% супермаркета для покупок до 50$. Программа генерирует случайное число
% между 0.01 и 50.00 и выводит на экран это число, как сумму оплаты
% покупки. Затем программа просит, чтобы пользователь ввел оплату одной
% купюрой 1$, 5$, 10$, 20$ или 50$. Если этот платеж меньше, чем сумма
% оплаты, то на экран выводится сообщение об ошибке. Если платеж
% достаточен, программа вычисляет сдачу, как список купюр и/или монет для
% сдачи, который должен быть составлен из наименьшего числа купюр и монет.
% Например, если сумма покупки составляет 2.33$, и купюра 10$ вводится, как
% платеж, то сдача-одна купюра 5$, две купюры по 1$, две четверти(по 25 
% центов), десять центов, один никель (5 центов) и два пенса.
format bank
a=0.01+(50-0.01)*rand
b=input('введите оплату одной купюрой, 1$, 5$, $10, $20, $50 ');
if a>b
    disp('ОШИБКА, сумма оплаты меньше стоимости покупки')
else
    c=[0.01 0.05 0.25 1 5 10 20];
    cc=(b-a)-c;
end

format bank
a=0.01+(50-0.01)*rand
b=input('введите оплату одной купюрой, 1$, 5$, $10, $20, $50 ');
if a>b
    disp('ОШИБКА, сумма оплаты меньше стоимости покупки')
else
    c=[0.01 0.05 0.25 1 5 10 20];
    cc=(b-a)-c;
    j=0;
    while cc~=[]
    cc(cc<=0)=[];
    i=length(cc);
    j=j+1;
    d(j)=c(i);
    end
end
d

Error using  ~= 
Matrix dimensions must agree.
 
Error in e6828 (line 21)
    while cc~=[]

c=[0.01 0.05 0.25 1 5 10 20];
a=5.451
b=10
if a>b
    disp('ОШИБКА, сумма оплаты меньше стоимости покупки')
else
    n=length(c);    
    cc1=a-floor(a/c(1))*c(1);
    cc=b-a+cc1+c(1)/2;
    for i=n:-1:1
        d(i)=floor(cc/c(i));
        cc=cc-d(i)*c(i);
    end
    c
    d
    str=sprintf('Скидка %.4f$',cc1);
    disp(str);
end