基于matlab完成一个判断年份是否为闰年并且输出日期为本年哪一天的程序

clear all;
clc;
y=input('year=?'); 
m=input('month=?'); 
d=input('day=?'); 
k=0;
if((rem(y,4)==0&&rem(y,100)~=0)||rem(y,400)==0)
   % disp(y);
   fprintf("%d年为闰年，2月份天数为29\n",y);
    k=1;
end
md=[0,31,59,90,120,151,181,212,243,273,304,334,365];
mdr=[0,31,60,91,121,152,182,213,244,274,305,335,366];
if(k)
    date=mdr(m)+d;
else
    date=md(m)+d;
end
fprintf("此天为本年第%d天\n",date);