选修MATLAB（CUMT）赵老师作业（部分）

选修MATLAB（CUMT）赵老师作业（部分）

2、用两种不同的方法求解线性方程组4a+ 2b-c +4d= 12 ，3a-b+ 2c-d = 20，12a+ 3b+ 5d= 8，6a+3c=15，。要求给出相应的代码和运行结果。

A=[4 2 -1 4; 3 -1 2 -1;12 3 0 5;6 0 3 0];
B=[12 20 8 15];
X=B*inv(A)

A=[4 2 -1 4;3 -1 2 -1;12 3 0 5;6 0 3 0];
B=[12 20 8 15];
X=B/A

x=-100:0.001:100;
y=(1-(2*(abs(x).^2))+(abs(x).^3)).*(abs(x)>=0&abs(x)<=1)+(4-8*(abs(x))+5*(abs(x)).^2-(abs(x)).^3).*(abs(x)>=1&abs(x)<=2)+0.*(abs(x)>=2);
plot(x,y);

t = 0:pi/50:10*pi;
st0 = sin(t);
ct0 = cos(t);
t1 =3:pi/50:10*pi+3;
st1 = sin(t);
ct1= cos(t);
T=[t;t1];
S=[st0;st1];
C=[ct1;ct1];
mesh(S,C,T);

6、求200以内的质数，以此为任务完成程序设计。要求:(1)不调用MATLAB的求质数的

函数，写出求质数的脚本文件; (2)列出200以内的所有质数; (3)求200以内所有质数的
和。[评分标准:脚本文件5分，列出200以内的所有质数3分，求200以内所有质数的和
2分]

function f= test6_is_prime_08193099(n)
    f = 1;
    n = uint16(n);
    if n<= 1
            return;
        else
            if n<=3
               f = 0;
               return;
            end
        end
        
        if( mod(n,2) ==0 || mod(n,3) ==0 )
           return;                
        end            
        i = 5;
        while(i*i<=n)
             if ( mod(n,i)== 0 || mod(n,i+2) == 0 )
                 return;
             end
            i = i + 6;
    end
    f = 0;   
end

水仙花数

m=100:999;
q=rem(m,10);
w=rem(fix(m/10),10);
e=fix(m/100);
t=find(m==q.^3+w.^3+e.^3);
s=m(t)