谁能帮我编写这个matlab程序？

>>x=-5:0.1:5

y=x

z=x.^4-2*x.^2.*y+x.^2+y.^2-2*x+5

[X,Y]=meshgrid(x,y)

Z=X.^4-2*X.^2.*Y+X.^2+Y.^2-2*X+5

subplot(121)

plot3(x,y,z)

subplot(122)

mesh(X,Y,Z)

t = 0:0.01:4*pi

x1 = 10*sin(t)

x2 = 6*abs(sin(2*t))

figure, hold on

plot(t,x1)

plot(t,x2,'--k')

clcclear

syms i a b c

b=[]c=1

for i=1:1000

a=mod(i,10)

if a^3==i

b(c)=i

c=c+1

end

end

b

%求出来的b就是所要找的数。

---