编写M文件

（编写M文件）2012-7-6笔记1、求圆的面积
1）%求圆的面积；
title('圆的面积');
%输入圆的半径；
%输出圆的面积；
r=10;
s=pi*r^2
2） %求圆的面积；
title('圆的面积');
%输入圆的半径；
%输出圆的面积；
r=input('r=');
s=pi*r^2
3）%求圆的面积；
title('圆的面积');
%输入圆的半径；
%输出圆的面积；
r=input('r=');
s=pi*r^2;
s=vpa(s,10) 要求有效数位为十位2、求一元二次方程的解
1）顺序语句的（不好）
%求解一元二次方程ax^2+bx+c=0；
%输入系数；
%输出方程的解；
a=input('a=');
b=input('b=');
c=input('c=');
d=b^2-4*a*c;
if d>0|d==0
x1=(-b+sqrt(d))/(2*a);
x2=(-b-sqrt(d))/(2*a);
end
x1
x2
输入 a=1
b=-3
c=-4
运行结果x1 =
4
x2 =
-1
2）顺序语句的程序
%求解一元二次方程ax^2+bx+c=0；
%输入系数；
%输出方程的解；
a=input('a=');
b=input('b=');
c=input('c=');
d=b^2-4*a*c;
if d>0|d==0
x1=(-b+sqrt(d))/(2*a)
x2=(-b-sqrt(d))/(2*a)
else x='方程无实数解'
end
3）条件语句的程序（较好）
%求解一元二次方程ax^2+bx+c=0；
%输入系数；
%输出方程的解；
a=input('a=');
b=input('b=');
c=input('c=');
d=b^2-4*a*c;
if d>0
x1=(-b+sqrt(d))/(2*a)
x2=(-b-sqrt(d))/(2*a)
elseif d==0
x=-b/(2*a)
else x1=-b/(2*a)+sqrt(-d)/(2*a)*i x2=-b/(2*a)-sqrt(-d)/(2*a)*i
end
输入a=1
b=1
c=1
运行结果
x1 =
-0.5000 + 0.8660i
x2 =
-0.5000 - 0.8660i
3、求前n项和
clear
%输入一个整数n；
%输出前n个整数的和；
n=input('n=');
s=0;
for k=1:n
s=s+k;
end
s
4、求n的阶乘
clear
%输入一个整数n；
%求n的阶乘；
n=input('n=')
s=1;
for k=1:n
s=s*k
end
s
5、求1到n的勾股数
clear
%输入一个整数n;
%输出n以内的勾股数；
n=input('n=');
m=1;
for a=3:n
for b=4:n
for c=5:n
if a^2+b^2==c^2&a<b&b<c A=[m a b c]
m=m+1;
end
end
end
end
6、蒙托卡罗法
求π值
1）clear
%输入一个整数n；
%输出pi的近似值;
n=input('n=');
m=0;
for k=1:n
a=rand(2);
x=a(1);
y=a(2);
if x^2+y^2<1
m=m+1;
end
end
pai=4*m/n
2）clear
%输入整数n;
%求π
n=input('n=')
a=rand(n,1);
b=rand(n,1);
c=find(a.^2+b.^2<=1);
d=length(c);
pai=4*d/n
向量的乘法加.a.^2+b.^2<=1
rand(m,n) m*n阶矩阵
rand(n) n*n阶矩阵
求在平面上图像y=x^2与y=12-x交的区域在第一象限x=12于y=9区域上的概率
1）clear
%求在平面上图像y=x^2与y=12-x交的区域在第一象限x=12于y=9区域上的概率；
%输入两个函数；
syms x y;
f=1;
y1=x^2
y2=12-x
ezplot(y1);
hold on;
ezplot(y2);
axis ([0,12,0,12]);
fdx=int(f,x,sqrt(y),12-y)
F=int(fdx,0,9)
2）clear
%求y=x^2 y=12-x y=0 围城的图形的面积；
n=input('n=')
a=rand(n,1)*12;
b=rand(n,1)*9;
d=0;
for k=1:n
if a(k)^2-b(k)>0&12-a(k)-b(k)>0&b(k)>0
d=d+1;
end
end
p=d/n*9*12
中国建模网
/mcm/。