如题，请高手指点。
压缩文件为m（new_main,newton,f1,df1）文件和TXT为m文件的代码。df1中需要求雅克比矩阵，程序写得有错误。
请求帮助，不胜感激！ qq:511921847.
输入初始值后运行出错！
>> S0=2.01;
S1=1.71;
S2=1.02;
S3=0.05;
t1=3.5;
t2=13;
t3=23;
>> new_main
??? Undefined function or variable 'S0'.

Error in ==> f1 at 6
f1=x/y*log((S0/S1)*((S1-y)/(S0-y)))+log((S1-y)/(S0-y))-z*t1;

Error in ==> new_ton at 9
x1=x0-f1(x0)/df1(x0);

Error in ==> new_main at 3
[x,n,data]=new_ton(x0);