aufgabe30a.m

来自「作业题 就是说matlab中的常微分方程问题」· M 代码 · 共 21 行

function aufgabe30a

x = colon(0,pi/4,2*pi); 

y = exp(i*x);

z=(-y+y.^2)./(-1/12+2*y./3+5*y./12);
        
fid = fopen('ausgabe30.m','w+');
        
fprintf(fid,'%12.8f  %12.8f\n',x,y);
        
fclose(fid);

x1 =  linspace(0, 2*pi); 

y1 = exp(i*x1);

z1=(-y1+y1.^2)./(-1/12+2*y1./3+5*y1./12);

plot(z1);