how to save this new images using dicomwrite

I have a programm where i take 4 sets of images, do certain equations on them and after that i would like to save the new set of images in my program (R) to a folder without changing the matrices of R the problem is with dicomwrite (either it is giving me errors and if i change a bit in my program, it gives me another matrices for R from type of UNI and not double). if someone good solve this i would really appreciate it note : i am saving the matrices resultes to my workspace without any problem thanks in advance

   for c=1:sizeOfStruct1(1,2)

      name1 = fullfile(folder_name1,filesStruct1(c).name);
      name2 = fullfile(folder_name2,filesStruct2(c).name);     
      name3 = fullfile(folder_name3,filesStruct3(c).name);
      name4 = fullfile(folder_name4,filesStruct4(c).name);

 % Read image from the first fcheckNplotolder, image from the second one,
 % convert their matrices to double and store them in A and B.    

      inv1{c} = double(dicomread(name1));
      inv1_pha{c} = double(dicomread(name2));       
      inv2{c} = double(dicomread(name3));
      inv2_pha{c} = double(dicomread(name4));

    alpha{c} = ((inv1_pha{c}.*pi)./(4096)) - (pi./2);
    beta{c} = ((inv2_pha{c}.*pi)./(4096)) - (pi./2);

      metadata = dicominfo(name1); 

 inv1_com{c} = inv1{c}.* exp(i.*alpha{c});
 inv2_com{c} = inv2{c}.* exp(i.*beta{c});

      N{c} = inv1_com{c}.* conj(inv2_com{c});
      D{c} = abs(inv1_com{c}.^2) + abs(inv2_com{c}.^2);
      R{c} = N{c}./D{c};

      dicomwrite(R, fullfile('Results',filesStruct1(c).name));
  end
  save('R');

[tarek abousaleh]

this is the whole code

function [ output_args ] = phase_equation( input_args )

% To delete the folder (Results) if it exists so we do not have to do it manually.

if exist('Results', 'dir')==7 system((sprintf('rm -r Results'))); end

% To create the folder (Results) containing our results.

if exist('Results', 'dir')==0 system((sprintf('mkdir Results'))); end

% To give the path of two folders which we want do the equation on.

folder_name1 = uigetdir('/home/gast/Desktop/'); folder_name2 = uigetdir('/home/gast/Desktop/'); folder_name3 = uigetdir('/home/gast/Desktop/'); folder_name4 = uigetdir('/home/gast/Desktop/');

% To read the contains of those folders we've choosen.

system((sprintf('ls "%s" -v | grep .IMA$ > list1.txt',folder_name1))); system((sprintf('ls "%s" -v | grep .IMA$ > list2.txt',folder_name2))); system((sprintf('ls "%s" -v | grep .IMA$ > list3.txt',folder_name3))); system((sprintf('ls "%s" -v | grep .IMA$ > list4.txt',folder_name4)));

% To put the contains of the folders in lists.

fid1 = fopen('list1.txt'); tline1 = fgetl(fid1); i=1; while ischar(tline1) filesStruct1(i).name=tline1; tline1 = fgetl(fid1); i=i+1; end fclose(fid1);

fid2 = fopen('list2.txt'); tline2 = fgetl(fid2); j=1; while ischar(tline2) filesStruct2(j).name=tline2; tline2 = fgetl(fid2); j=j+1; end fclose(fid2);

fid3 = fopen('list3.txt'); tline3 = fgetl(fid3); k=1; while ischar(tline3) filesStruct3(k).name=tline3; tline3 = fgetl(fid3); k=k+1; end fclose(fid3);

fid4 = fopen('list4.txt'); tline4 = fgetl(fid4); t=1; while ischar(tline4) filesStruct4(t).name=tline4; tline4 = fgetl(fid4); t=t+1; end fclose(fid4);

% To stop and exit the programm in case the number of the files in the % first folder does not math the number of the files in the second one.

sizeOfStruct1 = size(filesStruct1); sizeOfStruct2 = size(filesStruct2); sizeOfStruct3 = size(filesStruct3); sizeOfStruct4 = size(filesStruct4);

 if sizeOfStruct1(1,2) ~= sizeOfStruct2(1,2)
     msgbox('Number of images is not the same in the chosen folders!');
     exit(1);
 end

 if sizeOfStruct3(1,2) ~= sizeOfStruct4(1,2)
     msgbox('Number of images is not the same in the chosen folders!');
     exit(1);
 end

 % Here we do any equation we want on the data we took.

 for c=1:sizeOfStruct1(1,2)

      name1 = fullfile(folder_name1,filesStruct1(c).name);
      name2 = fullfile(folder_name2,filesStruct2(c).name);     
      name3 = fullfile(folder_name3,filesStruct3(c).name);
      name4 = fullfile(folder_name4,filesStruct4(c).name);

 % Read image from the first fcheckNplotolder, image from the second one,
 % convert their matrices to double and store them in A and B.    

      inv1{c} = double(dicomread(name1));
      inv1_pha{c} = double(dicomread(name2));       
      inv2{c} = double(dicomread(name3));
      inv2_pha{c} = double(dicomread(name4));

    alpha{c} = ((inv1_pha{c}.*pi)./(4096)) - (pi./2);
    beta{c} = ((inv2_pha{c}.*pi)./(4096)) - (pi./2);

      metadata = dicominfo(name1); 

 inv1_com{c} = inv1{c}.* exp(i.*alpha{c});
 inv2_com{c} = inv2{c}.* exp(i.*beta{c});

      N{c} = inv1_com{c}.* conj(inv2_com{c});
      D{c} = abs(inv1_com{c}.^2) + abs(inv2_com{c}.^2);
      R{c} = N{c}./D{c};

    %  dicomwrite(R,'Results',double, metadata);

end

save(Results,'R');

end

[tarek abousaleh]

if it does not work like this because i am using a cell array then in this one , ia m not storing the results in my workspace but with dicomwrite it is changing from double to UNI, which leads to a change in the images, here is the code like this:

   for c=1:sizeOfStruct1(1,2)

     name1 = strcat(folder_name1,'/',filesStruct1(c).name);
     name2 = strcat(folder_name2,'/',filesStruct2(c).name);
     name3 = strcat(folder_name3,'/',filesStruct3(c).name);
     name4 = strcat(folder_name4,'/',filesStruct4(c).name);

 % Read image from the first fcheckNplotolder, image from the second one,
 % convert their matrices to double and store them in A and B.    

      inv1 = double(dicomread(name1));
      inv1_pha = double(dicomread(name2));       
      inv2 = double(dicomread(name3));
      inv2_pha = double(dicomread(name4));

    alpha = ((inv1_pha.*pi)./(4096)) - (pi./2);
    beta = ((inv2_pha.*pi)./(4096)) - (pi./2);

      metadata = dicominfo(name1); 

 inv1_com = inv1.* exp(i.*alpha);
 inv2_com = inv2.* exp(i.*beta);

      N = inv1_com.* conj(inv2_com);
      D = abs(inv1_com.^2) + abs(inv2_com.^2);
      R = N./D

      dicomwrite(R, (strcat('Results/',filesStruct1(c).name)), metadata);

 end

[tarek abousaleh]

i would appreciate any help so much cause i have been stuck with this problem for a while