function result=main()
result1=0;
result2=0;
x=1;
y=2;
[Ntotle1,MMM,number1]=PRPSWP(x);
[Ntotle2,NNN,number2]=gu(y);
N=max(MMM,NNN);
R1=N^(number1);
R2=N^(number2);
result1=Ntotle1*R1;
result2=Ntotle2*R2;
result=(result2/result1)^(1/49);
diary 000result.m
sprintf('%e   %e   %e',result1,result2,result)
diary off