Creating Function that Divides of Mean Pixel value fo two images

I'm trying to create a function that divides the mean pixel values of two images. I tried to write a function, but the answer it's giving me is one which is incorrect because when I calculated the individual mean pixel values and divided them I got approx 1/2. Here is my code:

function [Beta] = betacalc(beta1,beta2)
 beta1im = imread(beta1);
 beta2im = imread(beta2);
 beta1mean = mean(beta1im(:));
 beta2mean = mean(beta1im(:));
 betadivide = beta1mean/beta2mean 
end

Answer 1

[Image Analyst]

The code is correct, except for the error that Beta is never assigned by your function.