matlab求隐函数方程并用插值法画图

方程：

(225+60*cos(100*t)-180*cos(y))^2+(80+60*sin(100*t)-180*sin(y))^2=120^2
我编的程序是这样：
t=linspace(0,0.02*pi,100);

y=subs(solve('(225+60*cos(100*t)-180*cos(y))^2+(80+60*sin(100*t)-180*sin(y))^2=120^2','y'));
但是，因为t相当于size（1 100）的矩阵，算式里的t和y又有平方，所以这里解出来的y变成两列的矩阵
而且两列元素各不相同，也不是互为相反数
我知道、这样是因为本来用(225+60*cos(100*t)-180*cos(y))^2，想得到的是(225+60*cos(100*t)-180*cos(y))这个矩阵式子里各个元素各自的平方，按说不能只用“^2”，要用".^2"点乘。用“^2”是矩阵相乘的平方、不等于每个元素的各自平方。
但是用“.^2"系统表示有错误。
用a=(225+60*cos(100*t)-180*cos(y)).^2，看了size(a)，是（1,100）；
于是把上面的"120^2"变成了size(1,100)的、元素全为14400的矩阵
但是，用.^2,时，仍然表示有错误
求大神指教O(∩_∩)O哈！

用其他程序解决这个方程也可以
只要还能画出拟合曲线

另，定义xc=225+60*cos（100*t)-180*cos(y),yc=80+60*sin(100*t)-180*sin(y)
画出xc-t，yc-t，xc-yc
求指教！！

ezplot('(225+60*cos(100*t)-180*cos(y))^2+(80+60*sin(100*t)-180*sin(y))^2-120^2')

用上面语句试试画图