matlab 中对变量赋值

syms x y f a b c d e f
f=100*((y-x^2))^2+(1-x)^2 ; %给定的原始多项式
a=diff(f,x); %求x一阶导数
b=diff(a,x); %求x二阶导数
c=diff(f,y); %求y一阶导数
d=diff(d,y); %求y二阶导数 %这一行应当是对c求导
e=diff(a,y); %求x，y二阶偏导数
f=diff(c,x); %求y，x二阶偏导数
S=[a;c]; %求一阶导数矩阵
Q=[b e;f d]; %求二阶导数矩阵
P=inv(Q); %求二阶导数逆矩阵
H=[x;y]; %给定初始的x（0）
K=H-P*S; %求x（1）

%%你的d的表达式出错。
x=2;
y=2;
subs(K)